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Abstract 

We review the present level of knowledge of the hadronic structure of the photon, 
as revealed in interactions involving quarks and gluons "in" the photon. The concept 
of photon structure functions is introduced in the description of deep-inelastic 
scattering, and existing parametrizations of the parton densities in the photon are 
reviewed. We then turn to hard 7p and 77 collisions, where we treat the production 
of jets, heavy quarks, hard (direct) photons, J/ip mesons, and lepton pairs. We also 
comment on issues that go beyond perturbation theory, including recent attempts at a 
comprehensive description of both hard and soft 7^ and 77 interactions. We conclude 
with a list of open problems. 
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1) Introduction 



The photon is the simplest of all bosons. It is the gauge boson of QED, which implies that 
it is massless and structureless (i.e., pointlike). Predictions for 7e interactions can be made 
with impressive accuracy, in some cases (e.g., g — 2 to better than one part in 10 8 . Indeed, 
for many physicists the great quantitative success of QED is one of the strongest arguments 
in favour of Quantum Field Theories in general. 

At first glance it might therefore be surprising, perhaps even a bit embarassing, that 
many reactions involving (quasi-)real photons are much less well understood, both theoret- 
ically and experimentally. However, a moment's reflection will show that this is really not 
astonishing at all. The uncertainty principle tells us that for a short period of time a photon 
can fluctuate into a pair of charged particles. Fluctuations into virtual two-lepton states 
are well understood, and are in fact a crucial ingredient of the quantitative success of QED. 
Fluctuations into quark-antiquark pairs are much more problematic, however. Whenever 
the lifetime of the virtual state exceeds about 10~ 25 sec, corresponding to a characteristic 
energy or momentum scale of about 1 GeV or less, the (virtual) qq pair has sufficient time to 
evolve into a complicated hadronic state that cannot be described by perturbative methods 
only. Even if the lifetime is shorter, i.e. the energy-momentum scale is larger, hard gluon 
emission and related processes complicate the picture substantially. 

The understanding of such virtual hadronic states becomes particularly important when 
they are "kicked on the mass shell" by an interaction of the photon. The most thoroughly 
studied reactions of this type involve interactions of a real and a virtual photon (e.g. in ej 
scattering); of two real photons (77 scattering at e + e~ colliders); and of a real photon with a 
hadron (e.g. 7p scattering). All these reactions allow to probe some aspects of the hadronic 
nature of the photon, and we will discuss them in turn. 

Why is it important to study such reactions? There are at least two answers. First, we 
can use data taken at present colliders to sharpen our predictions for (background) processes 
at future colliders. We will see that in some cases this may contribute significantly to the 
assessment of the potential of future colliders; this is true for high-energy linear e + e~ colliders 
that are now being discussed, and especially for the so-called 77 colliders. Secondly, while 
the hadronic structure of the photon certainly has nonperturbative aspects, we expect the 
photon to be a simpler system than any real hadron, like the proton. After all, no matter 
how complicated it is, the hadronic structure of the photon certainly originated from the 'jqq 
vertex; if this vertex vanished, the photon would not have any "hadronic nature" to speak of. 
No such (in principle) simple starting point can be defined for protons, constituent quarks 
themselves being highly complicated objects. One can therefore expect photonic reactions to 
be particularly well suited to study both perturbative and nonperturbative aspects of strong 
interactions, and everything in between. We should mention right away that the present 
theoretical description of, say, 7£> scattering is, if anything, even more complex than that of 
pp scattering; however, we feel that this does not refute our argument that 7p collisions are 
simpler at some fundamental level. At the very least reactions probing the hadronic nature of 
the photon give us an additional handle on most aspects of perturbative and nonperturbative 
QCD that are of relevance for collider physics. 

The remainder of this article is organized as follows. We start in Sec. 2 with a discussion 
of the scattering of a highly virtual (probing) photon off an almost real (target) photon, 7*7 
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or cy scattering. These were the first photonic reactions for which predictions were made in 
the framework of the quark parton model (QPM) and within QCD [[|. e7 scattering was 
also among the first of the "hard" photonic reactions, which can at least partly be described 
by perturbation theory, to be studied experimentally ||. Finally, the conceptual simplicity 
of these processes makes them ideal for introducing the notion of photon structure functions 
and the parton content of the photon. 

The commissioning of the ep collider HERA has opened a new era in the experimental 
study of hard 7p collisions, increasing the available jp centre-of-mass energy by about a 
factor of ten compared to previous experiments. In Sec. 3 we attempt to cover both the 
recent intense theoretical activity || that has been triggered by this prospect, and the 
relevant experimental data. 

Lately there has also been much progress in the field of hard scattering of two quasi-real 
photons. An important milestone was set by the AMY collaboration at the TRISTAN 
collider, who for the first time analyzed their data using an essentially complete leading order 
QCD event generator. Previous generators had omitted important classes of diagrams, and 
had therefore not been able to reproduce data taken at the older PEP and PETRA colliders. 
Recent developments in this area are summarized in Sec. 4. 

The emphasis in Sees. 2 to 4 is on hard processes, which are amenable to a perturbative 
treatment. In Sec. 5 we loosen this restriction and comment on issues that go beyond 
standard perturbative QCD. In particular, we discuss ■yp and 77 cross-sections, and related 
quantities, in the minijet model [0. This model has been introduced to describe purely 
hadronic (pp, pp) reactions; some modifications are necessary || before it can be used for 
7P and 77 scattering. Finally, Sec. 6 contains some concluding remarks, as well as a list of 
open problems and experimental challenges. 

This field is still very much in flux. Hardly a week goes by without a new experimental 
study of some reaction that probes the hadronic structure of the photon, a new or refined 
calculation of a relevant cross-section for given photon structure, and/or a new model for 
or parametrization of this structure. We nevertheless think it worthwhile to summarize the 
present status, partly to celebrate what has already been achieved, but mostly to highlight 
open questions and how to address them. We hope that this review, which updates ref.[(|, 
will be of some use both for those who have already worked on some aspects of the hadronic 
nature of the photon, and for those who contemplate starting such work. 



2) Photon Structure Functions 



In this section we introduce the concept of photon structure functions. To this end we first 
discuss in Sec. 2a deep-inelastic scattering, which is theoretically very clean, being fully 
inclusive; it is thus well suited to serve as the defining process for photon structure functions 
and the parton content of the photon. For reasons of space we have to be rather brief here; we 



refer the reader to ref . [lCj for a more pedagogical introduction to photon structure functions. 
In Sec. 2b we then describe existing parametrizations of the photonic parton distribution 
functions. 
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2a) Deep— inelastic e^j Scattering 

Studies of deep-inelastic electron nucleon scattering, 

eN -> eX, (1) 

where X is any hadronic system and the squared four momentum transfer Q 2 = —q 2 > 1 
GeV 2 , have contributed much to our understanding of strong interactions. The pioneering 



SLAC experiments [TT| played a key role in the development of both the quark parton model 
(QPM) fll^f , and, due to the observation of scaling violations and the concept of asymptotic 
freedom ||13|| , of QCD. Since then our understanding of QCD in general and the structure 
of the proton in particular has improved greatly. The most recent progress in this area has 



come from the analysis of data taken at the ep collider HERA; see ref.|L4j for a study of the 
impact of these data on the determination of the parton content of the proton. 

Formally deep-inelastic cy scattering is quite similar to ep scattering. The basic kinemat- 
ics is explained in Fig. 1. The differential cross section can be written in terms of the scaling 
variables x = Q 2 /(2p ■ q) and y = Q 2 /(sx), where a/s is the total available centre-of-mass 
(cms) energy: 



rf 2 cr(e7 — > eX) 
dxdy 



2 ™eV 



Q 4 



F2(x,Q 2 ) 



x 



Q 2 )} ; 



(2) 



this expression is completely analogous to the equation defining the protonic structure func- 
tions i<2 and Fl in terms of the differential cross-section for ep scattering via the exchange of 
a virtual photon.Q The special significance || of ey scattering lies in the fact that, while (at 
present) the x— dependence of the nucleonic structure functions can only be parametrized 
from data, the structure functions appearing in eq.(^) can be computed in the QPM from 
the diagram shown in Fig. 2a: 



_p7,QPM ( 



x,Q 2 



Qa 



Ti 



i \ -> 4 



x 2 + (1-x) 2 



W 2 

log — - + 8x (1 — x) 



- 1 



m„ 



(3) 



q K " u q 

where we have introduced the squared cms energy of the 7*7 system 



W 2 = Q 2 ( - - 1 ) . 



(4) 



The sum in eq.(|3|) runs over all quark flavours, and e q is the electric charge of quark q in 
units of the proton charge. An experimental test of this equation was thought to not only 
allow to confirm the existence of pointlike quarks, but also to measure their charges through 
the factor; both were topics of interest in the early 1970's when the study of ej scattering 
was first proposed 0. 

Unfortunately, eq.(^) depends on the quark masses m q . If this ansatz is to describe data 
Il| even approximately, one has to use constituent quark masses of a few hundred MeV here; 



constituent quarks are not very well defined in field theory. Moreover, we now know that 

*The Z exchange contribution to eq.(||) is negligible for Q 2 values that can be achieved in the foreseeable 
future. 
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QPM predictions can be modified substantially by QCD effects. In case of proton structure 
functions these lead to, among other things, scaling violations and a nonzero Fl. In case of 
e7 scattering, QCD corrections are described by the kind of diagrams shown in Figs. 2b, c. 
Diagrams of the type 2b leave the flavour structure unchanged and are therefore part of 
the (flavour) nonsinglet contribution to F^, while diagrams with several disconnected quark 
lines, as in Fig. 2c, contribute to the (flavour) singlet part of F]. 

The interest in photon structure functions received a boost in 1977, when Witten showed 
[[J that such diagrams can be computed exactly, at least in the so-called "asymptotic" limit 
of infinite Q 2 . Including next-to-leading order (NLO) corrections [16|, the result can be 
written as 

1 



F 2 7 ' asymp (x,Q 2 ) =a e 



-a(x) + b(x) 



(5) 



i; 



MQ 2 

where a and b are calculable functions of x. The absolute normalization of this "aymptotic 
solution is therefore given uniquely by a s (Q 2 ), i.e. by the value of the QCD scale parameter 
Aqcd- It was therefore hoped that eq.@ might be exploited for a very precise measurement 
of Aqcd- 

Unfortunately this no longer appears feasible. One problem is that, in order to derive 



eq.(p3), one has to neglect terms of the form ( S [S. 9 J , ] , where Ql is some input scale (see 

\<Xa{Qo)J 

below). Neglecting such terms is formally justified if a s (Q 2 ) <C a s (Q 2 ] ) and P is positive. 
Unfortunately the first inequality is usually not satisfied at experimentally accessible val- 
ues of Q 2 , assuming Qq is chosen in the region of applicability of perturbative QCD, i.e. 
«s(Qo)/ 7r ^ 1- Worse yet, P can be zero or even negative! In this case ignoring such terms 
is obviously a bad approximation. Indeed, one finds that eq.(|5]) contains divergencies as 
^0|1: 

a(x) ~ a;" ' 59 , b(x) ~ aT 1 . (6) 

The coefficient of the 1/x pole in b is negative; eq.([5]) therefore predicts negative counting 
rates at small x. Notice that the divergence is worse in the NLO contribution b than in the 
LO term a. It can be shown [|17]] that this trend continues in yet higher orders, i.e. the 
"asymptotic" prediction for F^ rapidly becomes more and more divergent for x — > as more 
higher order corrections are included: F% ~ x -4 ' 3 in 3rd order, ~ x -25.6 - m 4^ orc j er) anc [ 
so on. Clearly the "asymptotic" solution is not a very useful concept, having a violently 
divergent perturbative expansion.[] 

The worst divergencies in F^' asymp occur in the singlet sector, i.e. originate from diagrams 
of the type shown in Fig. 2c. It has been speculated jL8) that this hints at a nonperturba- 
tive solution. For example, if the invariant mass of the lowest qq pair in Fig. 2c is small, 
they might form a bound state. Traditionally, however, nonperturbative contributions have 
been estimated using the vector dominance model (VDM) ||19|| , from the diagrams shown 
in Fig. 2d: The target photon undergoes a transition into a nearly on-shell vector meson 
(p, uj, 0,...), so that e7 scattering is "really" ep, eu, . . . scattering, which should look 
qualitatively like ep scattering. In particular, the contribution of Fig. 2d should itself be 



^Notice that in the same "asymptotic" limit, nucleonic structure functions collapse to a 5— function at 
x — 0. While this is formally correct for infinite Q 2 , it gives obviously a poor description of the true proton 
structure at any finite Q 2 . 
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well-behaved, i.e. non-singular; it cannot cancel the divergencies of the "asymptotic" solu- 
tion! 

Indeed, the existence of the contribution shown in Fig. 2d demonstrates that we cannot 
hope to compute F^{x, Q 2 ) from perturbation theory alone. Moreover, even if we assume 
that the VDM correctly describes all nonperturbative contributions to F%, it seems essen- 
tially impossible to estimate them without making further assumptions. The problem is that 
the vector mesons p, u, <fi, . . . are much too short-lived to allow an independent measurement 
of their parton distribution functions.^] It therefore seems to us that the only meaningful 
approach is that suggested by Gliick and Reya pOf . That is, one formally sums the contri- 
butions from Figs. 2a-d into the single diagram of Fig. 2e, where we have introduced quark 
densities in the photon q](x,Q 2 ) such that (in LO) 



F2(x,Q 2 ) 



2xJ2elq7(x,Q 2 ) 



(7) 



where the sum runs over flavours, e qi is the electric charge of quark q^ in units of the proton 
charge, and the factor of 2 takes care or antiquarks. This is merely a definition. In the 
approach of ref . [^(J one does not attempt to compute the absolute size of the quark densities 
inside the photon. Rather, one introduces input distribution functions ql Q (x) = q](x,Ql) 
at some scale Qq. This scale is in principle arbitrary, as long as a s {Ql) is sufficiently small 
to allow for a meaningful perturbative expansion. In practice, is usually chosen as the 
smallest value where this criterion is assumed to be satisfied. We will come back to this 
point in the following subsection. 

Given these input distributions, the photonic parton densities, and thus F%, at different 
values of Q 2 can be computed using the inhomogeneous evolution equations. In LO, they 
read §11: 



dq&s( x >Q 2 ) 
dlogQ 2 

<fi?i{x,Q 2 ) 
dlogQ 2 

dG^{x,Q 2 ) 
dlogQ 2 



aem -kl s (x) + " s(g2 



2tt 



2tt 



a 



em / 7 



27T 

a s (Q 2 
2tt 



kl(x) + 



a s (Q 2 
2tt 



P» q ® (x, Q 2 ) + (P^ G ® G 1 ) (x, Q 



P Gq ® (a, Q 2 ) + (P^ G ® G 1 ) (x, Q 2 ) 



where we have used the notation 



[P®q)(x,Q 2 )= [ 1 ^-P(y)q(^Q 
Jx y v 



The Ptj are the usual (LO) j 



i splitting functions 



(8b) 



ic) 



(9) 



The inhomogeneous terms kj 



describe 7 — > qq splitting, i.e. the diagram of Fig. 2a; for one quark flavour, one has 



x 



3e 



x 



+ (1 



x 



(10) 



•''Of course, one can always define F2 to be the sum of the (divergent) "asymptotic" solution and a 
(divergent) "nonperturbative contribution" . However, nothing is gained by this as long as one cannot at 
least estimate the nonperturbative contribution. Our argument shows that the VDM is of no help here. 

§ Assuming that such distributions can even be defined for resonance states; e.g., should a p be treated as 
a qq resonance with two valence quarks, or a tttt resonance with four valence quarks? 
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The k] of eqs.(H) follow from eq.(]rD|) by taking appropriate sums or differences of quark 
flavours. Eq. ( |8aD describes the evolution of the nonsinglet distributions (differences of quark 
densities), i.e. resums only diagrams of the type shown in Fig. 2b, while eqs. (|§Tj,|c|) describe 
the evolution of the singlet sector (£ 7 = + qj), which includes diagrams of the kind 

shown in Fig. 2c. Notice that this necessitates the introduction of a gluon density inside the 
photon G 7 (x, Q 2 ), with its corresponding input distribution Gq(x) = G 7 (x, Ql). 

It is crucial to note that, given non-singular input distributions, the solutions of eqs.(||) 
will also remain [EO] well-behaved at all finite values of Q 2 . This is true both in LO and 



in NLO [ 23[ . On the other hand, by introducing a priori unknown input distributions, one 
clearly abandons the hope to make an absolute prediction of F%(x, Q 2 ) in terms of Aqcd 
alone. The solutions of eqs.([5D still show an approximately linear growth with logQ 2 ; in this 
sense eq.([|) remains approximately correct, but the functions a and b now do depend weakly 
on Q 2 (approximately like log log Q 2 ), and the x— dependence of b is not computable. In fact, 
only the nonleading Q 2 dependence (contained in the functions a and b), which corresponds 
formally to the scaling violations in the case of ep scattering^], can in this approach be 
used to determine Aqcd! a change of the normalization l/a s multiplying the first term can 
always be compensated by adding a constant to the input distributions. Notice that no 
momentum sum rule applies for the parton densities in the photon as defined here. The 
reason is that these densities are all of first order in the fine structure constant a em . Even 
a relatively large change in these densities can therefore always be compensated by a small 
change of the 0(a^ m ) term in the decomposition of the physical photon, which is simply 
the "bare" photon [with distribution function 5(1 — x)\. Formally this would manifest itself 
by the addition of 0(al m ) contributions to the inhomogeneous terms in eqs.(||), which are 
numerically negligible. [J] 

Before discussing our present knowledge of and parametrizations for the parton densities 
in the photon, we briefly address a few issues related to the calculation of As mentioned 
above, eqs.(^),(||) have been extended to NLO quite early, although a mistake in the two-loop 



7 — ► G splitting function was found |2(| only fairly recently. A full NLO treatment of massive 



quarks is now also available [27] for both and F2- A first treatment of small— a; effects 
in the photon structure functions, i.e. logl/x resummation and parton recombination, has 
been presented in ref . |28| ; however, the predicted steep increase of F% at small x has not been 
observed experimentally [Q. Finally, nonperturbative contributions to are expected to 
be greatly suppressed if the target photon is also far off-shell. One can therefore derive 
unambiguous QCD predictions f3(J in the region Q 2 ^> P 2 ^> A 2 , where the first strong 



inequality has been imposed to allow for a meaningful definition of structure functions| 



^Recall that in case of scattering, the QPM also predicts a logarithmic growth of F. 7 with Q 2 , see 
eq.©. " ; 

lift is sometimes argued that one can still use the "asymptotic" NLO prediction for to determine 
Aqcd, if one sticks to the region of large x where the influence of the x — > pole is supposed to be weak. 
However, this rests on the assumption that the nonperturbative contributions to F% are small at large x, and 
the second assumption that the terms that regularize the pole pj| vanish as x — > 1. The only way to test 
these assumptions seems to be to compare the value of Aqcd extracted in this way with other measurements 
p5[ of a s . In our opinion this shows that a measurement of at fixed Q 2 cannot be used to determine 
Aqcd unambiguously. Recall also that the perturbative expansion of the "asymptotic" prediction in yet 
higher orders in QCD is extremely problematic. 

**If Q 2 ~ P 2 , the use of fixed-order perturbation theory is more appropriate, since it includes terms 
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However, it has recently been pointed out [5I| that nonperturbative effects might survive 
longer than previously expected; an unambiguous prediction would then only be possible for 
very large P 2 , and even larger Q 2 , where the cross-section is very small. 



2b) Parametrizations of Photonic Parton Densities 

As discussed in the previous subsection, the Q 2 evolution of the photonic parton densities 
<P(x, Q 2 ) = (qj, G 7 )(x, Q 2 ) is uniquely determined by perturbative QCD, eqs.(|8p and their 
NLO extension. However, given the flaws of the "asymptotic" prediction for (p, one has to 
specify input distributions aft at a fixed Q 2 — Qq- The situation is therefore in principle quite 
similar to the case of nucleonic parton densities. In practice, however, the determination of 
these input distributions is much more difficult in case of the photon, for a variety of reasons. 

To begin with, no momentum sum rule applies for as discussed above. This means 
that it will be difficult to derive reliable information on Gq from measurements of F% alone: 
in LO, the gluon density only enters via the (subleading) Q 2 evolution in i 7 ^. In contrast, 
analyses of eN scattering revealed very early that gluons must carry about 50% of the 
proton's momentum, thereby fixing the overall scale of the gluon density in the nucleon. We 
will see below that existing parametrizations for G 1 still differ by sizable factors over the 
entire x range unless Q 2 is very large. 

Secondly, so far deep-inelastic scattering could only be studied at e + e~ colliders, 
where the target photon is itself radiated off one of the incoming leptons. Its spectrum is 



given by the well-known Weizsacker- Williams function |]32| 



Op 



2tcz 



2E 2 (l-z) 2 (l-cos9 n 



Ldz) = ^{ ! + (!-*)' log - ^ ~> \\ — -2(1-,) , (11) 



m 2 z 2 



where E is the electron beam energy and zE the energy of the target photon. In eq. (|TT|) 
we have assumed that there is no experimentally imposed lower bound on the virtuality 
p 2 = —P 2 of the target photon, so that P^ in = m 2 z 2 /(l — z); however, we have introduced 
an upper limit # max on the scattering angle of the electron emitting the target photon, in 
order to allow for antitagging. 

Eq.(|ll]) implies that the cross section from the measurement of which F£ is to be de- 
termined is of order da/dQ 2 ~ a^/ (vrQ 4 ) log (E/m e ), see eq.(@). [Recall that Fl is itself 



0{a em ).] The event rate is therefore quite small; the most recent measurements [15|, |29|, |33 
typically have around 1,000 events at Q 2 ~ 5 GeV 2 , and the statistics rapidly gets worse at 
higher Q 2 . This is to be compared with millions of events in, for example, deep-inelastic vN 
scattering [|4[ 



Another problem is that the e ± emitting the target photon is usually not detected, since 
it emerges at too small an angle. This means that the energy of the target photon, and 
hence the Bjorken variable x, can only be determined from the hadronic system. All existing 
analyses try to determine x from the invariant mass W, using eq.(|]). This is problematic, 
since at least some of the produced hadrons usually also escape undetected, e.g. in the beam 
pipes; the measured value of W, usually denoted by W v i S , is therefore generally smaller than 
the true value. (It can exceed the true W due to the finite energy resolution of real-world 

oc (P 2 /Q 2 ) n , which are not included in the usual structure functions. 
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detectors.) One has to correct for this by "unfolding" the measured W V i S distribution to 
arrive at the true W (or x) distribution. However, in order to do this one has to model 
the hadronic system X. In other words, a "measurement" of F£ only seems possible if we 
already know many details about multi-hadron production in 7*7 scattering! 

In practice the situation is not quite as bad, since one can check the assumptions made by 
comparing various distributions predicted by the Monte Carlo with real data. An iterative 
procedure can then be used to arrive at a model for X that should allow to do the unfolding 
reliably. However, one should be aware that this might lead to large uncertainties at the 
boundaries of the accessible range of x values. The reason is that in the event selection 
one usually imposes upper and lower cuts on W v i S , corresponding to small and large x, 
respectively. Relatively minor changes of the model for X can therefore significantly change 
the predicted efficiency for accepting events with x close to one of the kinematic boundaries. 
The region of large x (small W) is in any case plagued by higher twist uncertainties (e.g., 
production of resonances); however, the region of small x is very interesting, since here 
F2 is dominated by sea quarks whose density is closely related to the gluon distribution. 
It has been shown explicitly [|35|] that different ansatze for X can lead to quite different 



"measurements" of F% at small x. 



Most recent analyses use the "FKP model" of refs. |36[ as starting point of the unfolding 
procedure. Here F% is split into "soft" and "hard" components, depending on the virtuality 
of the quark exchanged in the t— or u— channel, see figs. 2a-c. The cut-off parameter to 
separating these two regions is to be fitted from the data. Contributions with \t\ < to are 
modelled using VDM ideas; in practice this means that a scaling (Q 2 — independent) ansatz 
is used, fitted from low— Q 2 data fll5| . Contributions with \t\ > t can be computed from an 



evolution equation in log|t|, with the boundary condition that this "pointlike" part vanish 
at t = t . This procedure is formally equivalent to the one suggested by Gliick and Reya 
[BO], if we identify Ql = t and impose the upper limit \t\ max = Q 2 for the t— evolution. 



However, as presently used [37, 29, 33 1, the unfolding procedure has several weaknesses: 



It uses a parametrization J37J] of the "pointlike" part of F% which includes some terms 
which are of next-to-leading order in a formal operator product expansion; however, 
not all such terms are included. In particular, it uses the actual kinematical maximum 
for the t— evolution, |t| max — Q 2 /x; as a result, the predicted tends to be larger 
at small and median x than what one would get from the usual Q 2 — evolution.^ On 
the other hand, this expression ignores sea quarks, i.e. g — > qq splitting; it therefore 
underestimates the true result for very small x. 

The Q 2 (or t) evolution of the "hadronic part" is ignored; at high Q 2 , this over- 
estimates the soft contribution to F^ at median and large x, and underestimates it at 
small x. 

While this procedure treats the effect of gluon radiation on the shape of F2 (approxi- 
mately) correctly, it does not include any parton showering in the MC event generator. 
Rather, the generator produces qq pairs whose px distribution in the 7*7 cms is either 
exponential (for the "hadronic part") or follows the QPM prediction (for the "point- 
like part"). This qq pair is then fed into a string-based fragmentation program. While 



*Very roughly, one replaces the overall growth oc logQ 2 by log (Q 2 /x). 
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string fragmentation can mimick the effects of parton showering to some extent, it is 
not able to produce additional jets. It is therefore not surprising that a recent analysis 
|33| found larger jet rates in ej scattering than the MC predicted. 



Some of these points have also been raised in ref.pHJ. We would like to emphasize that none 



of these problems is intrinsic to the formalism of ref. |36|]. Moreover, the input model used 
for the unfolding always turned out to be consistent within errors with the extracted if 
to is chosen in the vicinity of 0.5 GeV 2 . 

The point of this lengthy discussion is not to criticize our experimental colleagues. Rather, 
we hope that it might serve as a starting point for further work by people familiar with the 
design of MC event generators, which we are not. However, this discussion shows that present 
data on F% have to be taken with a grain of salt. In particular, the shortcomings listed 
above are usually not addressed in the estimates of the systematic error due to the unfolding 



procedure; this error only includes things like the choice of binning [29]. This might help 



to explain the apparent discrepancy between different data sets []15 |. Fortunately, new ideas 



for improved unfolding algorithms p9] are now under investigation; this should facilitate the 
measurement of K7 at small x, especially at high energy (LEP2). 

In spite of this, measurements of F% probably still provide the most reliable constraints 
on the input distributions <jo(x); they are certainly the only data that have been taken into 
account when constructing existing parametrizations of cp(x,Q 2 ). In the remainder of this 
subsection we will briefly describe these parametrizations. 



The simplest and oldest parametrizations p(| |41| are based on the "asymptotic" LO 
prediction || [4^]. While the theoretical basis for this prediction is weak, as described 
in the previous subsection, their simplicity makes these parametrizations useful for first 
estimates of reaction rates, as long as one stays away from very small x. Refs. |]4"0"|, ffl| only 
give parametrizations for Nf = 4 active flavours. More recently, Gordon and Storrow |43| 
provided different, more accurate fits for Nf = 3, 4 and 5. 

All other parametrizations involve some amount of data fitting. However, due to the 
rather large experimental errors of data on K7, additional assumptions always had to be 
made. In particular, quark and antiquark distributions (of the same flavour) are always as- 
sumed to be identical, which guarantees that the photon carries no flavour. This assumption 
is eminently reasonable. The 'jqq vertex treats quarks and antiquarks symmetrically, and 
we do not know of any effect that could destroy this symmetry. 



The DG parametrization p4| was the first to start from input distributions. At the time 
of this fit, only a single measurement of F% at fixed Q 2 ~ 5.3 GeV 2 existed. In order to 
determine three, in principle independent, input distributions (for nonsinglet and singlet 
quarks as well as gluons), two assumptions were made: All input quark densities were 
assumed to be proportional to the squared quark charges, i.e. u 1 = 4cf y = 4s 7 at Qq = 1 
GeV 2 ; and the gluon input was generated purely radiatively, i.e. 



2 

Po 

where po = 9 is the coefficient of the 1-loop QCD f3 function. This parametrization only 
exists in LO. Moreover, it treats flavour thresholds by introducing three independent sets of 
distributions for Nf = 3, 4, 5, so that the transition across a threshold is not automatically 
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smooth. Nevertheless it continues to describe data reasonably well, although a combined fit 
would probably give a higher \ 2 than for more recent parametrizations. 

The LAC parametrizations |{45| are based on a much larger data set. The main point 
of these fits was to demonstrate that data on F% constrain G 1 very poorly. In particular, 
they allow a very hard gluon, with xG 1 having a maximum at x ~ 0.9 (LAC3), as well as 
very soft gluon distributiuons, with xG 1 rising very steeply at low x (LAC1, LAC2). The 
LAC parametrizations only exist for Nf = 4 massless flavours and in LO. No assumptions 
about the relative sizes of the four input quark densities were made in the fit. LAC3 has 
been clearly excluded by data on jet production in ep scattering (Sec. 3a) as well as in real 
77 scattering (Sec. 4a); the experimental status of LAC1,2 is less clear. 

The recent WHIT parametrizations f46 | follow a similar philosophy as LAC, at least 



regarding the gluon input; however, their choices for Gq are much less extreme. In the 
WHIT 1,2, 3 parametrizations, gluons carry about half as much of the photon's momentum 
as quarks do (at the input scale Ql = 4 GeV 2 ), while in WHIT4,5,6 gluons and quarks carry 
about the same momentum fraction. These two groups of fits also have slightly different 
valence quark densities.]] The input gluon density is assumed to have the simple shape 
xGl oc (1 - x) c \ with c g = 3 for WHIT1,4; c g = 9 for WHIT2,5; and c g = 15 for WHIT3,6. 
Recall that the normalization is adjusted such that / xGl(x)dx is constant within each 
group of parametrizations; a larger c g therefore also means a larger G 1 at small x. The 
input distributions for the sea quarks are computed from the cross section for 7*g — > qq, 
regularized by light quark masses m = 0.5 GeV. The WHIT parametrizations only exist in 
LO, but great care has been taken to treat the (x— dependent) charm threshold correctly. 
This is much more important here than for nucleonic parton densities, since the photon very 
rapidly develops an "intrinsic charm" component from 7 — > cc splitting. 

The GRV parametrization f47| is the first NLO fit of q 1 ; a LO version is also available. 
This parametrization is based on the same "dynamical" philosophy as the earlier fits of 
protonic [48|] and pionic [49] parton densities by the same authors. The idea is to start from 



a very simple input at a very low Qq (0.25 GeV 2 in LO, 0.3 GeV 2 in NLO); this scale is 
assumed to be the same for p, n and 7 targets. The observed, more complex structure is then 
generated dynamically by the evolution equations. In case of the proton, only valence quarks 
were originally assumed to be present at scale Q%. While the gluon density does evolve fast 
enough to carry approximately half the proton's momentum at Q 2 of a few GeV 2 , it was 
found to be too soft in shape. The ansatz therefore had to be modified [SU| by introducing 
"valence-like" gluon and even sea quark densities already at the input scale, thereby giving 
up much of the original simplicity of the idea. Their pionic input distributions also include 
a "valence-like" gluon density, which is in fact strictly proportional to the valence quark 
density, but no sea quarks at scale Qq. 
In the photonic case, 

-7 / \ 47TQ; em _^ . . . x 

%,GRY\ X ) — K p, <7o,GRVl X J- \ 16 ) 
J o 



^The definition of "valence" and "sea" quarks used here differs from the more common "nonsinglet" and 
"singlet" distributions. The Q 2 evolution of the valence density is independent of G 7 , i.e. obeys eq.(8a), 
while a non-zero sea quark density is produced only through g — > qq splitting. In the absence of mass 
effects, the valence distributions are proportional to the squared quark charges, while the sea distributions 
are independent of the quark charge. 
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This input is motivated from VDM ideas, where / 2 determines the 7 — > p transition proba- 
bility (fp/Air ~ 2.2), and k has been introduced to describe contributions from heavier vector 
mesons (u, <p, . . .). In fact, k is the only free parameter in this ansatz; it was determined to 
be k = 2 (1.6) in LO (NLO). However, it should be clear that eg. (|I~3|) is an assumption that 



has to be tested experimentally; in particular, it is not obvious that the parton densities 
in a pion resemble those in a vector meson, or that QCD is applicable at such small input 
scales.^] On the positive side, the GRV parametrization ensures a smooth onset of the charm 
density, using an x— independent threshold. Moreover, care has been taken to split the NLO 
parametrizations for q 1 into LO and NLO pieces; only the former should be multiplied with 
NLO pieces of hard cross sections.| 



The GS parametrizations p3| were developed shortly after GRV, but follow a quite dif- 



ferent strategy. Problems with low input scales |53| are avoided by choosing Qq = 5.3 GeV 2 . 



This is certainly in the perturbative region, but necessitates a rather complicated ansatz for 
the input distributions: 

CgsO) = K^r^fa.Qo) + <?qpm(^<2o)- ( 14 ) 
Ip 

The free parameters in the fit are the momentum fractions carried by gluons and sea-quarks 
in the pion, the parameter k, and the light quark masses. In the GS2 parametrization, 



Gq is assumed to come entirely from the first term in eq. (|I~4"D , while in GS1 the second 



term also contributes via radiation, see eq. (|T2"D . The fit gives m u = m& = 0.29 GeV and 



k = 1.96. These values are not unreasonable. However, the ansatz (|14]) might be somewhat 
suspect: It holds in the perturbative region of Q 2 , but its form is not invariant under the 
evolution equations. On the other hand, for practical purposes it includes sufficiently many 
free parameters to allow a decent description of data on F^. This parametrization is now 



being updated |5^| . The new version uses a slightly reduced input scale Ql = 3 GeV 2 , and 



for the first time includes data on jet production in two-photon collisions (see Sec. 4a) in 
the fit; unfortunately this still does not allow to pin down G 1 with any precision. 

The AGF parametrization |f)5H is (in its "standard" form) quite similar to GRV. In par- 
ticular, they also assume that at a low input scale Ql = 0.25 GeV 2 the photonic parton 
densities are described by the VDM. There are some differences, however. First, AGF point 
out that in NLO the input densities are scheme-dependent, if physical quantities like 
are to be scheme-independent. GRV use their ansatz (p~3|) in their DIS 7 scheme, since it 
is perturbatively more stable than the more commonly used MS scheme. AGF point out 
that this treatment includes certain process-dependent terms in g 7 ; they therefore prefer to 
use the MS scheme, and define their input distributions to be the difference of a (regular) 
"VDM" term and a process-independent term containing a log(l — x) divergence, so that 
F^{x, Qq) is well-behaved in the limit x — > 1. Secondly, they include p — u — <fi interference 

•fit has been shown previously (5l[] that an ansatz like eq.([l3|) cannot be brought into agreement with the 
data if one insists on Qq > 1 GeV 2 . One way out [H would be to multiply the resulting F% with something 
like Q 2 /(Q 2 + /i 2 ), but this goes beyond the leading-twist partonic contributions. 

§This distinction is more important in the photonic case, since possible log(l — x) divergencies are not 
always regularized by parton densities falling like a power of 1 — x. However, even in the photonic case this 
problem is greatly ameliorated if one uses the "DIS 7 " scheme introduced (52) by the same authors. 
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when specifying the "VDM" input, so that 



<AGfO) = Ka e 



4 2 2 

I 2 2 

g W valence( X ) Qo) + g^seaO 3 ^ Qo) 



Gl AGF (x) = Ka em ^G*(x,Q 2 ), (15) 



where the pion distribution functions are taken from ref . [ 56| . (Recall that in NLO the in- 
put quark densities have to be modified by a subtraction term.) Perhaps the for practical 
purposes most important difference from the GRV parametrization is that the coefficient K 
is left free, i.e. separate fits are provided for the "anomalous" (or "pointlike") and "non- 
perturbative" contributions to q 7 , allowing the user to specify the absolute normalization 
(although not the shape) of the latter. 



Finally, two of the SaS parametrizations |38| are based on a similar philosophy as the 



GRV and AGF parametrizations, by assuming that at a low Qo — 0.6 GeV the perturba- 
tive component vanishes (SaSl). However, while the normalization of the nonperturbative 
contribution is taken from the VDM (including p, u and cf) contributions with fixed normal- 
ization), the shapes of the quark and gluon distributions are fitted from data. Although the 
SaS parametrizations are available in LO only, the authors attempt to estimate the scheme 
dependence (formally an NLO effect) providing a parametrization (SaSIM) where the non 
leading-log part of the QPM prediction for F2 has been added to eq.(|7|), while SaSID is 
based on eq.(0) alone. There are also two parametrizations with Q = 2 GeV; however, in 
this case the normalization of the fitted "soft" contribution had to be left free, and its shape 
is much harder than one what expects from hadronic parton densities. Again two sets of 
parametrizations are available, using different schemes (SaS2D, SaS2M). 

The SaSl sets preferred by the authors are quite similar to AGF; the real significance 
of ref. [ I is that it carefully describes the properties of the hadronic state X for both the 



hadronic and "anomalous" contributions, as needed for a full event characterization. We 
will come back to this aspect of their work in Sec. 5. 

In Fig. 3 we compare various LO parametrizations of F2 at Q 2 = 15 GeV 2 with recent 



data taken by the OPAL p9| and TOPAZ [|33| collaborations; present data are not able to 
distinguish between LO and NLO fits. In order to allow for a meaningful comparison, we have 
added a charm contribution to the OPAL data, as estimated from the QPM; this contribution 
had been subtracted in their analysis. We have used the DG and GRV parametrizations with 
Nf = 3 flavours, since their parametrizations of c 7 are meant to be used only if \ogQ 2 /m 2 c ^ 
1; the charm contribution has again been estimated from the QPM.|^| As discussed earlier, 
WHIT provides a parametrization of c 7 that includes the correct kinematical threshold, while 
LAC treat the charm as massless at all Q 2 . 

We see that most parametrizations give quite similar results for F2 over most of the 
relevant a;— range; the exception is LAC1, which exceeds the other parametrizations both at 
large and at very small x. It should be noted that the data points represent averages over 
the respective x bins; the lowest bin starts at x = 0.006 (0.02) for the OPAL (TOPAZ) data. 

^We have ignored the small contribution [p?] from j*g — > cc in this figure. 
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The first OPAL point is therefore in conflict J35J with the LAC1 prediction. Unfortunately 
there is also some discrepancy between the TOPAZ and OPAL data at low x. As discussed 
above, one is sensitive to the unfolding procedure here; for this reason, WHIT chose not to 
use these (and similar) points in their fit. (The other fits predate the data shown in Fig. 3.) 
This ambiguity in present low— x data is to be regretted, since in principle these data have 
the potential to discriminate between different ansatze for Gq. This can most clearly be seen 
by comparing the curves for WHIT4 (long dashed) and WHIT6 (long-short dashed), which 
have the same valence quark input, and even the same / xG^dx: WHIT4 has a harder gluon 
input distribution, and therefore predicts a larger at x ~ 0.1; WHIT6 has many more 
soft gluons, and therefore a very rapid increase of F^ for x < 0.05, not unlike LAC1. Finally, 
we should mention that the GS, AGF and SaS parametrizations also reproduce these data 
quite well. 

Discriminating between these parametrizations would be much easier if one could measure 
the gluon density directly. This is demonstrated in Fig. 4a, where we show results for xG 1 
at the same value of Q 2 ; we have chosen the same LO parametrizations as in Fig. 3, and 
included the LAC3 parametrization with its extremely hard gluon density. Note that, for 
example, WHIT4 and WHIT6 now differ by a factor of 5 for x around 0.3. The gluon 
distribution of WHIT6 is rather similar in shape to the one of LAC1, but significantly 
smaller in magnitude. Indeed, in all three LAC parametrizations, gluons carry significantly 
more momentum than quarks for Q 2 < 20 GeV 2 ; this is counter-intuitive | 43[ 



since m 



known hadrons, and hence presumably in a VMD-like low— Q 2 photon, gluons and quarks 
carry about equal momentum fractions, while at very high Q 2 the inhomogeneous evolution 
equations (H) predict that quarks in the photon carry about three times more momentum 
than gluons. Notice finally that GRV predicts a relatively flat gluon distribution. This 
results partly from the low value of the input scale Ql = 0.25 GeV 2 , compared to 1 GeV 2 for 
DG and 4 GeV 2 for WHIT and LAC1; a larger Q 2 jQ\ allows for more radiation of relatively 
hard gluons off large— x quarks. Recall also that their (pionic) input distribution includes a 
valence-like (hard) gluon density. 

Finally, Fig. 4b shows that some parametrizations also differ substantially in the flavour 
structure. Both DG and WHIT assume d 1 = s 1 at all Q 2 . Moreoever, DG assumes g^o oc e 2 
for the entire input quark distribution, while WHIT assumes this only for the valence input. 
Nevertheless the smaller value of Qq assumed by DG leads to a strangeness content quite 
similar to that predicted by WHIT: At very low x, sea quarks dominate, which have u 1 = 
(P = s 7 , so that s 7 /(w 7 + (P) ~ 1/2; at high x, q? oc e 2 t , so that s 7 /(w 7 + (P) ~ 1/5.0 Since 
sea quarks are produced by gluon splitting, the transition between the sea-dominated and 
valence-dominated regions depends on G 1 . 

In contrast, GRV assumes Uq = and Sq = at the input scale; for given F^, the 
former assumption increases u 1 + d 1 and the latter reduces s 7 , compared to the ansatz 
qj oc e 2 t . This explains the smallness of the strangeness content of the photon predicted 
by GRV, which persists to surprisingly large values of x. It is worth mentioning that AGF, 
which is otherwise quite similar to GRV [at least for the standard normalization of the 
nonperturbative contribution, K — 1 in eqs.flTS])], also assumes the input valence quark 

"The ratio of 4:1 between large— a: (valence) u and d quark densities is implicit to the structure of the 
WHIT fits, while in DG, it 7 and d 1 are parametrized independently; this explains the ~ 2% deviation 
between the two parametrizations shown in Fig. 4b at large x. 
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distributions to be proportional to the squared quark charges; its predictions for s 7 / '{u 1 + (P) 
are therefore quite similar to those of the WHIT group. GS falls between WHIT and GRV, 
since here a small but nonzero Sq is assumed, from the pionic sea quarks as well as the 
QPM part in eq. flU]) . On the other hand, LAC treats all four quark input distributions as 
completely independent quantities, without imposing any constraints between them. This 
results in the erratic behaviour of the strangeness content depicted in Fig. 4b; the ratio 
even exceeds unity for x ~ 0.1. Clearly LAC should therefore not be used when the flavour 
structure of the photon is important, e.g. in W and Z boson production at HERA. 



3) Resolved Photon Processes in jp Scattering 

In this section we discuss jp scattering reactions that are sensitive to the hadronic structure 
of the photon. Most of our numerical results will be for the ep collider HERA, with appro- 
priate (anti-)tagging conditions for the outgoing electron, in order to make sure that the 
virtuality of the exchanged photon is small; we will also make a few comments on fixed tar- 
get 7p scattering. We follow the terminology of ref. ^j7§ in distinguishing between direct and 
resolved photon contributions to a given process; the existence of these physically distinct 
contributions had already been emphasized in ref. f42|. The former are defined as reactions 



where the incident photon participates directly in the relevant hard scattering process; some 
examples are shown in Fig. 5a. In contrast, in resolved photon contributions the incoming 
photon takes part in the hard scattering via one of its constituents, a quark or gluon, as 
shown in Fig. 5b. Notice that in direct processes the entire photon energy goes into the 
hard partonic final state, while in resolved photon processes only a fraction x 7 of the photon 
energy is available for the hard scattering. This also leads to different topologies for the two 
classes of contributions. In LO, the direct contribution to jet production will give rise to 
two high— Pt jets and the remnant jet from the proton, see Fig. 6a; there will be little or 
no hadronic activity in the direction of the incident photon. In contrast, resolved photon 
contributions are characterized by a second spectator jet, which is formed when a coloured 
parton is "taken out" of the photon; this photonic remnant or spectator jet will usually 
go approximately in the direction of the incident photon, which coincides with the electron 
direction at HERA. This additional jet in principle allows distinction between direct and 
resolved photon contributions on an event-by-event basis. 

Occasionally one sees the erroneous statement in the literature that resolved photon 
contributions are NLO corrections to the corresponding direct process. However, by counting 
powers of coupling constants in Figs. 5a and 5b one can easily convince oneself that this is 
not the case. The direct contribution to di-jet production is obviously of order a em a s . 
The resolved photon contribution is of order q 7 ■ a 2 , where q 1 stands for any photonic 
parton density. We have seen early in Sec. 2 that these densities are of first order in a em . 
Moreover, we also saw that they grow logarithmically with the relevant scale Q 2 of the given 
process; in leading-log summed perturbation theory this log Q 2 has to be counted as a factor 
l/a s , see eq.(|^).[] Altogether the photonic parton densities therefore have to be counted as 
O (a em /a s ), so that q< ■ a 2 is again of order a em a s , just like the direct contribution. It 

* Recall that this expression remains approximately correct even beyond the "asymptotic" approximation, 
if we allow a weak Q 2 dependence of the functions a and b. 



15 



should be emphasized that LO QCD is always leading log summed, i.e. a s (Q 2 ) ■ logQ 2 is 
always counted as 0(1), not 0{a s ). To give one example, the resummation of such leading 
logarithms leads to the scaling violations of nucleonic structure functions; it should be clear 
that one is not performing an NLO analysis by simply using Q 2 — dependent parton densities 
when estimating pp cross-sections. 

Part of the confusion is caused by the fact that in NLO, direct and resolved photon 
contributions mix. This can, e.g., be seen from the Feynman diagram of Fig. 7, which 
shows an NLO contribution to direct jet production. If no restrictions are imposed on the 
transverse momentum of the outgoing antiquark, it could go in the photon direction and 
thus form a "remnant jet". In fact, the total contribution from the diagram of Fig. 7 to the 
inclusive jet (pair) cross-section will be dominated by configurations where the transverse 
momentum of the antiquark is small, due to the 1/t pole associated with the exchanged 
quark. The crucial point is that the contribution from this pole has already been included 
in the resolved LO qq' scattering process, where the exchanged quark is treated as being on- 
shell. In order to avoid double-counting, the contribution of the 1/t pole therefore has to be 
subtracted from the NLO contribution of Fig. 7; in other words, the corresponding collinear 
divergence is absorbed in the quark distribution function in the photon. We emphasize 
that this treatment is completely analogous to the calculation of NLO corrections to jet 
production in pp scattering. For example, the incident photon in Fig. 7 could be replaced by 
a gluon coming from the p. The above argument then tells us that qq' and gq scattering mix 
in NLO, i.e. a part of the NLO contribution from gq scattering has to be absorbed in qq' 
scattering; nevertheless nobody would consider one to be an NLO correction to the other.[] 

We do not attempt to split resolved photon contributions into those coming from the 
"anomalous", "pointlike" or "perturbative" part of q 1 and those due to the "hadronic" or 
"nonperturbative" part. We have seen in the previous section that it is not easy to separate 
these parts consistently; indeed, it should be clear that in reality there is a smooth transition 
from the one to the other. Nevertheless we will see later (in Sec. 5) that the existence of the 
"pointlike" contribution may have some impact on overall event characteristics. 

After these preliminaries, we are ready to discuss various hard jp reactions. We start 
with jet production in Sec. 3a, where both NLO calculations and high-energy data from the 
ep collider HERA are available. Open heavy quark production (Sec. 3b) has also been treated 
in NLO, and first HERA data have started to appear. Direct photon production (Sec. 3c) 
also also been treated in NLO, but no HERA data have yet been published. We then discuss 
the production of J/ij) mesons in Sec. 3d, and lepton-pair (Drell-Yan) production in Sec. 3e. 



3a) Jet Production in Collisions 

The production of high— px jets offers the largest cross-section of all hard jp scattering 
reactions. It was therefore the first such process for which resolved photon contributions 



were calculated [58], and also among the first reactions to be studied experimentally at 



tThe dia gram of Fig. 7 also has a divergence when the exchanged gluon goes on-shell; this is absorbed 
in the gluon density in the proton, i.e. in direct jg scattering. Finally, in principle both the exchanged 
quark and the gluon can be (nearly) on-shell. However, in this case none of the final state partons has 
large transverse momentum; this configuration therefore only contributes to soft processes, which cannot be 
treated perturbatively. 



16 



HERA |5m . Moreover, many properties of direct and resolved photon contributions 
to jet production carry over to the production of heavy quarks and direct photons, to be 
discussed in subsequent subsections. We therefore treat jet production in somewhat more 
detail than other photoproduction processes. 

Feynman diagrams contributing to jet production in 'jp scattering are sketched in Fig. 5. 
In leading order (LO), direct contributions (Fig. 5a) come from either the "QCD Compton" 
process (left diagram), or from photon-gluon fusion (right diagram); notice that, unlike 
DIS, this latter process is sensitive to the gluon content of the proton already in LO. The 
resolved photon contributions (Fig. 5b) involve the same matrix elements for 2 — > 2 QCD 
scattering processes that appear in calculations of jet production at purely hadronic (pp or 
pp) scattering. Note that, in contrast to the direct processes, many of these QCD scattering 
processes can proceed via the exchange of a gluon in the t— or u— channel; processes like 
99 ~ * QQ that only proceed via the exchange of a quark in the t— or u— channel, and/or the 
exchange of a gluon in the s— channel, contribute only little to the inclusive jet cross-section. 

In LO, the cross-section for the electroproduction of two (partonic) jets with transverse 
momentum p T and (pseudo)rapidities^] r)%, rj 2 can be written as: 

d3g( 77!T* ) = e fMUx /^t i&) ■ (is) 

dp T dr]idr]2 dt 

Here, i stands for a photon, quark or gluon, and j, k, I stand for a quark or gluon. If 
i — 7 (direct contribution), the function /j| e is just the Weizsacker- Williams photon flux 



TT|) ; otherwise it is given by 



fi\e(Xe) = I —f~f\e(z)fi\~f (— ) . (17) 



Z V z 



The resolved photon contribution to the cross-section ( |T5D therefore only depends on the 
product x e of the fraction z of the electron energy carried by the incident photon, and the 
fraction jz) of the photon energy carried by the parton in the photon. Of course, 

/j| 7 and fj\ p are nothing but the parton densities in the photon and proton, respectively. 
The (pseudo)rapidities of the jets are related to these Bjorken variables by: 

x P = \x T M{e^+e^)- (18a) 

xe= l xT ^¥s e ~ m+e ~ m ^ (i8b) 

with xt = 2pT/\fs, where is the ep centre-of-mass energy; note that in our convention 
positive rapidities correspond to the direction of the proton. Finally, the subprocess cross- 
sections cr for the direct |6l], and resolved photon |>2| contributions depend on the 
Mandelstam variables describing the hard partonic scattering, with s = x e x p s and 

(19) 




* There is no difference between the partonic rapidity and pseudorapidity in LO. 
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both solutions in eq.(JKS|) have to be included when evaluating eq. (|T^|) , which can be accom- 
plished by simply symmetrizing the subprocess cross-sections under t «-> u. 

Direct and resolved photon contributions to the photoproduction of jets were first com- 
pared by Owens [58] in 1979 for fixed target energies, where the resolved photon contributions 
were found to be subdominant, except at very small pt- Although a few theoretical analyses 
[[d3|, Q of the photoproduction of jets appeared in the first half of the 1980's, the importance 
of resolved photon contributions was fully appreciated only in 1987, when it was realized 
[ |5"7| , |5T)| ] that at HERA energies, they could exceed the direct contribution by as much as 



a factor of ten at px — 5 GeV, and remained dominant out to pr — 35 to 40 GeV. 

An update [67] of this result is shown in Fig. 8, where the ratio of resolved photon and 
direct contributions to the single jet inclusive cross-sectionjj] is plotted for the nominal HERA 
energy ^Js = 314 GeV. Unlike in ref. |66|, we have imposed some acceptance cuts, taken from 
a recent ZEUS analysis |68j: The jet has to fall in the pseudorapidity range —1 < 77j Ct < 2, 
and the "antitag" requirement that the outgoing electron is not seen in the main detector 
implies that the photon virtuality Q 2 < 4 GeV 2 . Most results in Fig. 8 have been obtained 
using the MRSD-' parametrization |)9| for the parton densities in the proton; comparison 



between the solid and the dotted curve, which is for the MRSDO' parametrization, shows 
that even the pre-HERA uncertainty of nucleon densities only leads to an uncertainty of a 
few percent in the ratio of Fig. 8, except at very small pt- HERA data on deep-inelastic 
scattering have since then improved our knowledge of the nucleon structure considerably; the 
impact of these data on parametrizations of G p (x) is still under investigation |70fl , but it is 
already clear that very soon the uncertainty from the parton densities in the nucleon will be 
negligible compared to the differences between predictions based on various parametrization 
of (p. In particular, HERA data clearly favour MRSD-' over MRSDO'; we therefore take the 
former as our standard choice. 

We see that implementation of the acceptance cuts reduces the region where resolved 
photon contributions are dominant to pr < 25 to 30 GeV; this is mostly due to the upper 
limit on 77j e t, which reduces the resolved photon contribution much more than the direct one 
(see below). Nevertheless, the former still exceeds the latter by a factor between 5 and 11 
at pt = 5 GeV; at present our lack of knowledge of q 1 does not allow us to predict this ratio 
more precisely. This dominance of resolved photon contributions at small pt can partly be 
explained by the fact that they get contributions from gluon exchange in the t— channel, see 
Fig. 5b; this enhances the squared matrix elements for resolved photon processes by a factor 
s/\t\ > 2, compared to those for direct processes. Colour factors generally also favour the 
former over the latter. Finally, we saw in Figs. 4 that for small x 7 , the parton densities in 
the photon can actually exceed l/a em substantially; for small x e , which contribute only at 
small pt, the integral in eq. (|i~7]) can therefore enhance resolved photon contributions even 
further .0 On the other hand, this convolution integral decreases more rapidly with increasing 
x e than the photon flux factor / 7 i e does; this explains the more rapid decrease of resolved 
photon contributions with increasing pt, and thus the shape of the curves in Fig. 8. 

Obviously parametrizations with sizable G 1 (LAC1, WHIT4) predict a considerably 
larger resolved photon contribution at small pt than those with smaller G 7 [WHIT1, WHIT3; 



^Recall that by definition, events with two accepted jets count twice here. 

^Recall that the factor of a em contained in q 1 is cancelled in the ratio by the explicit factor of a e 
appearing in the a for direct processes. 
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DG, GRV, AGF and GS2 (not shown) also belong to this class]. Note that the prediction 
from the WHIT1 parametrization exceeds that using WHIT3 even at the smallest px shown; 
this indicates that the pr distribution integrated over jet rapidities is not sensitive to very 
small Xj, where the WHIT3 gluon density exceeds that of WHIT1, see Fig. 4b. Finally, we 
should warn the reader that additional experimental cuts can change the ratio of direct and 
resolved photon contributions substantially. This is demonstrated by the dot-dashed curve, 
which has been obtained with the same parametrizations of parton densities in the proton 
and photon as the solid curve, but where we have demanded that the outgoing electron be 
detectable in the ZEUS luminosity monitor; this implies Q 2 < 0.01 GeV 2 and, more im- 
portantly, 0.2 < z < 0.75. The lower cut on the scaled photon energy z greatly reduces 
direct contributions at low pr, while the upper cut reduces resolved photon contributions 
at high p T more than direct ones; as a result, the pt~ dependence of the ratio of the two 
contributions becomes considerably steeper. 

In Fig. 9 we further split the resolved photon (9a) and direct (9b) contributions depending 
on whether the two high— p T partons in the final states are of two quarks, two gluons, or 
a quark and a gluon. Since the contributions from qq —* gg and gg — *> qq are very small, 
the curves in Fig. 9a can also be read as coming from qq, gg and qg initial states, while in 
Fig. 9b the qg and qq final states come from jq and ^yg initial states, respectively. We have 
used the WHIT1 parametrization for q 1 , and applied the same acceptance cuts as in Fig. 8. 

WHIT1 assumes a relatively small G 7 ; as a result, the gg final state is dominant only for 
Pt < 5 GeV. However, even the WHIT4 parametrization, which assumes a two times larger 
input distribution Gq, predicts the cross-section for the gg final state to be well below that 
for the qg final state for p? > 10 GeV. Note that in the majority of qg events the quark 
comes from the photon and the gluon from the proton; this is partly because the photon is 
assumed to be relatively poor in gluons, and partly because the cut r]^ et < 2 removes many 
events with a gluon in the photon in the initial state, as discussed below. We saw in sec. 2 
that the quark densities in the photon are much better known than G 7 , except at small 
Xj. Together with the result of Fig. 9a, this explains the rapid convergence of the curves in 
Fig. 8, although some difference between the LAC1 and WHIT predictions persists even at 
large pt (see also Fig. 3|) . Finally, the comparison of Figs. 9a and 8 shows that two quark 
final states dominate resolved photon contributions only for values of pt where the total jet 
cross section is already dominated by direct contributions. 

As expected from Fig. 8, the direct contributions shown in Fig. 9b have a considerably 
flatter transverse momentum spectrum than the resolved photon contributions. Unfortu- 
nately, photon-gluon fusion (the two quark final state) dominates the direct contribution 
only at relatively small pt', this is because the gluon density in the proton is softer (i.e., 
decreases more rapidly with increasing x p ) than the valence quark distributions. Since at 
small pr the inclusive jet cross section is dominated by resolved photon contributions, a 
direct study of photon-gluon fusion, which might allow to further constrain the gluon den- 



sity in the proton even at rather small x p []7I|, will be difficult unless the resolved photon 



contribution can be suppressed by additional cuts. 

As shown in Figs. 10, the (pseudo)rapidity distribution of the jet can be used to help 



^Notice, however, that all quark flavours contribute equally to the jet cross-section, while contributions 
to Fr[ are weighted by the squared charge. Two parametrizations can therefore have very similar i*! 7 and 
yet lead to different predictions for quark-initiated jet production at HERA. 
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disentangle direct and resolved photon contributions. In Fig. 10a we have used the present 
HERA energy y/s = 296 GeV, and applied the same cuts on the parton level that the ZEUS 
collaboration applied |j8| on their reconstructed jets. Unfortunately the px cut chosen is 
still too low to allow a direct comparison between partonic and jet cross sections, since 
a substantial part of the jets still comes from the "underlying event" (beam fragments, 
initial state radiation, and possibly multiple interactions producing minijets ]7^, |73|; see 
sec. 5); these effects are expected to become more important as one approaches the proton 
beam direction, which corresponds to i] = +oo in our convention, and can thus distort the 
jet rapidity distribution compared to the parton-level distribution shown in Figs. 10. We 
nevertheless expect that the differences between the curves shown in these figures will not 
be washed out by the underlying event. 

We see that in this single-differential cross section the direct contribution is always 
subdominant; however, one can enhance its relative importance by requiring ^j et < 0. Of 
more interest for us is the opposite region of positive, sizable i]^ et , which is sensitive to small 
x e , see eq.( |18b| ), and thus allows to probe the parton densities in the photon at small x 7 , 
see eq.([H]). Indeed, the figure shows that in this region one is most sensitive to differences 
between the various parametrizations of q 7 . We also observe again that these differences are 
much larger than even the pre-HERA uncertainty from nucleonic structure functions. 

We have already noted above that the jet rapidities only depend on the product x e = 
z ■ x 7 ; events with large i]j Ct can come from soft partons in hard photons (small x 1 , large 
z), but also from hard partons in soft photons (large x 1 , small z). Fortunately, ZEUS has 
demonstrated the ability to (approximately) determine z purely from the longitudinal 
momenta of particles in the main detector, without having to detect the outgoing electron 
in the luminosity monitor, which would reduce the accepted event rate by about a factor of 
four. This should allow one to increase the lower cut on z from 0.2 to 0.5, which enhances the 
relative importance of contributions with small x 7 , so that the differences between predictions 
based on different ansatze for q 1 become larger, as illustrated in Fig. 10b. (HI prefers to 
only use photoproduction events where the outgoing electron is tagged in the small-angle 
detector; z can then be determined from its energy.) Although this stronger cut on z reduces 
the resolved photon contribution by slightly more than a factor of two for 77j et > 1, it should 
still enhance the discriminative power of this measurement, given that the ZEUS data sample 
with the looser cuts contains almost 20,000 events with reconstructed jets. 

The integration over the rapidity of the second jet in single-jet inclusive cross sections 
leads to a substantial spread in x e and x p , see eqs.(f[8|)H; in order to further increase the 
sensitivity of the jet rate to the region of small x 1 in general and the gluon density in 
the photon in particular, one therefore has to study more differential cross-sections. In 
Figs. lla,b we show predictions for the triple-differential di-jet cross section, as given by 
eq.([T6|), for yfs = 314 GeV, px = 10 GeV and rji = r] 2 = r\. As in Fig. 9a we display resolved 
photon contributions with different final states separately, as predicted from the WHIT1 
parametrization. In Fig. 11a we have applied the antitag cut Q 2 < 4 GeV 2 on the virtuality 
of the photon, but we have not restricted the allowed range of the scaled photon energy z. 
We see that the direct contribution now dominates at 77 < 0, and remains sizable even for 
i] = 2; in the region of negative rj it is mostly due to photon-gluon fusion. Note that our 

^The integration over pr is less important here, since most events will have pr — pr.cut anyway. 
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choice 771 = r/2 implies t = u = — s/2, which minimizes the dynamical enhancement factor 
s/\t\ of most resolved photon contributions. Together with the slightly larger value of px 
this explains why direct contributions appear more prominent in Fig. 10a than in Fig. 9a. 

Contributions from the gluon in the photon peak at rj ~ +2. The reason is that a 
gluon usually only carries a rather small fraction of the photon energy, see Fig. 4a; for given 
s > 4p^, a large contribution to the total energy in the hard sub-process then has to come 
from the proton, leading to a rather strong boost of the high— partons in the proton 
direction. The quarks in the photon can be substantially more energetic, and therefore 
already contribute at rj ~ 0. 

Although according to the WHIT1 parametrization, the contribution from gluons in the 
photon is clearly enhanced at large 77, at rj = 2 slightly more than half of the total cross- 
section still comes from quarks in the photon or direct processes. As before, the sensitivity to 
the region of small x 7 can be further enhanced by imposing cuts on the scaled photon energy 
z. This is illustrated in Fig. lib, where we have required 0.3 < z < 0.8. The lower limit has 
been chosen such that the direct contribution vanishes identically for 77 > 0, see eq. ( |18b| ). 
(Recall that x e = z for direct contributions, while x e = zx 1 otherwise.) Moreover, the 
relative importance of contributions from gluons in the photon is clearly enhanced compared 
to Fig. 11a; in particular, the gg final state is now dominant for 77 — 2, where only ~ 20% 
of the total cross section is predicted to come from quarks in the photon. Recall that G 1 
might well be larger than assumed in the WHIT1 parametrization, so that the contributions 
from gluons in the photon might be even larger. It should therefore be possible to derive 
stringent constraints on G 1 even using jets with relatively large pt, by focussing on events 
with large rapidity and large z. Of course, the sensitivity to resolved photon contributions 
involving soft partons in the photon is even larger J7J] if one can measure a cross-section for 
fixed z; however, this will need a rather large event sample. 

If z and both jet rapidities have been measured, one can in principle reconstruct both x p 



and x~, using eqs.(18). In practice, HERA experiments use an estimator for x 7 to separate 



the events into "direct" and "resolved" samples; this estimator reduces to x 7 in LO QCD, 
but only uses the measured rapidities and transverse energies of the jets, and can thus also 
be computed in NLO (where it will differ from the partonic x 7 ). The ZEUS collaboration 
|]75| has used 284 reconstructed di-jet events from the first run of HERA to show that the 
x 7 distribution has a peak expected from direct contributions. The HI 

collaboration [76 1 has gone one step further and subtracted the contribution from quarks 
in the photon, estimated using the GRV parametrization. They find evidence for a non- 
vanishing contribution from gluons in the photon only for x 7 < 0.2, thereby ruling out the 
LAC3 parametrization (see Fig. 4a). HI reconstructs G 1 in the range 0.02 < x 7 < 0.2 at 
an effective scale /x 2 ~ 60 GeV 2 ; the result is in good agreement with GRV, but disfavours 
LAC1. However, a fairly sophisticated Monte Carlo analysis was necessary to extract G 1 
even using the LO formalism. For example, we saw in the discussion of Fig. 10a that parts 
of the "underlying event" contribute to the reconstructed jets, which obscures the relation 
between hard partons and jets. The HI analysis |76j depends on such details quite sensitively. 
It might therefore be somewhat premature to exclude the LAC1 parametrization on the basis 
of this evidence alone. 

The ZEUS collaboration has also published |77J an updated jet analysis, based on 12,000 
reconstructed di-jet events with Er(jet) > 6 GeV. They study da/dr] for 4,000 events with 
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— f]2 1 < 0.5, where rj is the average pseudorapidity of the two jets; this quantity is 
closely related to the double differential cross section d 2 a/ (drjid^) at rjx = r] 2 . They find 
reasonable agreement between their data and LO QCD predictions in the "direct" sample 
(events with large measured x 7 ), but most parametrizations of cp give too low predictions 
for the "resolved" sample. This again indicates that at this rather low value of Ex the 
"underlying event" plays quite an important role; as we will discuss in more detail in Sec. 5, 
this aspect is not well described by present standard QCD Monte Carlo event generators. 

Direct and resolved photon contributions are also expected to have different distributions 
in cos#*, where 8* is the cms scattering angle [[78]. Due to diagrams where a gluon is 
exchanged in the t— or w— channel, resolved photon contributions are more strongly peaked 
at small 8* than direct contributions. The increasing importance of the latter over the former 
at higher p T means that the cos#* distribution of di-jet events with large px, will be flatter 
than at low px- 

As remarked earlier, the most obvious distinction between direct and resolved photon 
events is that only the latter contain a photonic spectator (or remnant) jet, see Fig. 6. In 
Fig. 12 we show the average energy of this jet in the lab frame, as predicted from the WHIT1 
parametrization; on the parton level, this energy is simply given by E e ■ z ■ (1 — x 7 ). For 
negative rj (of the high— px jets), the spectator jet is rather soft, since both z and x 7 have to 
be large, so that 1 — x 7 is quite small. For rj ~ 0, events with a gluon from the photon usually 
have quite large z, but moderate x 7 , yielding a high spectator jet energy; it declines at large 
rj since the average z becomes smaller. Quark-initiated events typically have considerably 
larger x 7 and hence a softer photonic spectator jet. Obviously the average spectator jet 
energy will increase (decrease) if a lower (upper) cut on z is applied. 

In principle the results of Fig. 12 offer another possibility to enhance the contributions 
from gluons in the photon, by requiring the presence of an energetic remnant jet in the 
electron direction. However, in practice the energy of this jet cannot be measured very 
accurately, since some part of it will usually be lost in the beam pipe. It should nevertheless 
be emphasized that the first analysis of jet data [59| taken during the HERA pilot run found 
substantial energy deposition in the backward calorimeter even if all high— px jets have 
positive rapidities; this can be understood only if the photonic remnant jet is included in the 
MC simulation. This jet has recently been studied in more detail by the ZEUS collaboration 
f{J; we will discuss their results in Sec. 5. 

Recently new data from the Fermilab fixed target photoproduction experiment E683 have 



been published p0| . It uses a tagged photon beam with mean lab energy of 260 GeV, giving 
a mean -^/s of slightly over 20 GeV. A Monte Carlo analysis suggest that their di-jet sample, 
required to have two reconstructed jets with average px > 4 GeV, gets approximately equal 
contributions from direct and resolved photon contributions, in agreement with theoretical 
expectations. However, unlike at HERA, no direct experimental evidence for the existence 
of resolved photon contributions could be established (other than the overall event rate). 
In particular, the hadronic energy flow in the very forward direction was found to be quite 
similar for direct and resolved photon events, and even for di-jet events produced from a pion 
beam; this somewhat counter-intuitive result can be explained in terms of jet fluctuations. 

So far all our predictions have been computed in LO in QCD. As well known, the overall 
normalization of such predictions is uncertain, since in the leading log approximation used 
here, one cannot with certainty determine the values of the factorization scales in the parton 
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distribution functions and the renormalization scale appearing in the running QCD coupling 
constant. These scales have no physical significance; predictions would be independent of 
them if all orders of perturbation theory could be summed. One therefore expects reduced 
scale dependence already in next-to-leading order (NLO). Moreover, many quantities can 
be predicted meaningfully only if one allows at least three high— partons in the final state; 
these include the dependence of jet cross-sections on the jet definition, and the distribution 
in the transverse opening angle between the two (hardest) jets in events with (at least) two 
jets. 

The first step towards a full NLO calculation of jet photoproduction was taken in 1980 



with a calculation |yj of direct 2-^3 cross sections (e.g. — > qgg and 75 — > qqg). A first 
complete NLO calculation of the direct contribution, including virtual (1-loop) corrections, 
was performed |82| in 1986. These results were applied to jet production at fixed target 



energies in ||83|| , and at HERA energies in [jS3}| . An NLO prediction of resolved photon 
contributions to jet production become possible only after corrections to the hard partonic 
QCD cross-sections had been computed []55| . These results were applied to jet production 



from resolved photons at HERA in ref.|86|]. Finally, in refs. ||87|, pq , p9| |90| , complete NLO 
calculations for single-jet inclusive jet cross-sections were presented, including both direct 
and resolved photon contributions. 



Typical results are presented in Fig. 13, adapted from Bodeker et al. ||S9|| . We show 
the scale dependence of the predicted jet cross-section at HERA for Et = 25 GeV and 
r/jet = 1.5. Note that Et and pt are in general no longer identical in NLO, since now a jet 
might be made up of two partons. In Fig. 13 two partons have been merged into a single jet 
if AR = ^(Ar/) 2 + (A0) 2 < 0.7. 

The solid and short dashed curves have been obtained by setting all three scales (the 
renormalization scale, and the factorization scales for the photonic and nucleonic parton 
densities) equal to each other. As expected, the LO prediction exhibits a much stronger 
scale dependence than the NLO result. It is worth noting, however, that for the "natural" 
choice /i = pt (or Et in NLO), LO and NLO predictions almost coincide. Generally the 
difference between LO and NLO predictions was found to be quite small at HERA energies 
if fi = pt has been chosen and jets are defined with AR between about 0.7 and 1.Q This 
means that the results of Figs. 8-12 should not be affected too much by NLO corrections. 

As already discussed in the beginning of this section, in NLO the distinction between 
direct and NLO contributions is blurred. Diagrams like the one shown in Fig. 7 contribute 
to direct jet production in NLO, but they also contain a logarithmically divergent piece 
which has already been included in the LO resolved photon contribution; this piece there- 
fore has to be subtracted from the NLO direct contribution. This subtraction term grows 
logarithmically with the photonic factorization scale M 7 , which is also the scale appearing in 
the photonic parton distribution functions. In NLO the direct contribution (dotted curve) 
therefore decreases with increasing M 7 , while the resolved photon contribution (dot-dashed 
curve) increases. The sum of the two (long dashed) is nearly independent of M 7 ; the depen- 
dence does not cancel completely since in NLO the subtraction term is exactly proportional 
to logM 7 , while the photonic parton densities only increase approximately like logM 7 , as 

* There appears to be some discrepancy between the results of ref.|8^| and the earlier calculations |8(| p7[ ; 
this is now being sorted out (M. Greco, private communication). 
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discussed in Sec. 2. 

We saw above that di-jet cross sections that are differential in both jet rapidities are more 
powerful discriminators of photonic parton densities than single-jet inclusive cross sections. 
Unfortunately, in conventional NLO calculations going from single-jet to di-jet cross sections 
introduces considerable complications;^ as a result, only the direct production of jet pairs 
has been treated in NLO to date |H] . The urgent need for a full calculation is emphasized 



by the recent ZEUS data fTTH , which have the potential to discriminate between different 
parametrizations of photonic parton densities once the cross section can be predicted with 
some confidence. 

Recently the calculation of jet rates at HERA has been further refined by including two 



additional effects. First, eq. (|18a|) shows that (direct) jet production at negative rapidities 
probes the proton structure at quite low values of x p ; for p T < 10 GeV, partons with x p as low 
as 10~ 3 contribute. Small— x effects may then become important, and one may have to use 
32| the so-called kr factorization |p3R . Secondly, if one defines the photoproduction sample 



with a rather moderate (no-tag) cut on the outgoing electron, one includes contributions 
where the photon virtuality Q 2 may not be entirely negligible; recall, for example, that the 
ZEUS cuts |Q include events with Q 2 up to 4 GeV 2 . On the other hand, when using the 
simple Weizsacker-Williams approximation (|T7|), one assumes Q 2 to be small compared to all 
other scales in the problem. This is still a good approximation for the direct contribution in 
this case, but the parton densities in the photon become suppressed |4j once Q 2 > Aq CD . 
The simple factorization (|17]) then breaks down, but one can still define a "parton density 
in the electron", which will depend on the experimental cut on Q 2 [[JIJ. For the ZEUS cuts, 
this suppression only amounts to a few percent. 

Finally, as mentioned in the discussion of Fig. 10, so far jet production at HERA could 
only be investigated experimentally at rather moderate transverse momenta, where the "un- 
derlying event" can still contribute significantly to the reconstructed jets. The influence of 
the underlying event might be less problematic when one studies the production of high— p? 
particles pB| , |82"|] rather than jets. On the other hand, one now has to specify not only 
parton distribution functions, but also fragmentation functions, before a prediction can be 
made. Also, the cross-section falls off very rapidly with p?, forcing one to work at scales of 
only a few GeV where the convergence of the perturbative expansion might be rather slow. 
Moreover, one either has to experimentally identify different particle species (chiefly pions 
and kaons, if only charged particles are counted), or make assumptions about the relative 
abundances of these species. Recent measurements of the cross-section for the production 
of charged particles with p T > 1.5 GeV, by both the HI and ZEUS collaborations |J7| 



find good agreement with a theoretical NLO prediction |98[ as far as the pt spectrum is 



concerned; the pseudorapidity distribution measured by HI is less well described, but the 
experimental errors do not allow to draw a definite conclusion at this point. 

3b) The Photoproduction of Heavy Quarks 

The production of heavy quarks offers two theoretical advantages over the production of 
light partons (jets) discussed in the previous subsection. First, their large mass trq 3> Aqcd 

^This step should be much easier using the Monte Carlo method of ref.lss). 
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ensures that QCD perturbation theory is applicable in all of phase space, although non- 



perturbative corrections oc (AqcdA^q)™ - might not be negligible for charm quarks. In 
particular, the total cross-section without any cuts could now be predicted with some relia- 
bility if the values of certain parameters (mg, Aqcd and the parton distribution functions) 
were known precisely. Secondly, at least in leading order the number of contributing partonic 
processes is much smaller, making heavy quark production easier to analyze. Specifically, in 
LO only photon-gluon fusion contributes to QQ production from direct photons, while the 
relevant resolved photon processes are gg fusion and light qq annihilation.^ 

In LO, the QQ production cross-section in ep scattering is still given by eqs.(p^)- ([T8|) . 
However, 77^2 now have to be interpreted as true rapidities, which differ from the pseudora- 
pidity for massive particles; moreover, xt in eqs. (|T8|) is now given by l^jpj, + ttlq/ 'y/s , and 
eq.([l9|) has to be replaced by 



s 

t = - 
2 
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4 (m 2 Q + p\ 
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The partonic cross-sections a for 7^ 
PH , [1021 , respectively. 



QQ and gg, qq — > QQ can be found in refs.[|100|| and 



Resolved photon contributions to heavy quark production were first treated in ref.jH]], for 
the case of the top quark, whose mass was then believed to be in the vicinity of 40 GeV. We 
now know ||103| that the top quark is too heavy to be produced at HERA, but cc and bb pairs 
will be produced copiously. It was pointed out in refs. |]57|, 104 | that (for most parametriza- 
tions of q 1 ) resolved photon contributions to the total cross-sections are subdominant, but 
not negligible; e.g., they amount to ~ 20% for the DG parametrization. Resolved photon 
contributions to the production of heavy quark pairs are therefore considerably less impor- 
tant than for high— px jet production. The reason is that now the resolved photon processes 
also only involve gluon exchange in the s— channel or (heavy) quark exchange in the t— or 
u — channel; there is no enhancement factor s/\t\, unlike for jet production. Moreover, the 
two sub-processes involving the parton content of the photon have rather small colour fac- 
tors, again in contrast to the matrix elements appearing in resolved photon contributions to 
jet production. These analyses also showed that gg fusion is predicted to dominate over qq 
annihilation, by a factor of 10 (3.5) in case of the DG parametrization and cc (bb) production. 
The resolved photon contribution therefore offers a good opportunity to constrain G 7 , while 
(in LO) the direct contribution is proportional to G p . 

We saw in the previous subsection that at HERA, contributions from gluons in the photon 
are most important at sizable, positive rapidities, but are suppressed at negative rapidity. 
The rapidity distribution therefore offers a good handle for separating the two contributions 
to QQ production. In Fig. 14 we show the pr spectrum of c and b quarks at central rapidity, 



*It can be argued that at very high transverse momentum, pr ^> m,Q, a s logpr/mQ should be counted 
as O(l), rather than 0(a s ). In this case the "excitation" processes Qg — > Qg and Qq — > Qq also contribute 
in LO, since the Q— quark density in the photon grows logarithmically with the hard scale of the process. 
In these reactions the heavy quark jet is balanced by a light quark or gluon jet, while in QQ creation events 
two heavy quark jets occur with equal and opposite pt- These "flavour excitation" contributions are now 
under study Q. 
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2/1 = 2/2 = 0. We see that the resolved photon contribution (difference between solid and 
dashed curves) only amounts to 10-15% at small pr, and becomes entirely negligible for 
Pt > 5 GeV. In fact, as shown in ref.[|j(J, the resolved photon contribution to heavy quark 
production at HERA can always be suppressed to an insignificant level by requiring (at 
least) one of the two heavy quarks to emerge at y < 0.Q Note that we have chosen the 
WHIT4 parametrization in Fig. 14, which is characterized by a rather large and hard G 7 ; 
it predicts that resolved photon contributions amount to more than 20% (30%) of the total 
cc (bb) cross-section at HERA. This contribution is more important for bb production, since 
the direct contribution is suppressed by the small charge of the b quark; for parametrizations 
with rather hard G 1 (WHIT1,4, DG, GS2, GRV) this suppression is stronger than the relative 
reduction of the resolved photon contributon with increasing itlq, which is caused by the 
additional convolution (|i"7|) with the gluon density in the photon. 

Since resolved photon contributions to QQ production are insignificant at high p T (unless 
their importance is enhanced by specific cuts, as discussed below), the ratio of cc to bb 
cross-sections in the region px > 10 GeV simply reflects the ratio of their squared charges. 
This relative suppression of the bb cross-section means that b— tagging at HERA will be 
significantly more difficult than it is at pp colliders, where bb and cc cross-sections become 
equal at high p?. At pp colliders the harder fragmentation function of b— flavoured hadrons 
[ [105|1 means that inclusive high— p? muon production is dominated by bb events; this will not 
be true at HERA, however, except at very high pt where the cross-section is quite small. 
Charm quarks will also be a serious background to b— tagging by micro-vertex detectors at 
HERA. 

As mentioned earlier, the resolved photon contribution is concentrated at positive rapidi- 
ties. Just as in the case of jet production it can be isolated by a cut on the scaled incident 
photon energy z; e.g., requiring z > 0.3 at px = 10 GeV removes all direct contributions with 
Hi — Vi > 0. Fig. 15 shows that the remaining resolved photon contribution is indeed very 
sensitive to the gluon content of the photon. Even according to the WHIT1 parametriza- 
tion the qq annihilation contribution (short dashed) is considerably below the one from gg 
fusion (long dashed). The WHIT4 parametrization therefore predicts a considerably larger 
cross-section, but the shape of the rapidity distribution is similar to that predicted from 
WHIT1. In contrast, the LAC1 prediction differs in both normalization and shape. How- 
ever, we should warn the reader that without the cut on z, even at yx — yi — 2 the direct 
contribution would be at least ten times larger than the resolved one; the experimental im- 
plementation of this cut therefore has to be very efficient. It might even be necessary to 
require the presence of a photonic remnant jet to extract the resolved photon contribution; 
recall that this jet is expected to be quite energetic in events that originate from the gluon 
content of the photon. Recall also that there will be a large contribution from "charm ex- 
citation" |)9] if only one of the two high— pt jets is tagged as a heavy quark. Finally, the 
cross-section shown in Fig. 15 is quite small, even though we have not yet required any 
specific charm signal (e.g., a hard muon or reconstructed D* meson). Clearly HERA will 
have to accumulate significantly more data than the present 6 pb -1 (as of the end of 1994) 
to measure such triple-differential cross-sections even at lower p?. 

The predictions shown in Figs. 14 and 15 were computed in LO. NLO calculations of 

^This result does not hold for the LAC3 parametrization; fortunately, this parametrization is excluded 
by other data, as discussed in sees. 3a and 4a. 
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the photoproduction of heavy quarks exist |104j, |106|, |107|; unlike for jet production, even 



the fully differential cross-section is available in NLO ||108|| . It has been demonstrated ||109| 



that the direct contribution can be extracted reliably also in NLO, i.e. the sensitivity to 
the gluon content of the proton is not degraded. In view of their smaller size, extraction of 
the resolved photon contribution might prove more difficult. In particular, in direct events 
with an additional hard parton in the final state (jg — > QQg, r yq — > QQq) the heavy quarks 
can occur at large positive rapidity even after a cut on z has been applied; one may have to 
veto the occurence of additional high— p T jets, and/or require the two highest px jets to be 
back-to-back in the transverse plane, in order to efficiently suppress such backgrounds to 
the cross-section shown in Fig. 15. 

The considerable body of data on photoproduction of charm at fixed target energies 
(y/s < 20 GeV) is well described by NLO QCD calculations [ 110|| . However, at these low 



energies the resolved photon contribution is quite small; it is significant only on the backward 
direction (opposite to the incident photon) |pL08|| , where the experimental acceptance is poor. 
Very recently, first data on charm production at HERA have become available. The 



ZEUS collaboration [111] searched for fully reconstructed D mesons. They observe a 



signal of 48 ± 11 events within the acceptance region px(D*) > 1.5 GeV, \rj(D*)\ < 1.5; 
this corresponds to a(ep -> D±*X) = (32 ± 7 ty) nb at y/s = 296.7 GeV with Q 2 < 4 
GeV 2 and 0.15 < z < 0.86. They then attempt to estimate the total cc cross-section 
from this measurement; this, however, sensitively depends on the extrapolation of the cross- 
section into kinematic regions where it has not been measured, which introduces a strong 
dependence of the "measured" cc cross-section on the assumed parton distribution functions 
in the proton and photon as well as on m c . This not only greatly increases the quoted 
(systematic) error; even the central value depends on these assumptions. It does not make 
much sense to compare this "measured" cross-section with different theoretical predictions, 
since "measurement" and "prediction" depend on the same quantities! Notice also that 



the prediction for the total cc cross section suffers from large uncertainties ||112|| . Since 
cts{ m c) is still quite large, the perturbative expansion only converges slowly, which manifests 
itself in a rather strong dependence of even the NLO prediction on the factorization and 
renormalization scales. Moreover, the prediction is very sensitive to m c , decreasing by more 
than a factor of three when m c is increased from 1.2 to 1.8 GeV. Finally, small— a; effects 
might be sizable [ 112 , |113|| . All these sources of theoretical uncertainties are reduced once 
we require the charm quarks, or their fragmentation and decay products, to have significant 
transverse momentum. One should therefore directly compare QCD predictions for the 
cross-section in the experimentally accessible region with the data. 

The same remarks also apply to the as yet preliminary analysis of charm production by 
the HI collaboration [|114|| , which is based on events with a hard muon. They find 484 events 
where at least one muon satisfies pr(/-0 > 1-5 GeV and 30° < 9(fi) < 130°; some 280 of these 
events are expected to contain fake muons, or muons from 7r and K decays. This gives an 
accepted cross-section a(ep — > /i ± V) = (2.03 ± 0.43 ± 0.7) nb; about 95% of this signal is 
expected to come from cc events, the rest coming from bb production. 

We attempted to reproduce the cross-sections measured by the ZEUS and HI collab- 
orations with a parton-level MC generator based on LO QCD expressions. We take the 



renormalization and factorizations scales to be \jvn? c + p\ and m c = 1.6 GeV; as mentioned 
earlier, the px cuts greatly reduce the sensitivity to m c . However, these cuts also introduce 
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an additional difficulty in the theoretical treatment. They are sufficiently high so that frag- 
mentation effects will play a role. On the other hand, pt cannot safely be assumed to be 
much larger than m c here, so factorizing the result into a hard production cross section and 
a fragmentation function may not yet be a good approximation. We therefore ran our MC 
programs with two different options, using the standard Peterson et al. fragmentation func- 



tions [ 105" or no fragmentation at all. In the former case we also include contributions from 
the charm in the photon (qc — > qc and gc — > gc);| after all, the use of both fragmentation 
and structure functions rests on the factorization theorem, so it seems reasonable to treat 
them symmetrically. 

In our comparison with the ZEUS result we always include a factor of 0.26, which is [ 111 



the probability for a charm quark to fragment into a charged D* meson. We find that if 



we leave out both fragmentation and the contribution from the charm in the photon, we can 
reproduce the experimental cross section only if we assume a gluon distribution that increases 
rapidly at small x; MRSD-' works well, while the prediction of MRSDO' is too low by nearly 
a factor of two for all reasonable choices of momentum scale. Due to the rapidity cut, the 
contribution from resolved photon processes (chiefly gg fusion) only amounts to 20% or less. 
On the other hand, if we include fragmentation effects in the standard way, the average px of 
the charm quarks in accepted events increases by nearly a factor of two, thereby reducing the 
sensitivity to very small x; moreover, the total cross section is now actually dominated by 
contributions involving the charm content in the photon, chiefly eg scattering. As a result, 
we can reproduce the experimental cross section using either MRSDO' or MRSD-' partons 
in the proton, and LAC1 or WHIT partons for the photon; the DG parametrization now 
predicts a far too large cross section, since it assumes c 7 = w 7 , which is manifestly a bad 
approximation at these rather low momentum scales. 

The HI sample is less sensitive to the gluon in the proton at small x, due to both the Pri^) 
cut (which leads to a considerably higher mean transverse momentum for accepted charm 
quarks than in the ZEUS data sample), and the requirement that 8(fi) < 130°. The data 
therefore do not even allow to discriminate between MRSDO' and MRSD-' if we ignore both 
fragmentation and the charm content of the photon. In this case the resolved photon contri- 
bution is quite small, so that any (reasonable) combination of photon and proton structure 
functions is in agreement with the data, yielding a LO cross section of about 1.4 to 2.4 nb. 
If we include both charm fragmentation and the charm content of the photon, the sensitivity 
to soft gluons in the proton is reduced even further, and predictions using MRSDO' differ 
from those using MRSD-' by less than 0.1 nb. The LAC1 and WHIT parametrizations yield 
a predicted cross-section of about 1.5 nb, very close to the MRSDO' prediction without frag- 
mentation and without the contribution from charm in the photon. The DG parametrization 
gives a prediction of about 2.9 nb; as already discussed, this parametrization over-estimates 
the charm content of the photon, but even this high number is not inconsistent with the 
experimental result. 

3c) Direct Photon Production in ^yp Scattering 

Interest in the "deep-inelastic Compton process" •yp — > 7X dates back to the early days 
of the quark-parton model [ |115|| ; among other things, the cross-section for the simplest 

■^There is also a very small contribution of this kind from the charm in the proton, which we neglect. 
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contributing sub-process, yq — > 75, depends on the fourth power of the quark charge, and 
thus allows an independent determination of these charges. However, very soon after the 
introduction of the concept of the parton content of the photon it was realized ]|2| that other 
subprocesses contribute to the production of direct photonsf] at high p T : Not only resolved 
photon processes like qq — > #7 and gq — > qy have to be included in a consistent LO QCD 
treatment, but also fragmentation processes, where a high— parton fragments into a hard 
photon and an additional (nearly collinear) jet; note that the q — ► 7 and g — > 7 fragmen- 
tation functions [flip are 0(a em /a s ), just like the parton densities in the photon. A full LO 
treatment therefore has to include direct processes like 79 — > gq — > gqy and resolved photon 
processes like qq' — > qq' — > qq'y, these processes involve the same partonic matrix elements 
that appear in direct and resolved photon contributions to the production of high— px jets. 
Hence direct photon production is in general actually more complicated to analyze theoreti- 
cally than jet production, contrary to statements often found in the literature. On the other 
hand, photons are easier to study experimentally than jets; in particular, their energy can 
be measured considerably more precisely. 

In LO the triple differential cross-section for the production of a direct photon in pro- 
cesses that do not involve parton — > photon fragmentation is still given by eq.(lB); the 
cross-section for fragmentation processes contain a convolution with a fragmentation func- 
tion: 

d 3 a hag (ep -> ejjX) ^ r dz' 



dpTdrj^dr/j 



E J irk 

i,j,k,l 



Z — . 

' dt 



(21) 



Eqs. (P^8|) , (|T9|) also still apply here, but for the fragmentation contribution one has to replace 
Pt by p' T = Pt/z'. Moreover, unlike the case of di-jet production, the final state particles 



are now distinguishable. Note that eqs.fllq) allow a two-fold ambiguity for 771,2 in terms of 
the Bjorken— x variables; taking rji = r/ 7 , one has 



In 




(22) 



In case of jet production one can arbitrarily fix the sign in eq.(22), since this only corresponds 
to the definition of which of the two partons gives "jet 1" . However, in case of direct photon 
production the sign in eq.([22"D is correlated with the choice of sign in t, eq. (|T9|) . In particular, 
if t is defined as the momentum transfer from the incident photon, or parton in the incident 
photon, to the final photon, or parton fragmenting into the final photon, taking a + sign in 
eq. (fl9|) means \t\ > s/2, which implies that one has to take the + sign in eq. (ID; recall that 
we define the proton direction has having positive rapidity. 

The first quantitative estimates of cross-sections for the production of direct photons in 
7P scattering, including all sub-processes listed above, were presented in refs.jfTH], using 



very simple ansatze for parton distribution and fragmentation functions. The first partial 
NLO analysis for fixed-target energies was published in ref. ||41|| , where only NLO corrections 
to ■yq jq, as well as the direct NLO process 75 — > qqy, were included; all other contribu- 
tions were treated in LO. This was justified by the observation that at these (low) energies 
a rather modest pr cut on the outgoing photon suffices to greatly suppress contributions 

*As opposed to photons stemming from the decay of hadrons, mostly 7r° and r\ mesons. 
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involving photonic parton densities and/or fragmentation functions. Notice that only pro- 
cesses involving both the parton content of the photon and parton to photon fragmentation 
can proceed via gluon exchange in the t— or u— channel; one needs to convolute the cross 
section with two additional functions, compared to simple 75 — > 75 scattering, before the 
s/\t\ enhancement characteristic of these gluon exchange processes becomes available. It is 
therefore not surprising that resolved photon contributions are not quite as important in the 
case of high— p T direct photon production as they are for jet production. 

The analysis of ref . |OT was repeated independently, as well as extended to HERA energies, 
in ref. ||1 1 7|| . However, here the resolved photon contributions involving fragmentation (which 
have comparatively soft pt spectra) were omitted, while the contribution from 7$ — > 7$ (via a 
box diagram ||118|| ) was included; this last contribution was shown to be significant in certain 



regions of phase space (see below), even though it is formally of next-to-next-to-leading 
order (NNLO). 

In contrast, the emphasis of refs. ||1 19[| was on processes involving the parton content of 
the photon. In particular, it was shown that a clear signal from the gluon content of the 
photon (largely from gq — > 75) can be extracted from a sample of events with fixed energy of 
the incident photon. Due to the large Lorentz boost from the parton-parton cms to the lab 
frame these photons can be quite energetic if they emerge at large rapidity, even at rather 
small 

A further step towards a full NLO analysis was taken in ref. ||120|| , where corrections to the 
resolved photon contributions were included; however, all fragmentation contributions were 
still treated in LO. Notice that inclusion of NLO corrections to non-fragmentation contri- 
butions are necessary to reduce the dependence on the scale appearing in the fragmentation 
functions. In NLO, fragmentation and non-fragmentation contributions mix, just like direct 
and resolved photon processes do. For example, gq — > gq'j is an NLO correction to gq — > g7, 
but also part of the LO contribution involving q — > 7 fragmentation if the final state q and 
7 are (nearly) collinear. In order to avoid double counting, this collinear (divergent) con- 
tribution therefore has to be subtracted from the NLO correction. This subtraction term is 
proportional to the logarithm of the scale appearing in the fragmentation function. Since 
the subtraction term (obviously) appears with a negative sign in the final result it largely 
cancels the dependence on this "fragmentation scale"; the cancellation is not perfect, since 
the fragmentation function resums all orders of this leading logarithm, while the subtraction 
term does not. On the other hand, since the fragmentation contributions have only been 



treated in LO in ref. ||120|| , a rather strong dependence on the renormalization scale appearing 
in a s remains. 

In Fig. 16 we show some LO predictions for the rapidity dependence of the direct photon 
cross-section, adapted from numerical results of ref. [ 120|| . We see that the "Born" cross- 
section (from 75 — > jq) peaks at negative rapidity (in the direction of the incoming photon); 
even though the 1/u pole of the hard matrix element favours configurations where the final 
photon is emitted in the proton direction, negative rapidities are favoured since they probe 
the proton at small x, where the (sea) quark densities increase quickly. In the resolved 
photon process g 1 q p — > jq the final photon is also preferentially emitted in the proton 
direction; moreover, at large positive rapidity one becomes sensitive to the gluon density in 



''The same is true for very forward jets, of course. However, it should be significantly easier to detect a 
photon at small angle than to reconstruct a jet just a few degrees away from the proton beam. 
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the photon at small x, where it is (presumably) large. As a result, contributions oc G 1 are 



expected to dominate at large, positive r) 7 , as first noted in refs. ||119|| . Notice that resolved 



photon contributions involving parton to photon fragmentation have not been included in 
Fig. 16; however, we already know from the discussion of Sec. 3a that at large positive 
rapidity they are also dominated by contributions oc G 1 . Finally, the box contribution, 
7# — > 7<7, is important only at very negative ?7 7 , where one probes the gluon density in the 
proton at small x. 

The results of Fig. 16 have been computed using LO expressions. While NLO corrections 
are significant (e.g., they change the rapidity distribution of the direct contributons p] , |117|| ), 
they do not change the conclusion that direct photon production offers a good handle for 
constraining G 1 . More recently this conclusion has been challenged on different grounds 



[ [122| , |123| : The integral over the photon spectrum ( |Tl|) tends to smear out the distributions, 
which in Fig. 16 were shown for a fixed energy (10 GeV) of the incident photon. In particular, 
direct contributions from rather soft initial photons also populate the region of positive 
rapidity. These papers also complete the NLO calculations by including corrections to the 
fragmentation processes, but these corrections have little bearing on the question whether 
direct photon production is useful for constraining the gluon content of the photon. 

Ref. [|123j| introduces two additional refinements. It makes use of an NLO parametrization 
124j] of the parton to photon fragmentation functions, which updates the "asymptotic" ex- 



pressions of ref. [41 1 . More importantly, this analysis for the first time introduces an isolation 



requirement. At fixed target energies, backgrounds from ir and rj decays can be suppressed 



quite reliably, since they contain two photons with (usually) substantial opening angle ||125 
At higher energy, and higher pj,, this opening angle becomes much smaller, making it more 
difficult to detect both photons individually. At pp colliders a direct photon signal could 
therefore only be detected ||126|| if the photons were isolated, i.e. after events were discarded 
if more than some maximal (small) amount of energy was found in a cone around the photon. 
Such a cut reduces the fragmentation contribution considerably; in particular, the contribu- 
tion from g — > 7 fragmentation becomes very small, since here most of the energy usually 
goes into the accompanying jet rather than the photon. Nevertheless, at hadron colliders 
contributions from q — > 7 fragmentation can still be significant ||127|| . Strictly speaking, the 
necessity to impose an isolation cut at HERA has not yet been demonstrated. Indeed, the 
only background study we are aware of [12S] reaches the conclusion that even without isola- 
tion cut the background in the most interesting region of positive rapidity can be suppressed 
to the 25% level, which might be tolerable. However, this study ignores resolved photon 
contributions to the background (they are included for the signal); we saw in Sec. 3a that 
these contributions will increase the total jet cross-section (and hence also the cross-section 
for the production of high— px particles) by a large factor at positive rapidity. In the absence 
of a more complete background study we are therefore inclined to believe that an isolation 
cut will indeed be necessary at HERA. 

In Fig. 17, which has been adapted from numerical results of ref. ||123|| , we show the 
rapidity dependence of a(ep — > ejX) after requiring that the hadronic energy in a cone 



5 = y (Arj) 2 + (A0) 2 / coshr^ = 0.4 around the outgoing photon be less than 10% of the 
energy of that photon. No cut on the energy of the incident photon has been imposed. We see 
that, except at very large 7/ 7 , the cross-section is dominated by the direct contributions and 
resolved photon contributions involving the quark content of the photon only. In particular, 
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in sharp contrast to the "Born" cross-section of Fig. 16, after integration over the incident 
photon spectrum the direct contribution has a very flat rapidity distribution. The authors 
concluded that, while one might be able to further constrain the quark densities in the 
photon from this measurement (see the difference between the dotted curve, obtained from 
the GS parametrization, and the solid one, which is for the GRV parametrization), there is 
very little sensitivity to the gluon content of the photon. 

However, we saw in Sec. 3a that even a rather modest lower cut on the energy of the 
incident photon greatly suppresses direct photon contributions at positive rapidity; see e.g. 
Fig. 10. Such a cut would also reduce contributions oc q] substantially, but would have much 
less effect on contributions oc G 1 . The sensitivity to the gluon content of the photon would 
presumably be enhanced even more if the rapidity of the jet balancing the direct photon can 
be measured as well [|122|1 F1 We therefore believe that the conclusions of refs. [|122| , |123|] are 



probably too pessimistic. On the other hand, HERA will have to accumulate a large body of 
data before a substantial number of direct photon events with known (fixed) energy of the 
incident photon will become available; only then can predictions like those shown in Fig. 16 
be compared with experiment. 

So far no HERA data on direct photon production have been published. Data from a 
fixed target photoproduction experiment exist ||125|| ; they agree with theoretical predictions 



29|| , but the energy is too low to extract a signal for the resolved photon contribution. 
Before closing this subsection we briefly mention some work on closely related topics. 



In ref. [|130|1 the cross-section for c + 7 production has been computed (ignoring fragmenta- 
tion contributions, however), the main motivation being that this might allow to constrain 
the charm content of the photon. In ref. [|131|| it has been suggested that one might learn 



something about the parton content of polarized photons from the study of direct photon pro- 
duction. Indeed, longitudinally polarized electron beams should become available at HERA 
in a few years; part of this polarization will be passed on to photons emitted from these 
electrons, if they carry a substantial fraction of the electron's energy. However, one would 
also have to measure the polarization of the outgoing photon, which seems quite difficult. 
Finally, ref. ||132|| is a first study of the production of two direct photons. The cross- 



section is substantially smaller (by a factor oc e q a era /a s ) than the one for single direct photon 
production. Apart from fragmentation contributions, only the resolved photon process qq — > 
77 contributes in LO; processes like 75 — > 577 are formally of NLO, since q] oc a em /a s , 
and are indeed found to be subdominant numerically. Due to the large charge factor, this 
process might be useful for probing the up-quark density in the photon at small x. If pj. 
values as low as 3 GeV are experimentally accessible, one might even be able to extract the 
contribution from the box diagram, gg — > 77; since it is proportional to the gluon density 
in the photon, it should be concentrated at larger rapidities than the contribution from qq 



annihilation. This close analogue of the famous light-by-light scattering process ||133|| has 
yet to be studied experimentally. 

3d) J /if; Production in ^yp Scattering 

The inelastic production of J/ip mesons in (virtual or real) 7p collisions has long been 
regarded as one of the cleanest methods to constrain the shape of the gluon distribution in 

+ Unfortunately, no NLO calculation for 7+jet production exists as yet. 
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the proton 



1 134]. "Inelastic" here means that the quantity 



Z = 



(23) 



is significantly below unity In the proton rest frame, eq. (|23|) reduces to Z = Ej/^/E^; Z = 1 
therefore means that the incident photon transmits its entire energy to the J / if) meson. The 
cross-section for inelastic J /if) production is nowadays usually computed using the "colour 
singlet" model of ref.| 135|| , see Fig. 18. In this model one computes the cross-section for 
jg — > ccg, and projects out the contribution where the cc system is in a colour-singlet, 
s— wave, J = 1 state. One further assumes that the relative momentum of the c and c in the 
J /if> are negligible; the matrix element for J /if) production is then proportional to the wave 
function at the origin \I/(0), which can be determined from the leptonic decay width of J /if). 

This model seems to describe the Z distribution of fixed target data better | 136|| than 
the alternative "dual model" 



TUT 



The agreement in the high— Z region becomes even 
better if one introduces a small but nonvanishing relative momentum between the c and 
c quarks |138|| . Until recently the overall normalization of the cross-section was not well 
reproduced by theoretical calculations, i.e. a sizable "K-factor" had to be introduced. For- 
tunately NLO corrections to the direct diagram of Fig. 18 have recently been calculated 
[ |139|| . Together with the NLO corrections to the leptonic decay width 



p.40|| they give an 



overall inelastic cross-section in agreement with fixed target photoproduction experiments; 
however, the photoproduction cross-section extracted from fixed target leptoproduction ex- 
periments seems to be somewhat higher ||136|| . The calculation of ref. ||139|| predicts that the 
K-factor should be smaller at high (HERA) energies. 

As usual, resolved photon contributions are small at the energies of fixed target exper- 
iments, but they could be quite important at HERA [141, 142]. In addition to the process 
gg — ► J/tf>g, one also has to consider gg — > Xc and gg — > Xcg, with subsequent decay of the 



heavier \c states into a J /if) and a (soft) photon ||143[| ; the latter process has to be considered 
if one imposes a cut on the transverse momentum of the J /if). For pT^J/if)) > 3 to 5 GeV, 
the contribution from b — > J /if) decays becomes important, and eventually even dominates 
the direct contribution ||144|| . There might also be sizable contributions from eg — > J/ipc 
145|| . However, the typical momentum scale in most J /if) events is too small to reliably use 



charm distribution functions in the photon or proton; we therefore neglect this contribution 
here. 

As usual, resolved photon contributions are characterized by the presence of the photonic 
remnant jet, and by a rapidity distribution that peaks at large, positive values. In addition, 
direct and resolved photon contributions have very different Z distributions. The direct 
contribution is peaked at Z ~ 1, which corresponds to a small energy of the outgoing 
gluon. In contrast, for fixed x 7 the resolved photon contribution peaks at Z ~ a; 7 ; after 
convolution with the photon spectrum (|TTJ) this leads to a Z distribution that quickly rises 
with decreasing Z. 

Since only 



This is shown in Fig. 19, which we adapted from numerical results of ref. 1 146 



a mild cut on the px of the J / if) meson has been applied, the contribution from b decays is 
relatively small and has been neglected. LO expressions have been used everywhere; neither 
NLO nor non-relativistic corrections to the resolved photon contribution are as yet known. 



The Bl parametrization of ref. |[L47 ] has been used for the gluon density in the proton, and 
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DG for the photon. Finally, although the cross-section has not been multiplied with any 
branching ratio, it is clear that only the leptonic decays J/i/j —> e + e~, are detectable at 

HERA; the combined branching ratio for these modes is about 12%. The results of Fig. 19 
have been obtained under the condition that both leptons can be reconstructed by the ZEUS 
tracking system. Unfortunately this removes events at large rapidity, which greatly reduces 
the resolved photon contribution. Including the leptonic branching ratio, it only amounts 
to about 7 pb after integration over Z, compared to a direct contribution of about 110 pb.Q 
Requiring Z > 0.2 leaves a very pure direct sample, which should allow to determine the 
shape of the gluon density in the proton for 2-10 ~ 4 < x < 0.1 [[L46|| . 

Unfortunately the extraction of G 1 seems to be much less straightforward. The simulation 
of ref . ||146|| indicated that measurement errors would smear out the direct contribution into 
the region of small Z, totally obscuring the resolved photon contribution. It remains to be 
seen whether an unfolding procedure, and/or the application of additional cuts (e.g., tagging 
of the photon remnant jet) will improve the situation sufficiently to allow to extract new 
information on G 1 from J/ip production at HERA.[] 

Both HI and ZEUS have published first results ||148fl on J/ip production. However, in 
both analyses events were rejected if any particle in addition to the two leptons resulting 
from J / ij) decay was observed. This cut removes most if not all inelastic contributions, but is 
sensitive to elastic or diffractive J / ip production. This allows to test Pomeron-based models 
146| , |138| , |149|| , but teaches us nothing about the parton content of the photon. Very recently, 



ZEUS has announced [|150|1 preliminary results on inelastic J/ip production. However, they 
required Z > 0.2, which again removes most of the resolved photon contribution. 

Finally, we briefly mention associate J /ip + r y production. As first pointed out by Fletcher 
et al. ||142|| , in LO only the resolved process gg — > contributes, the direct contribution 

being forbidden by colour conservation. In principle this process therefore allows a clean 
determination of (the shape of) G 1 | |151| |. Unfortunately the cross-section at HERA is quite 
small, very roughly of order 0.1 to 1 pb after acceptance cuts and multiplication with the 
leptonic branching ratio. 



3e) Production of Lepton Pairs in Collisions 

In this subsection we discuss the production of lepton pairs, either due to the exchange 
of a virtual photon or from the decay of an on-shell W or Z boson. The corresponding 
cross-sections are quite small even at HERA energies, so we will be brief here. 

In the theoretical treatment of these reactions one has to distinguish the cross-section 
integrated over the transverse momentum of the lepton pair (not to be confused with the px of 
the individual leptons) from the cross-section for the production of a high— px lepton pair. In 
the former case, the only LO contribution comes from the resolved photon process qq — > / + /~, 
which produces a lepton pair with vanishing transverse momentum. The corresponding ep 
cross-section is 0(al m /a s ) since, as emphasized repeatedly, the parton distribution functions 
in the photon are (D(a em /a s ). In contrast, both the direct process jq — > l + l~q and the 

* Recall that these are LO predictions, and hence somewhat uncertain; the cross-sections also depend on 
the gluon densities in t he p hoton and proton, of course. 

^Previous analyses |l42| , [7^| , which had led to very optimistic conclusions, had used much milder accep- 
tance cuts on the leptons, or none at all. 
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resolved photon process gq — > l + l~q contribute in LO to the production of a lepton pair with 
sizable pr- The corresponding ep cross-sections are both 0(al m ). 

This distinction has been well understood in existing treatments of the "Drell-Yan" pro- 
duction of pairs via the exchange of a virtual photon; see ref.[ |141|| for a LO estimate, 
and [107, 152 , 153 ] for NLO calculations. Compared to direct photon production this process 
has the theoretical advantage that one does not have to worry about fragmentation contri- 
butions. Moreover, backgrounds are low, and the final state can be reconstructed cleanly 
even at rather low Pt{1)i allowing one to probe quite small x values | |153| |. As usual, the 



region of positive rapidity corresponds to small x 7 and moderate x p , while negative rapidi- 
ties correspond to large x 1 and very small x p . The production of lepton pairs at positive 
rapidity and sizable pt is sensitive to the gluon content of the photon ||152 |. Unfortunately 
ii -i daiep — > l + l~X) . 

the expected event rates are quite low: ~ 1 pb/GeV at iVf/+;- = 4 GeV 

after integration over the transverse momentum of the pair and before any acceptance cuts 
have been imposed ||153 . 



In contrast, the fact that the total photoproduction cross-section for W and Z bosons is, 
in LO, given by the resolved photon contribution only has not been appreciated in the existing 
literature.0 The resolved photon contribution has first been estimated in ref. [|155|1 , but here 



this contribution was considered to be an addition to the direct process even at vanishing p? 
of the heavy gauge boson. This is not correct, since the direct contribution has a u— channel 
singularity which has to be absorbed into the resolved photon contribution; simply adding 
both contributions implies double-counting. While current calculations ||156| , |154|| of the 
total W and Z photoproduction cross-sections are in our view not entirely satisfactory, since 
they mix LO and NLO contributions in an ill-controlled manner, they should predict the 
production of high— px gauge bosons quite accurately.^] This high— px region is sensitive to 
the form of the W + W~ r y vertex [|156|| . The total cross-sections for W + and W~ production 
at HERA amount to approximately 0.5 pb each. This cross-section is to be divided by 
another factor of five if only the clean ev e and /xz/^ final states are observable; indeed, a 
first study [ |157|| has concluded that the detection of hadronically decaying W and Z bosons 
at HERA is quite challenging, although perhaps not impossible. Finally, we mention that 
recently the HI collaboration has announced [|158|| observation of one event that can be 
interpreted as the production of a leptonically decaying W boson at high p T ; within the SM 
this interpretation is quite unlikely, due to the smallness of the predicted cross-section, but 
all other interpretations seem even less likely. Clearly much more data have to be analyzed 
before any definite conclusion can be drawn. 



4) Real 77 Scattering at e + e Colliders 

In this section we discuss hard processes with two (quasi-)real photons in the initial state 
and a hadronic final state. At present, and in the near future, such reactions can only be 



*The total W and Z cross-sections in ep collisions also get contributions where the heavy gauge boson 
is radiated off the electron line; in case of Z production at HERA this contributes roughly 50% of the total 



cross-section [154 



Un principle the resolved photon contribution from gq — > Wq' should be included in a full LO treatment, 
but it is strongly suppressed at HERA energies due to phase space constraints. 
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studied at e + e~ storage rings, where the effective photon flux is given by eq.([TT]). At future 
linear e + e~ colliders the photon spectrum might receive large additional contributions from 
"beamstrahlung" [|1 59|| , which is emitted when a particle is accelerated in the field produced 
by the opposite bunch; beamstrahlung photons are exactly on-shell. The produced spectrum 
sensitively depends on the size and shape of the electron and positron bunches; see ref. ||160| 



for a handy parametrization of the spectrum in terms of a few machine parameters. Finally, 
in recent years the possibility has been discussed to convert a linear e + e~ collider into a "77 
collider" by scattering laser photons off the incident beams |161| . The achievable luminos- 



ity is predicted to be comparable to the (geometrical) e + e luminosity prior to conversion of 

the beams; in contrast, at existing e + e~ storage rings one has £ 77 ~ {^f^ In ~ lCT 3 £ ee . 
The photon spectrum of such a photon collider again depends on the geometrical set-up, 
and also on the polarization of the incident electron and laser photon beams. Most of the re- 
sults presented here will be for present and near-future colliders; these are obviously of more 
immediate interest, and already allow to illustrate most physics principles. For comparison 
we will also give a few results for the more "futuristic" colliders. 

Since we now have two photons in the initial state, we have to distinguish three physically 
distinct event classes f42|, see Fig. 20. Following the notation of ref. [|162|| we call reactions 



where the entire energy of both photons goes into the hard subprocess "direct" , see Fig. 20a. 
If only one of the photons couples in a pointlike manner while the other participates via its 
quark and gluon content, as in Fig. 20b, the process is called "once resolved" ("1-res" for 
short). Finally, processes where both photons are resolved into their partonic constituents are 
called "twice resolved" ("2-res"); an example is shown in Fig. 20c. Recall that each resolved 
photon produces a remnant or "spectator" jet, which goes approximately in the direction 
of the incident e beams; the three event classes are therefore characterized by having zero, 
one or two of these photonic remnant jets. Since the parton distribution functions in the 
photon are 0(a cm /a s ), all three contributions are of the same order in coupling constants, 



and have to be treated on the same footing |42 



We saw in the previous section that direct and resolved photon contributions to photo- 
production processes mix in NLO QCD. Similarly, the three event classes contributing to real 
77 scattering mix once higher-order QCD corrections are included. Parts of the NLO direct 
(1-res) contributions have already been included in the LO 1-res (2-res) terms; these parts 
therefore have to be subtracted from the NLO contributions. Note that mixing between 
direct and twice resolved contributions only occurs in NNLO in QCD. 

Following our preceding discussion of 7£> scattering, we discuss different final states in 
separate sub-sections. Jet production is treated in Sec. 4a, open heavy flavour production 
in Sec. 4b, J/ip production in Sec. 4c, and direct photon production in Sec. 4d. 

4a) Jet Production in 77 Collisions 

As in case of 77? scattering, the production of jets offers the largest cross-section of all hard 
77 collisions that lead to hadronic final states. In LO, the cross-section can be written as 
[see eq.([U|)]: 

dM T 7 = e M^m ^'-f^ ■ (24) 

dp T drjidrj 2 .j^, dt 
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If i or j is a quark or gluon, fi\ e can again be obtained by convoluting the photon flux 
function / T | e with the quark or gluon density in the photon, see eq.(|17l). However, antitag 
requirements at current e + e~ collider experiments often allow larger photon virtualities than 
one has in photoproduction events at HERA; at the same time, the scale \t\ or p\ of the hard 
process is usually smaller. This means that the suppression of resolved photon contributions 
due to the reduced parton content of virtual photons is usually larger at e + e~ colliders than 
at HERA. In our numerical estimates to be presented below we have included this effect 
using the formalism of ref . [95] . 



The scaling variables in eq. (|24]) are related to the jet (pseudo)rapidities 771,2 by: 

Xl = ^( e ^ + e" 2 ); (25a) 

x 2 = y (e- m + e"" 2 ) , (25b) 

with xt = %Pt/ \/s as before. The Mandelstam variable s of the hard scattering sub-process 
is again given by s = Xix 2 s, and t is as in eq. (|19"D . Most of the sub-processes contributing 
to jet production in 77 collisions also contribute to jet production at HERA. In particular, 
if both i and j are partons (twice-resolved contribution), a is the same as in resolved 7/9 
collisions, and the case where either i or j is a parton while the other is a photon (single- 
resolved contribution) corresponds to the direct contribution at HERA. Finally, if i — j — 7 
(direct contribution), then k = q, I = q and one has 

da{n r qq) = 4^ ft *\ 

where e q is the electric charge of quark q in units of the proton charge, and the factor of 3 
is due to colour. 

The first quantitative estimate of high— p T jet production in (quasi-)real 77 scattering has 



been presented by Brodsky et al. [163]. However, gluon-initiated processes were omitted, 
and a rather crude parametrization for the quark densities in the photon was used.Q Gluon- 
initiated contributions have been included in refs. ||165f , but again very simple parametriza- 



tions for the parton densities in the photon were used. In spite of their shortcomings, these 
early studies clearly demonstrated that resolved photon contributions are quite important, 
and often even dominant, if xt < 0.2. Notice that jet production from two-photon collisions 
at present e + e~ colliders cannot be studied experimentally if xt > 0.4; the cross-section 
becomes too small, and annihilation backgrounds too large. Hence resolved photon contri- 
butions were predicted to be sizable for at least half the experimentally accessible range of 
x T . 

Nevertheless data [T5J on multi-hadron production from 77 collisions taken at the PEP 



and PETRA storage rings were usually only compared to the direct (QPM) contribution.^ 
Not surprisingly, a significant excess of data over MC prediction was observed. Given that the 



*Rcf. [163 also finds quite large higher twist contributions to the production of hi gh— pt mesons and 



jets. It was shown later [ 164 1 that the "constituent interchange model" used in ref. |163J over-estimates the 



normalization of these terms by as much as a factor of a thousand. 



'''The PLUTO collaboration |16£] attempted to include 1-res qg or 2-res qq final states in their analysis. 



However, they made many simplifying assumptions, some of which are incorrect. In particular, they assumed 
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importance of resolved photon contributions had been emphasized quite early on [ |163| , |165| . 
it is surprising that it took ten years before the origin of this excess of data over QPM 
prediction was clarified experimentally. 

In ref . [ |162|| it was pointed out that the characteristics of these excess events at least qual- 
itatively agreed with those expected from resolved photon contributions: Their spectrum 
is softer than that of the direct contribution, just as in case of photoproduction events. In 
addition, resolved photon events are less two-jet like, i.e. have smaller thrust T, than direct 
events; in fact, the PLUTO collaboration had shown [ |166[ | that the excess could be removed 
by requiring T > 0.9. Ref. [|162|| also contains the prediction that, for fixed pr, resolved 
photon events should rapidly become more important as the beam energy is increased. The 
reason is that raising yfs decreases xt, and hence x^2 in eqs. (f25[) ; due to the additional 
convolution with photonic parton distribution functions, the quark and gluon density in the 
electron is always considerably softer, i.e. grows faster with decreasing x, than the photon 
flux in the electron. Resolved photon contributions were therefore expected to play an even 
more important role at TRISTAN and LEP than at PEP and PETRA. This was confirmed 
soon afterwards by the AMY collaboration || , who showed that their data on multi-hadron 
production in antitagged 77 collisions could be reproduced by their QCD-based Monte Carlo 
program, in both shape and normalization, if and only if it included the full set of resolved 
photon contributions, including those with a gluon in the initial state. This was the first 
unambiguous observation of resolved photon interactions other than in deep-inelastic e7 
scattering, and the first direct experimental evidence for a nonzero gluon density in the 
photon. 

Early experimental analyses of multi-hadron production in 77 collisions , including 
the AMY analysis p) , were not amenable to direct comparison with theoretical calculations, 
since the experiments defined a "jet" as a thrust hemisphere. The production of actual, 
reconstructed jets (using a cone algorithm) has been studied only quite recently, first by 
the TOPAZ collaboration | 167 | and then by AMY [ 168|| . This is an important development, 



since it allows a much more direct comparison between theory and experiment. 

In the meantime theoretical estimates are also becoming more sophisticated, by including 
NLO corrections. This process was actually already startd in 1979 with two calculations |169| 
of 77 — > gg via a quark box diagram. NLO corrections to the direct process 77 — > qq (for 



massless quarks) were calculated soon afterwards [17C]. However, the result contains collinear 



divergencies, which have to be absorbed in the 1-res contribution; in NLO it therefore makes 



little sense to consider the direct contribution in isolation. This was recognized in ref. ||171 
which deals with the production of high— p? hadrons. Here the direct contribution was 
treated in NLO (including the contribution from 77 — > gg, which is formally NNLO), but 
at the time resolved photon contributions could only be included in leading order.| This 
shortcoming could be remedied only after the NLO corrections to the hard partonic sub- 
process had been calculated HS5fl . A full NLO calculation of the single-jet inclusive cross- 



that the cross-sections for these resolved photon contributions drops like the square of the 77 invariant mass 
W; in contrast, for fixed pr QCD predicts these cross-sections to increase with W. Moreover, PLUTO did 
not attempt to include 1-res and 2-res contributions simultaneously. 

■^Notice that one only needs the LO 1-res contribution from jq scattering in order to absorb the collinear 



divergencies of the NLO direct term. The procedure of ref. 1 171 should therefore be quite accurate at high 
xt, where the direct contribution dominates. 
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section has become available only recently []172|| ; it finds only modest NLO corrections to the 
total cross-section (< 25%, for a jet cone size AR = 1 used by both TOPAZ and AMY, and 
renormalization and factorization scale /z = pr), but the relative weights of direct, single- 
and double-resolved contributions can change by larger amounts. f\ Finally, in ref . ||173|| the 2- 



res contribution has been studied in NLO, with emphasis on constraining the gluon content 
of the photon. 

Since NLO corrections appear to be quite modest, and since many more LO parametriza- 
tions of parton densities in the photon exist, we will only present results from LO analyses 
here. As mentioned earlier, we do include the suppression of resolved photon contributions 
due to the finite virtuality of the photons, using the simplest ansatze in ref.[^]. (Quark and 
gluon densities have to be treated separately.) We remind the reader that, unlike ref. ||172 



these ansatze do not assume a power-like suppression of the "hadronic" contribution to the 
parton densities of virtual photons; they might therefore slightly under-estimate the (in any 
case only modest) size of the correction. 

In Fig. 21 we compare our LO calculation of the single-jet inclusive cross-sections as 
measured by the TOPAZ (a) and AMY (b) collaborations at TRISTAN (^i = 58 GeV). 
TOPAZ requires the jet to have pseudorapidity |^j e t| < 0.7, and rejects an event if it contains 
an electron or positron at an angle 6 > 3.2° relative to the beam pipe and with energy 
E > 0.25-Ebeam = 7.2 GeV. AMY accepts jets out to |?7jet| = 1-0; moreover, it can antitag an 
event only if it contains an electron or positron with E > 0.25E^ cam at an angle 9 > 14.1°. 

We see from Fig. 21 that the DG, LAC1 and WHIT1 parametrizations reproduce the 
TOPAZ data at < 4.5 GeV quite well, and also fit the AMY data over the entire px 
range. Notice, however, that we had to use Nf = 4 active flavours already at PrQet) = 
2.5 GeV in order to achieve this agreement. This is not really justifiable for the DG and 
LAC parametrizations, which treat the charm quark as massless; on the other hand, WHIT 
explicitly includes some charm mass effects, which leads to a greatly reduced contribution 
from charm in the photon, in agreement with expectations. (It makes up for the shortfall 
by slightly larger light quark and gluon densities.) The GRV parametrization seems to fall 
below the data at low pt\ this is mostly because we have used Aqcd(A/ = 3) = 0.4 GeV 
for DG and WHIT, as compared to 0.2 GeV for GRV and LAC, as is implicit in these 



parametrizations. NLO corrections have been found to be most important at low px ||172|| ; 
they might very well bring the GRV prediction into agreement with the data. In contrast, 
the WHIT4 prediction lies above the data at low pr, since it assumes a quite large and 
hard gluon density in the photon. This discrepancy might be reduced if the "hadronic" 
piece of the parton density of virtual photons is suppressed by a power of the virtuality, but 
it seems unlikely that this will restore full agreement with the data. All the other WHIT 
parametrizations seem acceptable, however. 

Fig. 21 shows that the LAC3 parametrization clearly over-estimates the cross-section, 
while the direct contribution by itself falls well below the data. Notice that the resolved 
photon contributions remain non-negligible out to the highest pt where data exist, although 
their importance clearly diminishes with increasing p?, as discussed above. Finally, we 



§ However, there seems to be a discrepancy between the LO result of rcf. [172 and our calculation. In 
particular, Aurenche et al. find that 2-res processes still contribute about 25% to the jet cross-section 
measured by TOPAZ at pr = 7.5 GeV, while in our calculation it amounts to less than 10% at this rather 
large value of xt (unless we use the LAC3 parametrization). 
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observe that our calculation falls somewhat below the TOPAZ data for px > 4.5 GeV. This 
is disturbing, since in principle our neglect of the charm mass should be better justified if 
p T S> m? c . Notice, however, that our calculation agrees very nicely with the AMY data. This 
hints at an experimental discrepancy between these two data sets. Indeed, from the increased 
acceptance in rapidity as well as the much looser antitag requirements used by AMY, one 
expects that the cross-section measured by AMY should exceed that measured by TOPAZ 



by a factor of about 1.7 to 1.8; here we have made use of the fact |p.67| , |168|| that da/drjj et 
is quite flat over the observed range, and have used eq.(|TTJ) to estimate the effect of the 
tagging criteria on the photon flux. On the other hand, the cross-section measured by AMY 
only exceeds that reported by TOPAZ by a factor of 1.36 ± 0.26 (1.29 ± 0.20) at p T = 4.55 
(6.5) GeV, where we have assumed all errors to be independent. Although this procedure 
probably over-estimates the true error, since part of the systematics should be common 
to both experiments, this comparison indicates that the discrepancy in any one px bin is 
not very significant statistically. Due to bin-to-bin correlations produced by the unfolding 
procedure, it is difficult to combine different bins to arrive at an overall significance of the 
discrepancy. At the lowest px bin, the ratio of the two experimental cross-sections is 1.55, 
which agrees very well with expectations.^] Moreover, both collaborations also publish di-jet 
cross-sections, although with somewhat larger relative errors; our LO calculation reproduces 
both these measurements over the entire px range. 

Note that the published cross-sections are actually partonic cross-sections, which have 
been obtained by unfolding the observed pr(jet) and ^(jet) distributions using a LO MC 
program. We will argue in Sec. 5 that such a procedure might be more reliable at TRISTAN 
than at HERA, since events with multiple partonic interactions should not pose much of a 
problem here. Nevertheless this procedure does produce some dependence on the details of 
the Monte Carlo program, e.g. via the predicted jet reconstruction efficiency. It might also 
be dangerous to directly compare these extracted cross-sections with NLO calculations, since 



the higher order corrections seem to change [172 the relative weights of the three classes 



of contributions compared to what has been assumed in the LO MGQ Both experiments 
also give the measured rate of jet events as a function of px\ however, at these rather low 
transverse momenta the jet reconstruction efficiency is still rather small and px~ dependent. 
A fair amount of MC work is therefore (unfortunately) necessary before a comparison with 
theoretical predictions can be made. 

The published data are based on an integrated luminosity of about 90 pb _1 for TOPAZ, 
and 27 pb _1 for AMY. The total available data samples are more than three times larger than 
this. Notice that the error bars in Fig. 21 are already quite small; if the slight discrepancy 
between AMY and TOPAZ can be resolved, the full data set should therefore allow a quite 
precise comparison between theory and experiment. 



^The effect of the looser antitag requirement of AMY is smaller at low pr, due to the dynamical upper 
bound P 2 < pj, that has to be imposed on the virtuality of the exchanged photons in order to meaningfully 
speak of parton or photon densities in the electron. 

*It is not obvious how the "theoretical" definition of the three classes of contributions used in ref. |172[ 
relates to a more experimental definition relevant for event reconstruction, which could e.g. be based on 
the presence of absence of (remnant) jets close to the beam pipes. For example, the relative weights of the 
three classes of contributions in the NLO calculation strongly depend on the factorization scale, which is 
an unphysical parameter. Unfortunately it is quite difficult to write an MC generator based on a full NLO 
calculation; no such generator for jet events produced in hadronic or photonic collisions exists as yet. 
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The larger data sets should also allow to investigate more differential distributions, using 
events with two reconstructed jets. As an example we show in Fig. 22 the triple-differential 
cross-section for events with two jets with equal (pseudo)rapidity 771 = 772 = 77. This distri- 
bution is quite flat for the direct contribution, at least over the range covered by TRISTAN 
detectors. For parametrizations with a relatively hard gluon density (WHIT1, DG) the total 
single-resolved contribution also depends only weakly on 77; for the given choice of p?, the 
decrease of the contributions oc q] {qg final state) is more or less balanced by the increase of 
those oc G 1 (qq final state). Increasing r\ means increasing x\ but decreasing x 2 in eqs.(^3f). 
Since the gluon content in the electron is always quite strongly peaked at small x while the 
photon flux function / 7 | e is relatively hard, the increase of fa\ e more than compensates the 
decrease of / 7 | e as 77 is increased. (Note that f q \ e increases much more slowly with decreasing 
x than fa\ e does.) If the gluon density in the photon is large but strongly peaked at small 
x (WHIT6, LAC1), the rise of the cross-section for 75 — > qq even leads to a total 1-res 
contribution that peaks at 77 ~ 2, rather than at 77 = 0. 

Turning to the twice-resolved contributions shown in Fig. 22b, we observe that the cross- 
sections for processes with two gluons in the final state, which most of the time also have two 
gluons in the initial state, is quite strongly peaked at 77 = 0, even for the comparatively hard 
gluon density of WHIT1. The same parametrization predicts the cross-section from qg — ► qg 
scattering to remain quite flat for \r)\ < 1, since the decrease of f q \ e (xi) is compensated by the 
increase of /G|e(^2)-Q Nevertheless, parametrizations with relatively hard gluon density still 
predict a significant decrease of the total 2-res cross-section as 77 is increased from to 1. In 
contrast, parametrizations with a large and soft gluon content (WHIT6, LAC1) predict the 
2-res contribution to rise with increasing 77, reaching a maximum at 77 ~ 1.3. This increase 
is entirely due to the qg final state, where the increase of fa\ e now over-compensates the 
decrease of f q \ e . 

As discussed in the next subsection, the TOPAZ collaboration has proven capable of 
detecting photon remnant jets with an efficiency of about 70%; this ability was already 
implicit in their observation [ I57 | of considerable energy flow at small angles relative to the 
beam pipes. They are now beginning to exploit this capability also in the jet analysis ||174|| . 
This shows that the 1-res, 2-res and direct contributions can indeed be studied separately. 
On the other hand, Fig. 22 shows that even the sum over all contributions should allow to 
discriminate between some existing parametrizations of the parton content of the photon. 
In particular, at px = 3 GeV, parametrizations with a large but soft gluon density predict 
a flat or even slowly rising jet cross-section as 77 is increased away from zero, in contrast to 
parametrizations with a small or hard gluon density. Gluon-initiated processes are expected 
to contribute more than 50% of the total 2-res cross-section, and still some 25 to 35% of the 
total di-jet cross-section at small rapidity. This fraction, and hence the difference between 
predictions using different parametrizations, will be smaller (larger) for px > (<)3 GeV. 

Experiments at the LEP storage ring should also be able to contribute significantly to 
our understanding of (almost) real 77 collisions. At LEP1 the vicinity of the Z peak means 
that multi-jet annihilation events are a much more severe background than at TRISTAN; 
however, a parton-level investigation [ |175( | concluded that this should not be much of a 
problem as long as the invariant mass Mjj of the two high— px jet system is below 15 

'Of course, there is also a contribution oc f q \ e (z2)fG\e( x i)) but it decreases very quickly with increasing, 
positive 77. 



41 



to 20 GeV. Indeed, some experimental results have already been published. The ALEPH 
collaboration repeated [|176|| the AMY analysis || of multi-hadron production. Although 
some of the details differ, as will be discussed in Sec. 5, ALEPH also finds that the description 
of the data is greatly improved once resolved photon interactions are included in the Monte 
Carlo model. 

DELPHI [177 reached similar conclusions regarding "minimum bias" multi-hadron pro- 
duction. In addition, they analyzed a subsample of events containing two jets with \r]\ < 1.0 
and pt > 1-75 GeV, where the jets were defined using a cluster algorithm. At these small 
transverse momenta the jet reconstruction efficiency is rather low [<t(2— jet) ~ O.ler(parton)], 
which introduces a strong dependence on the details of the MC program. DELPHI finds that 
DG falls below the data for pt < 3 GeV, while the LAC1 and GS parametrizations work 
well. Note that their event sample is mostly sensitive to direct events, as well as resolved 
photon events with large x, i.e. soft remnant jets, since they require W vis < 13 GeV and 
E v[s < 20 GeV; the calculation of these quantities includes information from the small-angle 
taggers, which cover angles down to 2.5° and should therefore see parts of the remnant jets. 
Fig. 3 shows that the DG parametrization does indeed fall well below LAC1 for large x and 
small Q 2 . 

More recently, DELPHI has published ||1 78|| a study of hadronic 77 events where either 
the e + or the e~ is tagged at a very small angle, corresponding to a photon virtuality 
P 2 ~ 0.1 GeV 2 . Since even stronger cuts were imposed on the hadronic system than in the 
first DELPHI analysis, the data sample was quite small (491 events); nevertheless it was 
sufficient to once again prove the necessity to include resolved photon contributions if the 
data are to be described by the MC program. This small angle tagging technique holds 
much promise, since the tag helps to reconstruct the kinematics of 77 events by determining 
the energy of one of the two photons. If (as at HERA) the energy of the second photon 
can be determined using the Jacquet-Blondel method, a measurement of the jet rapidities 
in di-jet events would allow to directly reconstruct the Bjorken— x variables of the partons 
in the photons. If, as TOPAZ results indicate, the DELPHI small angle detectors (not to 
be confused with the very small angle taggers used for the electron tagging) can be used for 
detecting the presence or absence of remnant jets, DELPHI (and similar detectors) should 
be able to perform quite detailed analyses of hard two-photon reactions. In order to fully 
exploit this potential, the strong upper limits on the visible energy and invariant mass used 
in refs. [ |177| , |178|| should be relaxed; this should certainly be possible at LEP2 energies, where 
the annihilation background is much smaller. 

The LEP energy is expected to soon be increased to y/s ~ 180 GeV; LEP experiments will 
then have the unique opportunity to study two-photon cross-sections over a wide range of 
energies. As illustrated in Fig. 23, the energy dependence can be quite strong. In this figure 
we show the triple-differential cross-section at p T = 5 GeV as a function of r)\ — r) 2 = r], 
imposing an antitag condition similar to that used by ALEPH ||1 76| . For ^/s = 90 GeV, xt 
is similar to the value used in Fig. 22, which leads to a similar shape of the pseudorapidity 
distribution; the kinks at rj ~ 1.5 occur because the ALEPH antitag becomes ineffective 
if the outgoing electron has less than half the beam energy, i.e for scaled photon energy 
z > 0.5. Raising ^ from 90 to 180 GeV increases the direct contribution at 77 = only 
by a factor of 1.45, while the 1-res and 2-res contributions increase by factors of 2.0 and 
3.2, respectively. This again demonstrates the strong dependence of the cross-section for 
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resolved photon processes on the available phase space, or, equivalently, on the Bjorken— x 
variables of eqs. (|25|) ■ Note that the single-jet inclusive cross-section at central rapdity grows 
even faster with energy, since the kinematical integration limits for the rapidity of the second 
jet also increase with y/s. Finally, we remind the reader that the shape and normalization 
of the solid curves in Fig. 22 could be quite different if the photon has a large but soft gluon 
component; this would also result in an even more rapid increase of the cross-section for 
resolved photon events with y/s. 

LEP will almost certainly be the highest energy e + e~ storage ring ever. In order to 
reach even higher energies, one will need to build a linear collider. In such a device each 
bunch can most likely only be used once; the repitition rate (number of bunch collisions 
per second) will therefore almost certainly be smaller than at LEP. At the same time the 
total luminosity must increase like the square of the beam energy in order to maintain a 
roughly constant rate of e + e~ annihilation events. The luminosity per bunch crossing will 
therefore have to be much larger than at existing storage rings, forcing one to use very dense 
bunches. The correspondingly large charge density gives rise to strong electromagnetic fields. 
When particles inside one bunch enter the field produced by the opposite bunch, they will 
be accelerated, and will therefore radiate real photons. This is known as "beamstrahlung" 

The flux of these beamstrahlung photons depends quite sensitively on the design char- 
acteristics of the collider; this is not surprising, since this radiation is due to the field of the 
entire bunch, not due to individual e + e~ collisions (unlike the only quasi-real "bremsstrahl- 
ung" photons we have dealt with up to now). Beamstrahlung can be reduced by using flat 
beams, and by splitting large bunches into "trains" of smaller ones. On the other hand, for 
a given class of designs, beamstrahlung increases quite rapidly with increasing beam energy. 

This is ilustrated in Figs. 24, which show the single-jet inclusive cross-section as a 
function of pr, where we have accepted jets with \rj\ < 2.0. We have used the WHIT1 
parametrization, and imposed the antitag requirement 9 < 10° when computing the brems- 
strahlung contribution to the total photon flux. The beamstrahlung spectrum has been 
calculated using the analytical expressions of ref. [|160|| , for the JLC design as specified at 
the 1993 international linear collider conference | 179|| . Designs for linear colliders are still 



evolving; the results of Fig. 24 should therefore be considered as indicative only. 

At \fs =0.5 TeV (Fig. 24a) the beamstrahlung spectrum is considerably softer than 
the equivalent bremsstrahlung spectrum (jTT|). This enhances the relative importance of the 
direct contribution, since the cross-section for resolved photon contributions increases with 
the two-photon invariant mass W while that for the direct contribution decreases. For 
Pt > 100 GeV, the total cross-section is dominated by directly interacting bremsstrahlung 
photons. Notice that this design leads to a luminosity of about 50 fb _1 per year; the two- 
photon cross-section should therefore remain measureable for jets with p? well above 100 
GeV. Moreover, one expects of the order 10 8 events per year with a jet with > 5 GeV. 
This sounds like a large number, but corresponds to a trigger rate of 10 Hz or less, which 
should be easily manageable. 

Increasing \fs to 1.0 TeV (Fig. 24b) greatly increases the flux of beamstrahlung photons, 
and also makes it harder. This enhances the relative importance of resolved photon contri- 
butions, as can be seen from the xt value where direct and resolved photon cross-sections 
are equal. Notice also that now beamstrahlung increases the cross-section by about a fac- 
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tor of 5.5 even at px = 200 GeV. Finally, comparing Figs. 24a and b, we see that the jet 
cross-section increases by a factor of about 7 (20) for px = 5 (100) GeV; without beamstrahl- 
ung, the corresponding factors would "only" have been 3 and 5, respectively. We remind 
the reader that these results depend on the specific machine design; see refs. ||180| , |181|| for 
further discussions of this point. 

As mentioned in the beginning of this section, one might be able to convert a linear e + e~ 



collider into a 77 collider by back-scattering laser photons off the e + and e~ beams ||161 |. 
The spectrum and luminosity of such a collider depend quite strongly on details such as the 
polarization of the laser photons and incident electrons beams. The production of jets at 



such a collider has first been discussed in ref. [|182| . In Fig. 25 we present a result for the 



simple case of unpolarized beams and small distance between the conversion and interaction 
points. This leads to a photon spectrum which peaks at z = 0.828, where it also cuts off; the 
two-photon luminosity is then quite flat over a wide range of W. Most modifications that 
have been discussed in the literature ||161|| give even harder photon spectra. By comparing 
Fig. 25 with Fig. 24a we see that, as expected, the much harder photon spectrum of the 
77 collider has greatly increased the relative importance of resolved photon contributions. 
Moreover, the hard photon spectrum allows to efficiently access soft partons in the photon 
via 1-res contributions at large rapidity. As a result, the cross-section increases with 7] even 
for parametrizations with rather modest gluon content, like WHIT1; this is to be contrasted 
with the situation at present and future e + e~ colliders, see Figs. 22 and 23. This also explains 
the large difference between the predictions from the LAC1 and WHIT1 parametrizations 
close to the kinematical maximum of rj, inspite of the rather large value of xt- 

Finally, we should warn the reader that from LEP2 energies onwards, multiple interac- 
tions could substantially increase the true jet cross-section, compared to the simple parton- 
level estimates of Figs. 23 to 25; we already saw in Sec. 3a that this phenomenon seems 
to play an important role in jet production from resolved photons at HERA. This will be 
discussed in more detail in Sec. 5. 



4b) Heavy Quark Production in 77 Collisions 

Apart from the production of jets discussed in the previous subsection, the production of 
heavy quarks is the only hard QCD process that has been studied experimentally in two- 
photon collisions. As discussed in Sec. 3b, the main advantage of heavy QQ pair production 
is that perturbative QCD should be applicable over the entire phase space. However, as 
we already saw in the discussion of the photoproduction of heavy quarks, predictions for 
total cc production rates are at present still quite uncertain, due to unknown higher order 
corrections (which lead to a strong scale dependence) and the uncertainty in the value of m c 
to be used here. 

Progress in the theoretical treatment of heavy quark production in two-photon processes 
has been relatively slow. The first complete LO calculation, including resolved photon pro- 
cesses, was only performed in 1989 ||162fl . The contribution from twice resolved processes was 
found to be very small at TRISTAN energies, but the 1-res contribution (from 7$ — > cc) is 
quite sizable. Further LO predictions, for newer parton densities in the photon, were pub- 



lished in [ 4H|| . A full NLO analysis of the direct and 1-res contributions has been performed 
in ref. [|183j| . 
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In Fig. 26 we show updated predictions for the total cc cross-section in the PETRA to 
LEP2 energy range. We have included NLO corrections to the direct contribution, using 
a simple parametrization|] In ref. [183] the direct contribution to the total cross-section is 



written as 

a dir ( 7 7 - QQ{g)) = ^ (4? + 4™ sC M) . (27) 
Here c^ describes the well-known [|184|1 tree-level (QPM) prediction, while cj^ can be 



77 

parametrized as: 

c - = |-^(^-i)' r<2 - 637 

= 0.35r~ a3 , r > 2.637 (28) 

where r = W^/lAmq) — 1. This parametrization is exact at threshold r — > 0, and describes 
the full NLO result to better than 10% accuracy for all r < 100. On the other hand, the 
predictions for the 1-res contribution shown in Fig. 26 have been computed in LO only. One 
reason is that here NLO corrections are considerably smaller than for the direct contribution 
183|| , largely because the threshold ( "Sommerfeld" ) corrections are negative for a colour- 



octet QQ state. Moreover, the uncertainty of the prediction is much larger than for the 
direct cross-section, making it less important to include 10% corrections. 

In Fig. 26 this uncertainty is given by the width of the bands defined by two curves with 
the same pattern. Following ref. ||162j| we have in all cases required W 77 > 2m D = 3.74 GeV, 



but the "dynamical" charm quark mass appearing in the expressions for the hard cross- 
sections is less well defined. In case of the direct cross-section we have varied m c between 
1.3 GeV (upper dotted curve) and 1.6 GeV (lower dotted curve). Another uncertainty arises 
from the choice of scale in a s and (for the 1-res contribution) in G 7 . In Fig. 26 we have 
estimated this uncertainty by varying this scale between M cg /4 (upper curve) to M c£ (lower 
curve). In case of the direct contribution the combined uncertainty only amounts to slightly 
over 20%, almost independently of \fs. Note, however, that other authors [|112| , |183|| prefer 
to use an even wider range of values for m c . 

Unfortunately the uncertainty for the prediction of the 1-res contribution is considerably 
larger than this. One reason is that now a s already appears in the tree-level cross-section, 
leading to a stronger dependence on the renormalization scale; there is also a factorization 
scale dependence in this case. As usual, the inclusion of NLO corrections should reduce 
these scale uncertainties. However, the biggest uncertainty comes from the choice of m c , and 
of the minimal allowed M c5 . In case of the direct contribution one always has W 77 = M c5 , 
at least in leading order. On the other hand, 1-res events have W 77 > M cg ; it is then not 
clear whether one has to require M c5 > 2m n in order to describe open charm production, 
or whether it is sufficient to have W 77 > 2m d- If M c5 < 2m d < W 77 , some energy has to 
be transferred from the remnant jet to the "hard" cc pair; it has to be remembered that 
in any case colour needs to be exchanged between these two systems in the hadronization 
step. It seems unlikely to us that this soft energy exchange can exceed 1 GeV, however. We 
have therefore used m c = 1.4 GeV, M c5 > 2m c and scale \i = M c5 /4: in order to estimate the 

*Note that here direct and 1-res contributions remain well-defined even in NLO, since there is no LO 
1-res contribution from the quarks in the photon. 
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upper limit of the 1-res contribution, and m c = 1.6 GeV, M c5 > 2m£> and /x = M cg for the 
lower limit. 

Fig. 26 shows that this uncertainty makes it impossible to distinguish between the DG 
and WHIT1 parametrizations based on the total cc production rate alone. Note that the 
uncertainty is larger for LAC1, which has a very soft gluon density and therefore reacts 
very sensitively to changes of the lower bound of M c5 . Even for the more conservative 
parametrizations, the uncertainty amounts to roughly a factor of two. Nevertheless, some 
sensitivity to the parton densities does remain; we will come back to this point shortly. 
Finally, we note that even at LEP2 the 2-res contribution amounts to at most 5% of the 
183 1; we have therefore not shown it in this figure. 



1-res one 



In the last few years several experimental studies of charm production in 77 collisions 
have been published. The first analysis, by the JADE collaboration [|185|1 , looked for fully 
reconstructed charged D* mesons in single-tag events. They reported a considerable excess 
over QPM predictions, but this has not been confirmed by the TASSO ||186|| or TPC/27 ||187 



collaborations, who took data at similar energies. TASSO also measured the cross section 
for exclusive D°D° production. The total cc production cross-sections derived by these two 
experiments agree with expectations [ 183 1 . 

The first study of charm production in two photon collisions at TRISTAN, by the TOPAZ 
collaboration ||188| , also searched for fully reconstructed D* mesons. Unfortunately the 
reconstruction efficiency is quite low, leading to rather poor statistics (a few dozen events per 
experiment). TOPAZ therefore also published a study |189| based on D^* — > nfD° decays, 
where only the soft pion needs to be detected; here the signal after background subtraction 
consists of 372 ± 54 events. VENUS ||190|| and TOPAZ [ |191|| have also published analyses of 
hadronic two-photon events containing an electron or positron (electron-inclusive analysis); 
after subtracting backgrounds (mostly from Dalitz decays and photon conversions), these 
experiments extract a charm signal with (9(100) events each. Finally, very recently ALEPH 
presented | 192 | results of an analysis based on 33 fully reconstructed D^* mesons. 

Unfortunately these experimental results are somewhat contradictory. All TRISTAN 
experiments find some excess of events with high p?] e.g., TOPAZ ||189|| reports a 2.9 a 
excess of events with pt(D*) > 3.6 GeV. This is actually not all that surprising, given that 
we had to include contributions from the charm in the photon in order to reproduce the 
production of central jets with p? > 2.5 GeV; such "charm excitation" contributions are not 
included in present MC programs.^] On the other hand, ALEPH ||192|| finds a cross-section 
for the production of charged D* mesons with px > 2.0 GeV that agrees with the lower 
range of predictions; no excess is visible here. 

Most of these experimental analyses are only sensitive to charmed hadrons with significant 



transverse momentum. The notable exception is the TOPAZ inclusive electron analysis [|191 



electrons with momentum as low as 400 MeV are accepted, so that even charm quarks at 
rest can contribute to the signal. Moreover, TOPAZ used their forward calorimeter, which 
covers angles down to 3.2° with respect to the beam pipe, to look for the presence of photon 
remnant jets in the event; their MC predicts a jet tagging efficiency of 73 ± 2%. This allows 
them to study direct and resolved photon contributions separately. The direct contribution 



''The NLO correction to the direct process should describe 7c — > gc scattering accurately at these energies; 
however, it is not clear whether the parametrized form of the NLO corrections used by TOPAZ [188, |18S[ 
treats such contributions properly. Moreover, 2-res excitation contributions are not included at all. 
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agrees roughly with the upper range of NLO predictions. The resolved photon contribution 
agrees with NLO predictions using the LAC1 parametrization, while the DG prediction 
seems to be at least two standard deviations below the data. This demonstrates that this 
kind of analysis can lead to significant constraints on the gluon density in the photon. 

In order to compare experimental results on charm production with theoretical calcu- 
lations, one has to model the quark to hadron transition. This is more difficult here than 
in the more familiar high-energy e + e~ annihilation events, since for the low values of 
or Pt{c) relevant for present two-photon data the effect of hard gluon radiation can not 
be absorbed into fragmentation functions; this is only possible if gluons are predominantly 
emitted collinear to the quarks, which requires px 3> m c . In this context it is interesting to 
note that the CLEO collaboration ||193|| , working at values of M cS only slightly above the 
upper range probed by present two-photon experiments, found that popular fragmentation 
functions could describe the observed D* meson spectrum only if the emission of hard gluons 
was allowed explicitly. In particular, it was not possible to describe the spectrum by simply 
convoluting the tree-level cc cross-section with the Peterson et al. fragmentation function 
10 .511 . Other fragmentation functions gave a better fit, but only if the fragmentation param- 



eters were chosen differently from those used for data at higher energy; this is not surprising, 
since fragmentation functions are scale dependent, just like parton distribution functions. It 
is therefore encouraging to note that an NLO event generator for heavy quark production 
in two photon collisions, written by J. Zunft, is now publicly available |[L94j| .f| Once hard 
gluon radiation has been included explicitly, the fragmentation function only has to include 
nonperturbative effects, which should indeed factorize to good approximation. Studies like 



that by the TOPAZ collaboration ||191|| , which are sensitive to charmed hadrons at rest, can 
play an important role in testing this formalism, since by definition fragmentation functions 
can only change the spectrum, but not the total cross-section. Once the fragmentation pro- 
cess has been fully understood, direct cc pair production might be the best way to measure 
the value of m c to be used in perturbative QCD calculations, since the scale uncertainty is 
rather small here. A good knowledge of m c would help to sharpen predictions for 1-res cc 
production, as well as photo- and hadro-production of charm. 

Fig. 26 shows that the total cc pair production cross-section is expected to grow quite 
rapidly with energy. However, at higher energies it might be necessary to impose quite 
stringent cuts in order to extract a charm signal; this will reduce the detectable cross- 
section significantly. For example, in ref. ||183|| it was found that requiring one of the charm 
quarks to have rapidity \y\ < 1.7 and pt > 5 GeV reduces the cross-section by almost a 
factor of 50. This still leaves us with at least 5,000 cc events in 500 pb _1 of data; however, 
at this point we have not yet required anything that would actually identify these events as 
being due to charm, e.g. a hard lepton or a reconstructed D* meson. 

At a 500 GeV e + e~ linear collider even the DG parametrization predicts |[L81| , |195|| the 
total cc pair cross-section to reach a value between 5 and 50 nb, depending on the amount 
of beamstrahlung generated. However, requiring the event to contain at least one muon with 
rapidity \y\ < 2 and Pt(aO > 5 GeV reduces this cross-section by at least a factor of 2,000. 

As in case of jet production, cross-sections for the production of heavy quarks can be 
boosted considerably by converting an e + e~ collider into a 77 collider by means of laser 

■'This generator only allows the emission of a single hard gluon; the generator used by CLEO to model 
e + e~ annihilation events included 2^4 matrix elements. 
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backscattering. However, at such a collider QCD processes will likely be regarded primarily 
as backgrounds to "new physics" signals. In particular, one advantage of a 77 collider 
is that it can produce a Higgs boson as an s— channel resonance. If m# < 150 GeV, 
this Higgs boson will predominantly decay into into a bb pair [ 196 . The direct (tree-level) 



background from 77 — > bb can be reduced significantly by chosing the two incident photons to 
be (predominantly) in a J z = state [197|. There is still a direct background from 77 — > bbg, 
but this appears to be manageable [198] - However, at least if one uses a collider with a broad 
photon spectrum, which allows to search a wide range of Higgs masses simultaneously, the 
dominant background will actually come from resolved photon events. This has been pointed 



out in ref . [|199|1 , where the 1-res — > bb contribution is studied. As emphasized in ref . [|200 



the exact size of this background will also depend on the extent to which a gluon in a 
polarized photon is itself polarized. 

The to date most complete study of resolved photon backgrounds to the Higgs signal at a 
broad-band 77 collider has recently been performed by Baillargeon et al. ||201|| . They ignore 
polarization effects on the parton densities in the photon (which amounts to the assumption 
AG 7 = 0), but include all possible processes, including b and c excitation contributions; 
these are very large in the relevant kinematical region. They include production of c quarks, 
as well as light partons, in order to estimate the effects of imperfect 6— tagging. 

An example is shown in Fig. 27, which shows the signal for a 120 GeV Standard Model 
Higgs boson [ |196|| (assuming a bb invariant mass resolution of 5 GeV), as well as various 
backgrounds. The electron beams have 175 GeV energy, and 90% polarization; 100% polar- 
ized laser beams have been assumed. This suppresses the direct backgrounds by about an 
order of magnitude. Moreover, a px cut of 30 GeV has been applied. We see that the largest 
contribution to the single— b inclusive cross-section comes from 2-res b excitation processes 
(labelled as "&X 2 _ res "). Depending on details of the b— tagging efficiency, it is therefroe often 
advantageous to require both high— p T jets to be tagged as b quarks, or at least as heavy 
flavours. 

Another conclusion of ref. | |201[ | is that Higgs searches at a 77 collider become easier with 
increasing Higgs mass and decreasing beam energy. The reason is that, as we have seen 
several times already, cross-sections for resolved photon processes increase rapidly with y/s, 
but decrease equally quickly when the invariant mass of the hard system is increased for fixed 
yfs. In particular, for the case shown in Fig. 27, even with a b— tagging efficiency comparable 
to present LEP detectors, a 7 o Higgs signal could be extracted from 10 fb _1 of data, which 
roughly corresponds to one year's running; for the same parameters, one could at best hope 
to extract a 3.5 a signal at a 500 GeV collider. Since one here probes the photon at very high 
momentum scales these conclusions do not depend very sensitively on the parametrization 
chosen ||201 |. Nevertheless, Fig. 27 clearly shows that a good understanding of heavy quark 
production in resolved photon interactions is mandatory if we ever want to look for Higgs 
bosons with mass below 150 GeV at a 77 collider. 



4c) J ftp Production in 77 Collisions 

J/ijj production in principle offers a clean method for constraining G 1 in two photon colli- 
sions ||162|| . At least in the framework of the colour singlet model |pL 35|| and to leading order 
in a s , only resolved photon processes can contribute to 77 — > J/^+hadrons. Moreover, the 
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2-res contribution is expected to be very small, just as in case of open heavy flavour pro- 
duction discussed in the previous subsection. In LO the cross-section is therefore essentially 
proportional to G 7 ; the hard sub-process, 75 — > J /if) + g, is the same as in direct inelas- 
tic photoproduction of J /if) mesons, see Sec. 3d. As mentioned in that subsection, NLO 
corrections to J /if) production in 75 fusion have recently been computed ||139|| . We use the 
results of this calculation to update our previous LO predictions | |162| | for J /if) production 
at TRISTAN energies, and also extend the calculation into the LEP energy range. 

In order to simplify this calculation, we again use a parametrized form of the NLO 
corrections. According to ref. |p.39|| the total cross-section for 75 — > J/if> + X, integrated over 
the region Z < 0.9 where Z is the inelasticity defined in eq. (|23"l) , can be written as:[] 



a{ 19 ^ J/if) + X) 



a em a 2 Se 2 c \if)(0)\ 
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, (29) 



where r = s/m 2 j^ — 1. The wave function at the origin |?/>(0)| 2 can be derived from the 
measured leptonic partial width ||140| , |139|| . The function c^°\r) describes the tree-level 



cross-section [|135|j ; we do not impose any cut on Z for this part, since it remains finite in 
the limit Z — > 1, unlike the NLO correction, which diverges in this limit. These corrections 
are given by the two functions c^'(r), c^(r), which can be parametrized as: 
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Note that the NLO corrections (|30|) remain finite as r — > 00, while the tree-level cross-section 
drops like 1/r 2 in this limit. 

In Fig. 28 we show results of a partial NLO calculation of J/if> production at current 
e + e~ colliders based on eqs. (|29|) and (^). We did not include the small NLO contribution 



from 7^ scattering [1313]; moreover, the NLO direct contribution from 77 — > J /if) + gg is 
not yet available. Note also that if)' production with subsequent if)' — > J /if) decay should 
increase the total J /if) yield by another 15% or so. The results of Fig. 28 are for a no-tag 
situation, and include the suppression of G 1 due to the virtuality of the photon |p5" |. We 
checked that our parametrization of the NLO cross-section shows significantly reduced scale 
dependence, as expected, at least as long as the scale /i 2 > 1.5m 2 ; the cross-section drops 



quickly for even smaller \i. The same behaviour has been observed in ||139 | for the case of 
photoproduction of J /if). We therefore estimate the scale dependence by varying jj 2 between 
m 2 (lower curves) and 2m 2 (upper curves); the cross-section slowly decreases again for larger 
values of fi. In addition, we have varied m c in eq.(|29|) from 1.6 GeV (lower curves) to 1.3 
GeV (upper curves). Since we are now explicity producing a colour singlet state already in 



*In ref. 1391 the factor , 



has been written as l/m^.. However, in the framework of the colour singlet 

model it seems more natural to us to have the partonic cross-section scale like rather than l/m? c . In 

the latter case the uncertainty related to the value of m c discussed below would obviously be significantly 
larger. 
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the hard scattering process we have always required the 7^ cms energy to exceed mj/^ } even 
if 2m c < mj/f.f\ 

Fig. 28 shows that the uncertainty estimated in this way is still quite sizable; most of 
this uncertainty, however, is due to the variation of the factor l/m 2 c in eq. (|29|) . Notice also 
that the LO estimate (dotted curve, for scale /i 2 = 2m 2 c and m c = mj/^/2) is quite close to 
the lower range of NLO predictions. This indicates that the perturbative result is relatively 
stable. Recall also that ref. ||139|| finds good agreement of their prediction with low-energy 
photoproduction data, although leptoproduction data favour larger cross-sections. However, 
there might be sizable corrections to the non-relativistic approximation intrinsic to the colour 



singlet model ||138|| ; this has not been included in our estimate of the theoretical uncertainty. 

We see that the total expected cross-section for J / ip production at TRISTAN lies in the 
range from about 0.4 to 2.5 pb, depending on the parametrization of G 1 '. This corresponds 
to several hundred events per experiment in the total TRISTAN data sample. However, one 
will almost certainly have to demand that the J/ip meson decays into an e + e~ or pair; 
this reduces the rate by about a factor of 7.5 even before any acceptance cuts have been 
applied. The cross-section increases quite rapidly with energy, possibly exceeding 10 pb at 
LEP2; however, the acceptance at these higher energies might also be poorer. A detailed 
MC study is necessary to decide whether the study of J / ip production in two photon events 
at current colliders is feasible. 

As usual, the cross-section at future linear e + e" colliders is expected to be even larger, 
partly due to contributions from beamstrahlung photons. Since the invariant mass of the 
produced final state is quite small, the cross-section is now largest for designs giving a large 
number of rather soft photons. For example, at the 500 GeV TESLA collider (we again use 
the 1993 design ||179| 1) one expects a total J/ip production cross-section between 0.3 and 2.0 



fib, or some 10 7 events per year. At the same collider the T(ls) cross-section is expected to 
lie between 0.25 and 0.85 pb, which still gives around 10,000 events per year. It is at present 
not clear what fraction of these events would be detectable in the considerably "dirtier" 
environment of such linear colliders. 



4d) Direct Photon Production in 77 Scattering 

The production of a direct photon in the collision of two (quasi-)real photons, 77 — > 7 X, is 
another reaction that in LO only receives contributions from resolved photon processes: The 
direct 2^3 process 77 — > qqj is 0(a^ m ), and thus of next-to-leading order compared to the 
1-res contribution from 7g — > 75, as well as the 2-res contributions from gq — > ^q and qq — > 
#7, which are formally 0{a^ m / a s ); the relevant hard scattering cross-sections can be found in 
ref. ||202||. There are additional LO contributions involving parton — > photon fragmentation; 



as discussed in Sec. 3c, the corresponding fragmentation functions are 0(a em /a s ), which 
compensates for the relative factor of a s /a em that appears in the hard sub-process cross- 
sections for these contributions. However, these fragmentation contributions will in general 
lead to softer photons, which are accompanied by a jet. 



^In ref.[ |l39[ r has also been written as a function of m^, rather than a function of m 2 j^. However, this 
results in negative cross-sections for m c > It seems to us that the m c dependence of the c functions 

should be treated on a par with corrections to the non-relativistic approximation used in the colour singlet 
model. 
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The only existing calculation |202|| for this process omitted the fragmentation contribu- 



tions; this roughly corresponds to imposing stringent isolation requirements for the direct 
photon. Ref . ||202| | uses LO expressions; an NLO calculation is now under way [54 . 

In Fig. 29 we show updated LO estimates^. We have required > 1.5 GeV, since QCD 
perturbation theory would be very unreliable for even softer photons. We have also applied 
the acceptance cut |r/ 7 | < 1 on the (pseudo) rapidity of the emitted photon. In case of the 
WHIT1 parametrization, 2-res and 1-res contributions are shown separately, together with 
the sum; for all other parametrizations only the total prediction is given. We see that WHIT1 
predicts the 2-res contribution to be very small in the PEP to TRISTAN energy range, and 
to remain sub-dominant even at the highest LEP2 energy. On the other hand, according to 
the LAC1 parametrization the 2-res contribution should be dominant for y^i > 130 GeV, 
due to the rapid rise of the cross-section for gq — > jq, which profits from the steep gluon 
density assumed in this parametrization. In contrast, the GRV prediction for the 2-res 
contribution at y^i = 200 GeV is about a factor of 1.6 below the WHIT1 result, largely due 
to the smaller value of Aq CD that has been assumed for GRV; note that we are probing QCD 
at rather small momentum scales, where a s depends quite sensitively on Aqcd- 

The similarity of the various predictions shown in Fig. 29 is therefore somewhat mis- 
leading; clearly one could gain more information if the 2-res and 1-res contributions can be 
separated experimentally, e.g. by remnant jet tagging. Due to the large charge factor, the 
1-res contribution mostly probes the u— quark density in the photon. Its measurement might 
therefore yield valuable information about the flavour structure of the photon at rather low 
momentum scales, where nonperturbative contributions to the photonic parton densities are 
still important. Whenever the 2-res contribution is sizable, it is dominated by gq scattering; 
it could therefore give us an additional handle on the gluon content of the photon. 

Unfortunately the cross-section is not large even if we allow jP T to be as small as 1.5 GeV; 
the NLO calculation now being performed |53| will show how well QCD perturbation theory 
converges at such low momentum scales.[| Furthermore, the Bjorken— x of the partons (or 
photon) "in" the electron can only be reconstructed if the rapidity of the jet balancing the 
photon is also known. However, requiring |^j e t| < 1 reduces the cross-section by another 
factor of 1.7 (2.0) at y/s = 30 (200) GeV; moreover, the reconstruction efficiency for such 
a soft jet might be quite low. On the positive side, Fig. 29 shows that each TRISTAN 
experiment should have seen about 100 77 events with an isolated hard central photon 
in the final state; at LEP2, several hundred such events should be accumulated by each 
experiment. They are well worth looking for. 



5) Beyond Perturbation Theory 



In this section we will address several topics that go beyond QCD perturbation theory. We 
will be rather brief here, due to lack of both space and expertise on our part. Nevertheless, 



*As discussed in Sec. 3c, care has to be taken when computing the rapidities in this case. This has not 



been described properly in ref. 1 202 , although the numerical results of this paper are correct. 

^When calculating the LAC1 prediction shown in Fig. 29 we had to "freeze" the parton densities, i.e. 
use min(4 GeV 2 , (p^) 2 ) as momentum scale, since this parametrization produces negative parton densities 
if used at scales below the input scale Qq = 4 GeV 2 . 
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we feel that we should at least comment on these issues. For one thing, when comparing per- 
turbative QCD calculations with data, one nearly always has to correct for nonperturbative 
effects, usually with the aid of Monte Carlo event generators; we have already commented 
on the MC dependence of certain results in the preceding sections. Secondly, as remarked 
in the Introduction, the study of reactions involving real photons might provide us with 
new insight into the interplay of soft and hard QCD; given that perturbative QCD is by 
now generally accepted as the theory of hard hadronic interactions, much future research in 
strong interactions will presumably be focussed on this transition region. 

nonperturbative effects clearly play an important role in so-called minimum bias events. 
This name derives from the (idealized) definition that in a minimum bias event sample, each 
inelastic scattering event should be included, without any experimental (trigger) bias. This 
is very hard to achieve in practice, in particular for final states with low particle multiplicity 
and small energy deposition in the detector. Measuring total inelastic cross-sections is 
therefore not easy. 

One usually distinguishes three classes of events in jp scattering. In quasi-elastic events, 
7P — > Vp, the proton remains intact while the photon is transformed into a single vector 
meson V; usually only p, lu, <fr, and occasionally J/ip, are included here. In diffract ive events 
either the photon or the proton (or both) gets broken up into several hadrons, but no colour 
is exchanged between the two systems. In such events there is usually a gap in rapidity space 
between the two diffractive systems (in double diffractive events), or between the hadronic 
system and the single p or V (in single diffraction). In contrast, in non-diffr active events 
colour is thought to be exchanged between the photon and the proton. Both get broken up 
(in case of direct interactions, the photon is absorbed), and soft particles usually fill rapidity 
space more or less uniformly. 

All these components have to be modelled accurately for a complete understanding of 
minimum bias events, and in particular if one wants to extract the total 7p cross-section 
from the observed event rate. Both HERA experiments have performed P03 , 204 ] such 
measurements. HI | |204| | has so far only published an analysis of 1992 data, while ZEUS 
has also published ||203|| results based on the much larger 1993 sample. Both experiments 
use different MC models to describe the various event classes, and then try to fit the rel- 
ative weights to the data. For example, HI determines that about 26% of all events are 
quasi-elastic or diffractive; their result for the total jp cross-section at ^Js = 195 GeV isQ 
<7^ t (Hl) = (171 ± 7 ± 22) fib, where statistical and systematic errors are listed separately. 
The 1993 ZEUS data at yfs ~ 180 GeV favour a slightly larger quasi-elastic + diffractive 
event fraction (~ 35%), but a lower total cross-section: ct^ZEUS) = (143 ± 4 ± 17) pb. 
Very recently, HI announced ||205|| preliminary results from a special 1994 run with a min- 
imum bias trigger. Their analysis indicates an even larger quasi-elastic + diffractive event 
fraction (around 40%), and <r™(Hl, prelim.) = (172 ± 3 ± 10) ph. 

Note that the systematical uncertainties of these measurements are significantly larger 
than the statistical errors. This indicates that our understanding of minimum bias photo- 
production events still needs to be improved. Note also that the best description of non- 
diffractive events found by ZEUS [ p03|| , based on a superposition of "soft" and "semi-hard" 
events (see below), still has x 2 /d.o.f. > 2. The detection efficiencies, and hence the extracted 

*We have increased the published number by 8%, since HI used an inaccurate expression for / 7 | e when 
converting ep into 7p cross-sections 
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value of 0"*°*, determined from the same MC program may therefore also not yet be entirely 
reliable. 

So far little effort has been made to understand diffractive events in the framework of 



QCD. A possible connection ||206|| might be derived from the observation of jets in single- 



diffractive pp events by the UA8 collaboration ||207|| , as well as the recent observation of 
di-jet events with rapidity gap between the jets by the DO collaboration [ |208| |. Very recently 
the ZEUS collaboration announced [|209|| preliminary results indicating the presence of such 
di-jet events with gap between the jets at HERA; ZEUS interprets them as resolved photon 
events with large x 1 . 

Much more effort has been devoted to the understanding of non-diffractive inelastic 
reactions. This is necessary even for the study of hard processes, since the "underlying event" 
in events with high— jets is closely related to minimum bias events without any (obvious) 
hard interactions, nonperturbative effects again play an important role in the description of 
such events; however, it by now seems quite likely that semi-hard QCD interactions, with 
partonic transverse momenta of order 1 to 2 GeV, are also important here. Such interactions 
are said to lead to "minijets" : Even though a perturbative (hard) scattering took place, the 
resulting "jet" might be too soft to pass experimental jet identification cuts. 

The idea that such minijets might play an important role in minimum bias hadronic 
physics dates back to 1973, when Cline et al. proposed that the rapid increase of the 
inclusive jet cross-section with energy might be the root cause of the observed increase of 
total hadronic cross-sections. This jet cross-section is given by 



cr 



jet, 
ab \ 



j dp T f dx X f dx 2 ^2 fi\a(Xi)fj\b(x 2 ) daiJ -" kl ^ , 

Jip^/s JAp^/(x 1 s) i jkl dPT 



where subscripts a and b denote particles (7, p, . . .) and i, j, k, I are partons. s = X1X2S as 
usual, and a are hard partonic scattering cross-sections. Note that da/dpx oc p^ 3 ; the cross- 
section defined in eq.([JT]) therefore depends very sensitively on pT,mm, which is supposed to 
parametrize the transition from perturbative to nonperturbative QCD. 

If y/s ^> PT,min, the integral in eq.Q3"T|) receives its dominant contribution from xi j2 <C 1. 
The relevant parton densities can then be approximated by a simple power law, / oc x~ J . 
In case of pp or pp scattering, a = b and the cross-section asymptotically scales like [ p!0|| 




aJetoc ^lv— J log v— ' (32) 

if J > 1. For J ~ 1.3, as measured by HERA, the jet cross-section will therefore grow much 
faster than the total pp cross-section, which only grows oc log 2 s (Froissart bound ||211|| ), or, 
phenomenologically ||212|| for yfs < 2 TeV, oc s 08 . Eventually the jet cross-section ( |3~T1) will 



therefore exceed the total pp cross-section. 

This apparent paradox is solved by the observation that, by definition, inclusive cross- 
sections include a multiplicity factor. Since a hard partonic scattering always produces a 
pair of (mini-)jets, we can write 



0$ = (n jet pair Kf , (33) 
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where {n^ et pa ir) is the average number of (mini-)jet pairs per inelastic collision, a 
then implies (n- sct pa i r ) > 1, which means that, on average, each inelastic event contains more 
than one hard partonic scatter. The simplest possible assumption about these multiple 
partonic interactions is that they occur completely independently of each other, in which 
case ^jet 

p a i r obeys a Poisson distribution. At a slightly higher level of sophistication, one 
assumes these interactions to be independent only at fixed impact parameter b; indeed, it 
seems natural to assume that events with small b usually have larger n- ]et pair . This leads 



to the eikonal formalism [213 1 , where one writes the cross-section as an integral over the 



two-dimensional impact parameter b: 
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(34) 



Here Pa b d is the probability that both initial particles are in a hadronic state when they 
interact (see below), and A ab describes the transverse overlap of the partons. In realistic 
T4H it is necessary to introduce both nonperturbative (soft) and hard contributions 



analyses 
to the eikonal Xab- 

where a^ h is given by eq. 
evident from the identity 



Xab{s)=a^(s) + a^(s), (35) 
The connection to the Poisson distribution then becomes 
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(36) 



where the dots stand for the second line in eq.(^). Each term in the sum corresponds to 
the cross-section for events with exactly n hard partonic scatters; the second line gives the 

cross-section for events without any hard scatter [probability = exp ^— Aab ptld ab 1 , for fixed 

impact parameter b], but with a soft interaction. Note that the inclusive jet cross-section 
(BTp can be recovered from eq. (BB1) by introducing a multiplicity factor n in the sum: 
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(37) 



since the function A ab describing the transverse overlap of the partons in the two projectiles 
is by definition normalized such that J d 2 bA ab (b) = 1. 

We have already emphasized that, for fixed impact parameter, the ansatz ( |3~3| ) assumes 
hard scatters to occur independently of each other. This cannot be strictly true, since 
energy conservation implies that the maximal Bjorken— x for the partons in the second hard 
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interaction is less than 1, but it might still be a good approximation at high energies. A 
second assumption is that the dependence of the parton densities on Bjorken— x and on the 
impact parameter b factorizes; only in this case can the dependence on b be parametrized by 
a single function A a f){b). At present we are not (yet?) able to derive these properties from 
first principles. The ansatz (]3lf) therefore goes beyond perturbation theory, even though it 
postulates that perturbative interactions play an important role in minimum bias events; it 
clearly needs to be checked against experimental data to assess its validity. 

It has been shown ||214|| that eq. (|3"4]) with P^ d = 1 can be brought into agreement with 



data on total pp cross-sections, with pT,min around 1.5 GeV. In these calculations A p p(b) is 
usually computed from the Fourier transform of the charge form factor of the proton; this 
amounts to the assumption that colour charges track the distribution of electric charge in 
the nucleon. Moreover, the soft cross-section is usually parametrized as 



off {*) 



i °\ 



(38) 



where o"o,i are constants. Including pT,min, one therefore has three new free paramters (beyond 
those describing high— px jet production) to fit total and inelastic cross-sections (the elastic 
cross-section can be calculated in this model using the optical theorem); the success of such 
a fit is not entirely trivial. In particular, the rise of the total cross-section appears very 
naturally here, although the rate of increase is determined by PT,mim an d is thus fitted rather 
than predicted. These fits do have some weaknesses, however. They basically ignore the 
existence of diffractive events; more exactly, if the diffractive cross-section is supposed to be 
part of <r soft , eq. (|38|) , then eq.([3(J) leads to the prediction that the fraction of diffractive events 
(which are part of the second term, without hard scatters) decreases quickly with energy. 
A remedy within the framework of the minijet model is possible ||215|| , but only at the cost 
of introducing additional free parameters. Secondly, existing analyses typically assume that 
parton densities only increase slowly with decreasing x, at least at scale jj 2 ~ Px m i n j HERA 
DIS data favour a much steeper behaviour [14|. It is not clear whether the faster increase of 
0-jet p rec li c ted by such parton densities can be compensated by re-fitting the free parameters 
of the model. We should also keep in mind that total cross-section data can just as well be 
fitted in a more conventional Pomeron picture [ 2T2" |. 

Fortunately there is more direct evidence that an ansatz like eq.([34j) can describe some 
features of hadronic interactions. Using eq.(^), and the standard machinery to describe 
parton — > hadron transitions, the idea of having multiple partonic interactions within a 
single pp reaction can be built into an event generator |216| , [2171 , 218|| . In ref. [|216|1 it was 
shown that many features of pp scattering as seen at the CERN SpS can be understood in 
such a picture. These authors did not try to fit total cross-section data; instead, PT,mm wa s 
determined from the measured multiplicity distribution. It is encouraging that this again 
leads to a value around 1.5 GeV. Multiple interactions then naturally lead to the observed 
broad (non-Poissonian) multiplicity distribution. The perhaps greatest success of this ansatz 
is that it reproduces the "jet pedestal" effect. It had been observed that events with hard jets 
(Et > 15 GeV or so) also contain more hadronic activity far away from the jet cores (in the 
"pedestal") than minimum bias events do. Eq.(|3~4|) predicts such a behaviour, since events 
with hard jets are likely to have small impact parameter b, which significantly enhances 
the probability of having additional semi-hard interactions (minijets) in the same events. 
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Further qualitative and quantitative successes of this ansatz are described in refs. [ 2TB , 2T7 . 
However, it should be noted that these analyses again assume a rather flat small— x behaviour 
of the parton distribution functions. 

The most direct confirmation of eq.(|3~4]) would obviously be the observation of events 
with (at least) two independent jet pairs due to multiple interactions. These jets should be 
pairwise back-to-back in the transverse plane; the transverse opening angle should have a 
flat distribution for these events. Both properties are quite distinct from QCD 2 — > 4 events 
(usually called "double bremsstrahlung" ) . The cross-section for the production of at least 

two independent jet pairs can be computed by inserting the multiplicity factor ^ ^ ^ 

the sum in eq . (|36[) , and summing all n > 2: 
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(39) 



Energy conservation will again necessitate a slight modification of this simple expression. 

The first experimental search for multiple interactions was performed by the AFS collabo- 
ration plj at the CERN ISR (pp collisions at = 63 GeV). They required £ E T (jet) > 28 
GeV, which implies J2i x i — 0.45. It is doubtful whether multiple interactions at these rather 
large Bjorken— x can be treated as independent; indeed, implementing energy conservation 
here reduces the expected event rate by about a factor of 4. Moreover, no complete cal- 
culation of the QCD 2 — > 4 background processes was available at the time, so that this 
background had to be modelled using a leading logarithmic approximation in a region of 
phase space where no large logarithms occur. The large signal reported by the AFS collab- 
oration, corresponding to <j e fj ~ 5 mb, therefore should be taken with a grain of salt. 



More recently the UA2 collaboration ||220|| found some indication for the existence of mul- 
tiple interactions, but due to its rather small statistical significance (~ 3 standard deviations) 
they prefer to only quote the lower bound <j e fj > 8.3 mb. Finally, the CDF collaboration 
[ [2 2 1[ | found a 2.7 o signal, corresponding to <j c g ~ 12 mb, in the range expected by model 
calculations. 



In the first theoretical studies |222| of minijet production in 7p collisions, eq.(|3"4|) with 
P^ d = 1 was assumed to hold also for the case of photoproduction. In this case eikonalization 
effects are very small even at HERA energies, and a rather rapid increase of the total cross- 
section was predicted. However, as first pointed out by Collins and Ladinsky ||, P^ d 
should be 0(a em ). This enhances the exponent in eq. (|3~4D by a factor of order l/a em , so that 
eikonalization does become relevant at HERA energies if pT,mm < 2.5 GeV. The necessity 
to have P^ d ~ 0(a cm ) can most easily be seen from eq.(^), once we recognize that is 
0(a em )] if P^ d = 1, the production of additional minijet pairs would therefore be suppressed 
by powers of a em , even though only strong interactions are involved once the photon has 
been transformed into a hadronic state. Notice also that cancels out in eq. (|37|) . 

In ref.||, P^ d was taken to be the 7 — > p transition probability ~ 1/300. In ref . |223 



it was suggested to instead estimate it as the momentum fraction carried by partons in the 
photon at scale ji 2 ~ Pr,minj this gives slightly larger values P^ d ~ 1/200. One can (roughly) 
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reproduce the HERA measurements with either number. 

One might argue that P^ d should really be of order a em /a s ; after all, in eq. d36|) one wants 



the production of a second jet pair to be suppressed by a factor a 2 , not just a single power of 
a s J^\ However, it should be recognized that P^ d cannot be discussed independently of A lp {b). 



In refs.|223, 224 , 225 ] A 1P {b) was computed from the Fourier transform of the pion form 
factor; it is not at all clear whether this describes the transverse distribution of partons in a 
vector meson properly, and it is certainly not applicable to the perturbative ( "anomalous" ) 
component of the photon structure functions. Mathematically, eq.(|54]) only depends on the 



combination A ab /P^ b d ] one can see this by using the substitution b' = b ■ J Pab d '- 



a 



Z\s) = / d 2 b' [1 - exp (-A' ab (b') Xab (s))} , (40) 



with A' ab (b') = A ab [b'/^P^j /P^ d . Note that Jd 2 A(b) = 1 implies Jd 2 b'A'(b') = 1 as 

well; nevertheless b' cannot be interpreted as the physical impact parameter. Eq.(f4~0"D shows 
that one can always compensate an increase in P^ b d by a steeper decrease of A ab {b) at large 
b, which implies an increase of A(b) at small b. One can therefore simply fix P~~ to some 
value and fit A lp {b) to data, or vice versa; this approach was followed in ref. [ p25|| . 

However, in case of photons the ansatz ( p4|) might in any case be too simple. As already 
discussed in the Introduction, the hadronic structure of the photon is generally assumed to 
have both perturbative and nonperturbative components. The corresponding contributions 
to the parton densities have quite different x and Q 2 dependence. It therefore makes sense 
to split the cross-section ( |54"D into (at least) two terms, characterized by different values of 
P had , different parton densities [and hence different a^ et , see eq.(|31|)], and presumably also 
different overlap functions A(b). Yet another term has to be introduced to describe direct 
interactions; these cannot be eikonalized, since here the photon is "used up" after the first 
interaction. 



Two different ansatze of this kind have been proposed so far. Both Honjo et al. [226 



and Schuler and Sjostrand ||227| , |38f split the hadronic photon into a discrete sum over vector 



mesons, and a perturbative ( "anomalous" ) contribution; symbolically 

|7) = l7>bare + £ ^-\V) + JP^q^ (41) 



fv 



where |7)bare is the "direct" photonF], e/ fv are the 7 — > V transition amplitudes, \V) is a 
vector meson state, and \qq) is a state that develops from a hard 7 — > qq splitting. The 
sum in eq. flHD is assumed to be incoherent, so that each term also contributes separately to 
the parton densities and to the total 7^ cross-section. In refs. [|226|, |227|| the vector meson 



contributions to the parton content of the photon were described in terms of pion structure 
functions, while in ref. ||38|| the shapes of these contributions were fitted from K7 data (SaS 



parametrization; see Sec. 2b). The perturbative contribution to eq. fl4l|) is really a continuum 



^Recall that even the resolved photon contribution to is 0(a cnl a s ), since the parton densities in the 
photon are 0(a cm /a s ). 

■''Strictly speaking the coefficient of the first term in eq. (|4T|) should differ from unity by an amount of 
order a om ; this effect can safely be ignored. 
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of states; its contribution to the parton densities can be written as [3~8~| 



/■pert 
Ji\-y 1 



(42) 



Note that the perturbative 7 — > qq transition has been factored out here; it gives rise to 
the factor dk 2 /k 2 . The "state distributions" fi\ q q therefore obey homogeneous evolution 
equations in the factorizatioon scale /i 2 , with the (leading order) boundary conditions f38 
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fq\qq(x, k ,k ) — fq\qq(x, k ,k ) — — 
fq'\qq( X i k , k ) = fc\qq(x, k , k ) = 0. 
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(43a) 
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At this point the treatment of refs. ||226| and ||227| , p8[ diverges. Schuler and Sjostrand do 
not attempt to predict the total ^p cross-section. Rather, the emphasis of their work is on the 
properties of photoproduction events, both minimum bias events and events containing hard 
jets. They therefore assume a parametrization of the total cross-section as in ref. [212] . The 
overall normalization of the nonperturbative contributions to the photon structure functions, 
given by the fy, is fixed by low energy data. This also determines the contribution of this 
component to the total cross-section, since the Vp cross-sections are assumed to be identical 
to the 7ip cross-section. Moreover, they assume that the perturbative component (f42|) only 

contributes to hard scattering events, which are very rare at Js(jp) 



10 GeV. However, 

their choice of the fy gives a contribution to al^ which falls about 20% below the data at 
these low energies. The difference then has to come from direct interactions. This forces 
them to include direct jp scattering with partonic transverse momentum as low as 0.5 GeV.| 
The same value is also used for the cut-off parameter fco in eq . fl42|) ; this is in accordance 
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with ref. [£26] , as well as (approximately) with the GRV [ 47 1 and AGF [|55J parametrization. 

Ref. |[227|| devotes much attention to the description of quasi-elastic and diffractive events, 
which are entirely due to contributions from the \V) states. These states also contribute to 
soft and semi-hard (minijet) non-diffractive 7£> interactions, where the cut-off pr,mm — 1-3 
GeV at -y/i = 200 GeV is fixed from multiplicity distributions of pp events at the same 
energy. Schuler and Sjostrand allow for multiple interactions in this sector, but assume 
them to be completely independent of each other; in contrast, the ansatz (|3~4D assumes 
independent scattering only at fixed impact parameter b. We have seen above that the 
correlation between the presence of hard jets and small b offered a natural explanation of the 
jet pedestal effect. Assuming completely independent interactions means that every partonic 
interaction occurs with probability 



PjetiPl 



d°vp 



a 



nd 

Vp 



dp T 



(44) 



where the superscript "nd" stands for non-diffractive. Finally, in this model no multiple 
interactions are allowed to originate from the perturbative contribution (f4"2"D to the photon 



§In ref. p27| Schuler and Sjostrand modify the proton structure functions for scales fj? < 5 GeV 2 and/or 
small x, in order to enforce a smooth transition between DIS and photoproduction. However, in ref. (38) they 
use standard leading twist QCD to describe for momentum scales down to about 1 GeV 2 . 
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structure functions. In order to maintain the assumed s 08 behaviour of the total 'jp cross- 
section it then becomes necessary to increase pr,mm (in this sector only!) linearly with ^fs, 



so that p T ,min ^ 2.2 GeV for HERA (■ s Js(jp) ~ 200 GeV).Q 

The great advantage of the model of refs. | |227| , is that is comes pre-packaged in 
the successful LUND/PYTHIA Monte Carlo event generator ||228|| . This allows to test many 
aspects of the model against data from HERA, as well as from experiments at lower energies. 
On the positive side, not only does the total 7p cross-section seem to follow the "universal" 
s o.o8 behaviour, but the relative sizes of the contributions from quasi-elastic, diffractive, and 
non-diffractive events also match more or less the model predictions. The recent observation 
by the ZEUS collaboration [79 that the transverse momentum kx of the photon remnant 
jet in events with at least two high— p? jets is much harder than the Gaussian with width 
~ 0.4 GeV characteristic for hadronic remnant jets can also be regarded as a success of this 
model; indeed, it has been predicted more than ten years ago [ 223 ] that the perturbative 
contribution should have a kx distribution oc dk\jk\ for large kx-, which is directly related 
to the ansatz (|42"1). Moreover, the observation |)7| that in the region below about 1 GeV, 
the shape of the pt— distribution of charged particles at HERA closely matches that of pp 
events at the same energy indicates the presence of a component in the hadronic photon 
that behaves similarly to other hadrons. At higher p?, the spectrum becomes harder, in 
agreement with perturbative calculations |98| that include direct contributions as well as 
contributions from the perturbative part of the photon structure. 

There also appear to be some problems with this model in its present form, however. We 
already mentioned that the properties of the ZEUS minimum bias sample [ 2"U3 ] cannot be 
described properly by PYTHIA. Moreover, when trying to extract partonic cross-sections 
(and, ultimately, G 1 ) from the measured di-jet cross-section, the HI collaboration observed 
[ [70|| that switching on multiple partonic interactions in PYTHIA is not sufficient to describe 
the energy flow in resolved photon events|] It is at present not clear whether this problem 
can be solved by replacing the simple ansatz ( f4"4| ) by something like eq.(|34]) (for the contri- 
bution from the \V) states), or whether a more substantial modification of the model will 
be necessary. 

Finally, we find the introduction of three different scales that are meant to separate 
perturbative and nonperturbative interactions not very appealing. One might argue that 
PT.min ought to be larger for the pertubative \qq) component than for the \V) components, 
since the state that develops from a hard 7 — > qq splitting may have a smaller transverse 



size (as we will see shortly, this is assumed in ref.[ [226|| ), so that a larger momentum transfer 
becomes necessary to resolve individual colour charges. (A similar argument is used in 

*In an alternative version of this model, j>T,min is held fixed in the perturbative sector, while the /y are 
assumed to decrease with energy. Clearly this can at best be a temporary solution, since eventually the 
minijet contribution from the perturbative sector alone will exceed the assumed total 7^ cross-section unless 
it is unitarized in some way. Schuler and Sjostrand therefore favour the variant of their model with constant 
fv- They also discuss a few other versions, not of all of which are meant to be potentially realistic. 

t Since including multiple interactions improved agreement with the data. HI implicitly assumed that the 
remaining difference is due to even more partonic interactions; since these produce a jet pedestal that is 
uniform in (f>, they can then simply subtract this pedestal from the measured of the jet to arrive at the 
partonic (up to showering and fragmentation effects, which are presumably described adequately by the 
model). However, this treatment would not give the correct answer if the additional energy flow was (partly) 
due to, e.g., enhanced initial state radiation, which is not uniform in </>. 
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ref. [ 225 1 for the photon as a whole.) However, this does not explain why ko, which is also 
the minimal px of direct interactions, is so much smaller than even the value of pr.mm used 
for the semi-hard interactions of the \V) components. Alternatively, one might argue that 
quark exchange (which dominates direct interactions and, obviously, 7 — > qq splitting) is 
perturbative down to significantly smaller momentum scales than gluon exchange (which 
dominates minijet production in all resolved photon interactions). This could explain why 
k is smaller than pr.min, but the difference between the PT,mm values used for the \V) and 
I qq) components cannot be explained by such an argument. The different energy dependence 
of the two pr,min values is also difficult to understand intuitively. 

In this respect the model of ref . ||226 1 seems somewhat more appealing. Here a similar 
value for the cut-off k in eq.( [42]) is used, but all minjet production, from direct photons 
as well as the different classes of resolved photons, is assumed to be regularized by a single, 
energy independent parameter pT, m in- A fit to low energy data, which show a small but 
significant increase of the total 7^ cross section between 10 and about 18 GeV (the highest 
energy that had been reached in fixed target experiments), gives pr,min — 15 GeV. Since this 
is significantly larger than fc , Honjo et al. also allow the \qq) states to have nonperturbative 
interactions, which are assumed to scale oc 1/k 2 ; the idea here is that the "hardness" k 2 
of the 7 —>■ qq splitting determines the transverse size of the state that develops from this 
splitting. Minijet production from this state is also eikonalized, where the typical transverse 
size [which determines A(b)\ again scales like 1/k 2 ; that is, eikonalization a la eq.(|34D occurs 
for each value of k independently. As written in ref. [[2261 the model predicts a total 7p cross- 
section at ^/i = 200 GeV between 190 and 250 fib, the lower end of which is compatible at 
least with the HI measurement |[204|| ; presumably the agreement could be improved if one 
allowed the x and/or b dependence of the partons in the \V) states to differ from those in 
the pion. However, given the assumptions involved, we do not think that total cross-section 
data can be used to constrain parton distribution functions. Note that this model does 
predict a jet pedestal effect; it might even be stonger than in pp collisions, since for the 
perturbative contribution to the photon structure functions a very hard interaction not only 
implies small b, but also (at least on average) a larger value of k 2 and hence a narrower A(b), 
which increases the chances for additional semi-hard interactions. Unfortunately this model 
has not yet been built into an event generator. 

We finally mention two slightly different approaches to 7p scattering. In ref. [ 230 it was 
observed that, at least in a single-component model of the photonQ, where P^ d and A^ p (b) 
are determined from "canonical" VDM ideas, it is difficult to simultaneously describe the 
rise of 0"*°* between yfs = 10 and 18 GeV, and the rather small value of the total cross- 
section measured at HERA (or at least the smaller ZEUS result [[203|| ). These authors 
therefore allowed the soft contribution to the cross-section to grow oc s°- 058 ; similar to the 



behaviour found earlier in ref. ||232|| . Adding a minijet contribution estimated using the DG 
parametrization with pT,mm — 3.0 GeV (where eikonalization effects are still quite small for 
HERA energies) then gives a^p(y/s = 200 GeV) ~ 160 fib. A similar approach is taken in 



single component description was used here since in ref. [ |231| the perturbative component of the photon 
structure functions had been found to be too small to have much impact on minijet prod uctio n. However, 
comparison with the explicit calculation of ref. (3^] shows that the simple estimate of ref. [231], which used 
the double scaling limit of QCD, greatly under-estimates the perturbative contribution to G 7 at small x 
and scale /x 2 ~ p\ n 



,mm ' 
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ref.| |233| |. This work attempts a comprehensive description of ^p scattering, including quasi- 
elastic and diffractive events, in the framework of the Dual Parton Model | 234 |; it is now 
being implemented in the PHOJET event generator. Since the rise of the soft cross-section 
is now also fitted from data, in these models one can no longer claim to explain the rise of 
total cross-sections in terms of minijets. Moreover, pT,mm — 3.0 GeV seems rather large to 
us; there is evidence from J/ifj decays, DIS, and e + e~ annihilation that perturbative QCD 
is applicable at lower momentum scales. This is also supported by analyses of 77 data, to 
which we turn next. 

Once the description of the photon has been fixed, e.g. by fitting the free parameters of 



models like those of refs. [|226|, |227|, |233| from jp data, the properties of 77 interactions can in 
principle be predicted unambiguously ||235|| . However, two-photon experiments have needed 
comprehensive event generators well before the recent sophisticated models of 7p scattering 
were developed. We already discussed in Sec. 2b that the determination of Bjorken— x in 
deep-inelastic e7 scattering necessitates a model of the hadronic final state; in Sec. 4 we 
mentioned that the event generators that were used to model quasi-real 77 scattering prior 
to the pioneering AMY study || did not include resolved photon processes at all. 

The situation has certainly improved a great deal since then; however, there are still 
some problems. At present most experiments use the same basic ansatz to describe real 
77 scattering, which contains three separate contributions. Soft interactions, characterized 
by an exponential pr spectrum, are modelled using VDM ideas.0 The second component is 
estimated from the QPM. At high px it coincides with the direct contribution introduced in 
Sec. 4, but the px — > divergence is here regularized by constituent quark masses (typically 
m u = m<i ~ 300 MeV, m s ~ 500 MeV), rather than by a momentum cut-off. This QPM 
contribution therefore extends to (partonic) p? = 0, and is hence not entirely perturbative.^] 



In ref . ||236|| it has been shown that data on total 77 cross-sections and F2 can be fitted from 
these two components only. However, the description of multi-hadron and jet production in 
77 collisions is only possible if one also introduces single and double resolved contributions, 
which form the third component of current 77 event generators. 

So far all currently used generators agree. There are some differences in the details, 
however. For example, the TRISTAN experiments |6], |167| , |168|| , as well as DELPHI [[1771 , |178 
use parameter values for the description of the soft ("VDM") component that have been 
determined by experiments at lower energies. This might be dangerous, since these earlier 
analyses did not allow for resolved photon interactions; there is a strong correlation between 
the assumed size of the soft component and the parameter pT,mm, which essentially fixes the 
size of the resolved photon contribution to the minimum bias event sample. The ALEPH 
collaboration [|176|| therefore preferred to fit both p^min and the parameters describing the 



soft contribution from their own data. This might explain the rather large value of PT,mm, 2.5 
GeV, obtained by ALEPH for the DG parametrization, compared to 1.45 GeV for DELPHI 
and 1.6 to 2.0 GeV for AMY and TOPAZ. The expression for the total 77 cross-section 



§This "VDM component" is not to be confused with the nonperturbative contribution to photon structure 
functions, which is often also estimated from the VDM. As emphasized in refs. [227, 235 , these "VDM 
partons" also participate in hard resolved photon interactions. 



^Recall that in the model of refs. [227, 235 1 direct interactions with partonic pt as low as 500 MeV are 



allowed; in practice this should give a similar contribution as the QPM prediction using constituent quark 
masses. 
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derived by ALEPH also differs significantly from conventional VDM expectations. It should 
be noted that the jet cross-sections derived by TRISTAN experiments ||167| , |168|| are based 
only on high— px events; these analyses are therefore not sensitive to the assumed value of 
PT.min- Furthermore, the experimental definition (trigger conditions) of the "minimum bias" 
sample differs quite significantly between the various experiments, with ALEPH having the 
loosest cuts and largest visible cross-section, while AMY uses quite stringent requirements. 
Of course, a complete event generator should still give the same values (within errors) for its 
free parameters even when quite different data samples are used for the fit. The discrepancies 
between the values of PT,mm determined by different experiments therefore indicates that 
(some) present generators are not complete. 

One problem of these generators is that they do not include initial and final state radiation 
(parton showers). ALEPH states fL76| that generators like PYTHIA |228| and HERWIG 
|[237|| with parton showering switched on did not describe the data well. However, QCD 
tells us that such radiation must exist at some level. It might, for example, change the 



Pt~ dependence of the partonic cross-sections extracted by TOPAZ 



and AMY 



although the string fragmentation scheme used in these analyses may mimick some of these 
effects; a preliminary TOPAZ study [Tffl] finds that their MC program gives a reasonable, 
but not perfect description of the distribution of the transverse opening angle <fr between the 
jets in di-jet events. Note that initial state radiation should be treated differently for photons 
than for real hadrons ||227| , [238|1 , due to the presence of the hard ^qq vertex. Whenever one 
reaches this vertex in the standard backward evolution [|239|1 of initial state radiation, the 
shower has to be terminated, even if the parton's virtuality is still quite high at this point. 
This feature is included in the latest versions of PYTHIA and HERWIG, at least in an 
average {x— independent) sense. 

Another shortcoming of present 77 event generators is that they do not allow for the 
hard transverse momentum distribution of the remnant jets that has been predicted by 

This seems to have little impact on the 



ill- 



theory ||229|| and seen experimentally by ZEUS 
Pt— spectrum of the hard jets J79]] , but might change the predicted efficiency for detecting 
remnant jets, e.g. in the TOPAZ detector [[240| , |167| , |191|| . 



Finally, these event generators at present do not allow for multiple partonic interactions. 
Even if only the contribution from the nonperturbative part of photon structure functions is 
eikonalized, multiple interaction effects will play a significant role from LEP energies onward 
H|. In a single-component model of the photon, eq.(^|) with P^ d ~ a^ m , four jet events 



due to multiple parton scattering might already be detectable at TRISTANf], and ought to 
become a prominent feature of two-photon events at LEP2 energies ||241| |. The study of 
77 collisions at the highest available energies should therefore tell us how to unitarize the 
contribution to minijet production from the perturbative ("anomalous") component of the 
photon structure. 

Apart from its intrinsic interest, a good understanding of minimum bias 77 events is 
necessary for a realistic evaluation of possible hadronic backgrounds at future linear e + e~ 
colliders. In ref . ||242|| it was pointed out that in designs with strong beamstrahlung, several 
hadronic 77 collisions might occur in each bunch crossing, thereby effectively giving rise 
to an underlying event. This conclusion has been criticized [ [243|| on the grounds that the 

*Even in this case they are not expected to significantly affect the prop erties of events with one or two 
high— px jets; they do therefore not invalidate the analyses of refs.[167, 168 
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minijet contribution had not been eikonalized in ref . [|242| | . However, in refs. [fL80| , |244 



is 



has been shown that eikonalization does not significantly change predictions for the next 
generation of linear colliders, planned to operate at y/s ~ 500 GeV. The "hadron crisis" has 
not been solved by improved calculations, but by improved machine designs, which reduce 
beamstrahlung by employing a larger number of flat bunches. It should be emphasized that 
hadronic backgrounds are not the only, and perhaps not even the most severe, problem of 
designs with high beamstrahlung; they also suffer from large backgrounds from soft e + e~ 
pairs []245 |, as well as from a poorly defined beam energy, which complicates the study of 
thresholds [[2 46| . 

So far only one Monte Carlo study of the properties of these soft hadronic backgrounds has 
been published, by Chen et al. |[244|| . They assume that the total 77 cross-section tracks the 
measured pp cross-section at high energies. It should be pointed out that at present even the 
value of the total 77 cross-section at low energies is only poorly determined experimentally. 
Chen et al. then modelled minimum bias events using ISA JET [247] , which contains a 
parametrization of the underlying event as determined in pp and pp collisions. Using the 
same model for 77 collisions may not be a bad first approximation, given the similarities 
of some features of HERA and pp events discussed earlier. However, Chen et al. did not 
attempt to compare their model with actual data. They also used a rather large value for 
PT,min, 3.6 GeV, this being (approximately) the smallest partonic pr leading to observable jets 
at the UAl experiment [|248|| at the SpS. However, TRISTAN experiments have demonstrated 
|167| , |168| , |174|1 that much softer jets can be reconstructed in 77 collisions. The analysis by 
Chen et al., which uses the simplified unitarization scheme (f4"4"D, therefore probably under- 



estimates the effect of minijets on event characteristics.^] Chen et al. conclude that the 
average 77 event only deposits about 3.5 GeV of hadronic energy in the central detector 
(| cos 6*| < 0.9) at a typical design for a 500 GeV e + e~ collider. It would be interesting to 
repeat their analysis, using an event generator that has been tuned to describe real 77 data. 

Finally, we briefly mention that eikonalization effects [243] do certainly play an important 
role if the next e + e~ collider is converted into a 77 collider, which greatly increases the 
average 77 cms energy. The analysis of ref. [|244j| indicates that even here one can stay well 
below one hadronic event per bunch crossing, at least at a 500 GeV collider; the problem 
rapidly gets worse at higher energies, mostly due to the need to increase the luminosity oc s 
in order to achieve a constant rate of hard events. Since at those colliders 77 scattering 
occurs at energies well beyond the reach of current colliders, and in view of the incomplete 
description even of these accessible events, any prediction of event properties at future 77 
colliders should be taken with a grain of salt. 



6) Outlook 

In this article we have reviewed the present knowledge of resolved photon interactions. Great 
progress has been made in the last few years, both experimentally and theoretically. We now 

^As noted above, ISAJET simply fits the properties of the underlying event in pp collisions, including 
effects that can be explained in terms of minijets. Howe ver, i t is known that 72? events are more "jetty" 
than pp events, and hard partonic reactions are expected |23ql to play an even more prominent role in 77 
scattering. 
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know for a fact that partons "in" (quasi-)real photons play an important role in both ^p and 
77 interactions. Theory tells us that the hard ^qq vertex should lead to (a component of the) 
photonic parton densities that have a much harder x— dependence than the more familiar 
nucleonic parton distribution functions, and which should grow approximately linearly with 
the logarithm of the momentum scale at which the photon is being probed. Both properties 
have been confirmed experimentally, at least for the quark densities. Studies of photonic 
remnant jets also indicate that resolved photons do not always act like ordinary hadrons. 

At the same time many properties of minimum bias photoproduction and two photon 
events do resemble those of purely hadronic collisions at similar energies. Moreover, the total 
7P cross-section seems to show approximately the same energy dependence as total hadronic 
cross-sections. We therefore have good evidence for the existence of both perturbative and 
nonperturbative contributions to the hadronic structure of the photon. 

Qualitatively speaking this summarizes the present level of understanding of resolved 
photons. In the preceding sections we have reviewed the various pieces of information that 
contributed to this understanding; we have also described attempts to cast this understand- 
ing in a quantitative form. It seems rather futile to try and summarize this summary. We 
therefore decided to conclude this article with a list of open problems, which might become 
the foci of future work. 

The most pressing problems in deep-inelastic &y scattering (Sec. 2) have to do with the 
description of the final state: 

• The measurement of Bjorken— x has traditionally relied on a reconstruction of the 
invariant mass W of the hadronic final state from the measured quantity W v i S . An 
improvement of this procedure has recently been suggested [3J|. It might also be pos- 



sible to determine the energy of the target photon using the Jacquet-Blondel method, 
as is done by HERA experiments. Detailed Monte Carlo work, using realistic detector 
resolutions and acceptances, is needed to decide how well these ideas work in practice. 

• There are some discrepancies between existing data on F^- Improved reconstruction 
methods might be of some help here. It is not clear whether older PEP and PETRA 
data can be re-analyzed in this way, but new measurements and/or new analyses of 
F2 at small x 7 might help to resolve the discrepancy between recent results by the 
TOPAZ and OPAL collaborations. Such studies have the potential to discriminate 
between existing parametrizations of photonic parton densities. 

A list of open problems in hard ip scattering (Sec. 3) includes: 

• No complete NLO treatment of di-jet production exists. Recall that one needs to mea- 
sure the rapidities of both high— pt partons/jets in a hard event in order to reconstruct 
the Bjorken— x variables. 

• The measured jet cross-sections should be extended both in rapidity and in p?. The 
former increases the sensitivity to the interesting region of small x 7 , while the latter 
should allow to test theory cleanly, since a detailed understanding of the underlying 
event (see below) is less crucial at high p?, and differences between parametrizations 
of photonic parton densities are small at large x 1 and large momentum scale. 
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It might be interesting to try to correlate properties of the photonic remnant jet with 
those of the high— px jets. In the usual "nonperturbative + anomalous" description of 
the hadronic photon, the nonperturbative component should always have a remnant 
jet with very small kx', this component is also characterized by soft parton densities. 
In this picture one therefore expects nontrivial correlations between x 1 and kx- 

Studies of heavy flavour production hold great potential. We do not think it very 
interesting to try to derive total cross-sections from measurements covering only a 
limited region of phase space, which contains only a small fraction of all produced 
heavy quarks. It might be more fruitful to attempt to extract the resolved photon 
contribution, which is sensitive to the as yet poorly constrained gluon density in the 
photon. At high px, "excitation" contributions from the charm in the photon have 
to be taken into account. An important open problem is the fragmentation of rather 
soft (low— px) charm quarks, which contribute most to the total charm cross-section. 
Theoretical predictions are more reliable for b production, but it might be difficult to 
find a clean signal. 

The production of direct photons is by now quite well understood, although an NLO 
calculation of photon + jet production would certainly be welcome. Realistic back- 
ground studies are also needed, but can presumably only be performed by members of 
HERA experiments. 



In spite of the recent NLO calculation [13S] of direct J/ip production in the colour 



singlet model, much needs to be done here: The resolved photon contribution is only 
known to leading order in the colour singlet model. Nothing is known about the con- 
tribution from the colour octet component of the wave function of the J /if) |[249|| , which 



is accessible to resolved photons already in LO. In addition, there are contributions at 
high px coming from charm and gluon fragmentation ||250|| . 



• It is important to test our understanding of the hadronic photon in as many different 
channels as possible. The production of Drell-Yan lepton pairs, two photon final states, 
and associate J /if) + 7 final states are all plagued by rather small cross-sections, but 
this should at least partly be compensated by the cleanliness of the final states. 

We should emphasize that many of the measurements proposed here need more luminosity 
than HERA has delivered so far. Fortunately new data are being collected at an ever 
accelerating pace. 

Some open problems in hard 77 scattering (Sec. 4) are: 

• An NLO calculation of di-jet production is still lacking. 

• There is some discrepancy in the high— px end between the partonic cross-sections 
published by AMY and TOPAZ. Also, it is not clear how these cross-sections, which 
have been extracted with the aid of leading order event generators, can be compared 
to NLO calculations. 

• More detailed tests of theory should be possible if the kinematic reconstruction of 
events with high— px jets can be improved. One possibility might be to determine the 
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Bjorken— x variables using small angle electron taggers and/or the Jacquet-Blondel 
method. A study of the energy distribution of the remnant jet(s), and possible corre- 
lations with the hard jet(s), might also prove rewarding. All these measurements are 
probably quite difficult, however. 

• Measurements of the cross-sections for the production of high— p T particles should be 
relatively straightforward; they test similar aspects of the theory as studies of jet final 
states do. 

• The excess of high— p T charm events seen by TRISTAN experiments might well be 
explainable in terms of charm excitation from the photon; in the kinematical regime 
that is being probed by these experiments, this should be well described by the NLO 
event generator that has become available recently. Recall, however, that ALEPH does 
not see any excess. 

• J/ip and direct photon production has so far only been treated in LO. No data on 
these final states exists as yet. 

Finally, some of the most challenging problems go beyond perturbation theory (Sec. 5): 

• We have to understand the energy flow in resolved photoproduction events. This is 
crucial for the comparison of measured jet rates with theoretical (partonic) calculations, 
especially for presently accessible values of px- 

• The modelling of minimum bias photoproduction events needs to be improved. Among 
other things, this is necessary for an accurate determination of total ■yp cross-sections. 

• The total hadronic 77 cross-section is at present only determined very poorly. An 
accurate measurement can probably only be made if at least one of the outgoing 
electrons is tagged at a small angle. This is important for the assessment of hadronic 
backgrounds at future linear e + e~ colliders, which may play a role in the decision which 
of several competing designs should be pursued. 

• The role of multiple partonic interactions needs to be clarified. They will probably 
(help to) explain the energy flow in hard 7p events mentioned above. Theoretically 
least understood is the treatment of contributions from the perturbative component of 
photon structure functions. Studies of multi-jet production in high energy 77 collisions 
might yield important clues here. This also has ramifications for a theory (as opposed 
to parametrization) of total hadronic cross-sections. 

• 77 event generators in present use are known to be incomplete: They do not include 
parton showers, and assume too soft a k T distribution of the remnant jets. At higher 
energies multiple interactions might also affect the properties of the underlying event 
substantially. Finally, some standardization might help to decide whether putative 
differences in published results signal real discrepancies in the data, or are merely a 
reflection of different MC generators. 
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• It is important to test the (potentially) complete event generators for jp and 77 scat- 
tering that are now being developed in as many different ways as possible. Not only 
should things like the energy flow, multiplicity distributions, particle pt spectra, and 
strangeness production be investigated, there might also be interesting correlations 
between some of these quantities. 

This list is almost certainly not complete. Moreover, new questions are likely to arise 
as soon as current problems are solved. Much needs to be done before we can say with 
confidence that we understand the structure of light. 
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Figure Captions 



Fig. 1 Generic Feynman diagram for deep-inelastic e7 scattering: A probing photon with 
large virtuality Q 2 = — q 2 scatters off the (quasi-)real target photon (with virtuality 
P 2 ~ 0) to produce a hadronic final state X. 



Fig. 2 Feynman diagrams that contribute to 7*7 scattering, i.e. the blob in Fig. 1. (a) shows 
a tree-level (QPM) contribution, while (b) and (c) are examples of perturbative QCD 
corrections in the flavour non-singlet and flavour singlet sector, respectively, and (d) 
shows nonperturbative contributions as estimated in the VDM. All these contributions 
are formally summed in (e), where the quark densities in the photon q] are introduced. 



Fig. 3 Predictions of various existing parametrizations of photonic parton densities for the 
structure function K7 at Q 2 = 15 GeV 2 are compared with data from the OPAL |29 
(squares) and TOPAZ |33] (diamonds) collaborations. 



Fig. 4 Predictions of various existing parametrizations of photonic parton densities for the 
gluon content in the photon (a), and for the ratio of strange to non-strange light quark 
densities (b). Notice that we have used a logarithmic scale for the y— axis in (a), while 
in (b) x has been shown on a logarithmic scale. Both figures are for Q 2 = 15 GeV 2 . 



Fig. 5 Feynman diagrams contributing to the photoproduction of high— jets. Direct con- 
tributions from the QCD Compton scattering and photon-gluon fusion are shown in 
(a), while (b) gives examples of resolved photon contributions involving the scattering 
of two quarks, a quark and a gluon, or two gluons. 



Fig. 6 The topology of photoproduction events. In (a) the incident photon couples directly, 
leading to a final state with two high— p? jets (in leading order), and a remnant jet 
from the proton. In resolved photon contributions (b) the photon gives rise to a second 
remnant jet. 



Fig. 7 Example of a diagram that can be interpreted as either a leading order direct or resolved 
photon contribution, or an NLO direct contribution, depending on the virtualities of 
the exchanged quark and gluon; see the text for a detailed discussion. 



Fig. 8 The ratio of resolved photon and direct contributions to the inclusive jet cross-section 
at HERA, for jet pseudorapidity between -1.0 and +2.0. We have used the MRSD-' 
parametrization for the proton, and various parametrizations for the photon as indi- 
cated. These are leading order predictions, for momentum scale /1 2 = p\. While most 
curves have been computed using the simple antitagging requirement Q 2 < 4 GeV 2 , the 
dot-dashed curve has been calculated demanding Q 2 < 0.01 GeV 2 and scaled initial 
photon energy z between 0.20 and 0.75. 
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Fig. 9 Resolved photon (a) and direct (b) contributions to the inclusive jet cross-section at 
HERA. Contributions from qq, gq and gg final states are shown separately, where q 
stands for a quark or antiquark of any flavour. We have used the same antitagging 
and jet rapidity cuts as in Fig. 8. 

Fig. 10 The rapidity distribution of photoproduced jets with transverse momentum exceed- 
ing 8 GeV at HERA. The dotted curves show the contribution from directly inter- 
acting photons only, while the other curves show the total prediction using various 
parametrizations of photonic parton densities. The requirements for the scaled initial 
photon energy z differ in (a) and (b); requiring larger z clearly increases the differences 
between the predictions. 

Fig. 11 The triple differential di-jet cross-section at HERA, for the case of equal jet pseudo- 
rapidities. The direct and three subclasses of resolved photon contributions are shown 
separately. In (a) no cut on z is applied, while in (b) only events with 0.3 < z < 0.8 
are accepted; in particular, the requirement z > 0.3 removes the direct contribution 
for rji = r/2 > 0.2. 

Fig. 12 The average energy of the photonic remnant jet for resolved photon events at HERA 
with px = 10 GeV, as a function of the common pseudorapidity of the high— px jets. 
No cuts on the incident photon energy have been applied. The contributions from 
qq, qg and gg final states are again shown separately. 

Fig. 13 The scale dependence of the inclusive jet cross section at HERA, for jets with trans- 
verse energy Et = 25 GeV and pseudorapidity = +1.5, as predicted using the GRV 
parametrization for the photon and MRSD0 for the proton. The dotted, dot-dashed 
and long dashed curves have been obtained by varying the photonic factorization scale 
M 1 only, while the short dashed and solid curves are for the case where both factor- 
ization scales and the renormalization scale are identical. Adapted from Bodeker et al. 

IE- 
Fig. 14 The transverse momentum spectrum of cc and bb pairs, for the case where both heavy 
quarks are produced centrally (rapidity y — 0). Note that at high p T the single charm 
(or bottom) inclusive cross-section will receive substantial "excitation" contributions 
involving the heavy flavour content of the photon, which are not included here. Even 
though we have used a parametrization with a hard G 1 , the resolved photon contribu- 
tion is clearly very small. 

Fig. 15 The triple differential cross-section for the photoproduction of cc pairs at HERA, for 
Pt{c) = 10 GeV, as a function of the rapidities of the two heavy quarks, which are 
taken to be equal. These are LO results; MRSD-' has been used for the nucleonic 
parton densities. For the WHIT1 parametrization contributions from gg fusion and 
qq annihilation are shown separately, while only the sum is shown for the WHIT4 and 
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LAC1 parametrizations. Note that for the range of rapidities shown here, the cut on 
the scaled incident photon energy z removes the direct contribution completely (in 
leading order). 



Fig. 16 



Fig. 17 



Fig. 18 



Fig. 19 



Fig. 20 



Fig. 21 



The rapidity dependence of the direct photon cross-section for incident photon energy 
E 1 = 10 GeV and transverse momentum of the final photon pj, = 5 GeV, as predicted 
from a leading order calculation [ |12C|| using the ABFOW parametrization ||121|| for 
the proton and AGF for the photon. Note that the rapidity r? 7 is measured in the 
7P centre-of-mass frame here. The contribution from the gluon in the photon clearly 
dominates for rj 7 > 1. 



The rapidity dependence of the direct photon cross-section at HERA for pj. = 5 GeV, 
as predicted from an NLO calculation ||123|| using the GRV parametrization 



50 for 



the proton. Note that rj 1 is measured in the HERA lab frame. No cut on the energy of 
the incident photon has been applied; the contribution from the gluon in the photon is 
sub-dominant everywhere. Some sensitivity to the quark densities in the photon still 
remains, as demonstrated by the difference between the solid (GRV) and dotted (GS) 
curves. 



A Feynman diagram contributing to the photoproduction of J/ if) mesons in the colour 
singlet model. For the resolved photon contribution, the incident photon has to be 
replaced by another gluon. 



The J/i/j production cross-section at HERA is shown as a function of the elasticity 
parameter Z of eq. (|23"|) , as predicted from a leading order calculation [|146|| . Acceptance 
cuts have been applied on the transverse momentum of the J/ ip meson and on the angle 
of both leptons that originate from J/ip decay, but the leptonic branching ratio has 
not been included. 



Typical Feynman diagrams (left) and topologies (right) of the three classes of contribu- 
tions to the production of high— px jets in (quasi-)real 77 collisions: direct (a), single 
resolved (b), and double resolved (c). Note that each resolved photon gives rise to a 
remnant jet. 



Leading order predictions for the inclusive jet cross-section at TRISTAN. The lower 
dotted curve in (a) shows the direct contribution only, while the other curves show 
the total prediction using various parametrizations of photonic parton densities, as 
indicated. The cuts on the pseudorapidity of the jet and on the maximal scattering 
angle of the outgoing e ± differ, (a) corresponding to the cuts used by TOPAZ and (b) 
to those used by AMY; note that the antitagging cut is only effective if z < 0.75. Data 
by these two experiments [ |167| , |168| are also shown. 
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Fig. 22 The triple differential di-jet cross-section at TRISTAN, as predicted in leading order, 
(a) shows the direct and 1-res contributions, while (b) shows 2-res contributions. For 
the WHIT1 parametrization the contributions from different final states are shown 
separately, while only the sum is given for the DG and WHIT6 parametrizations. The 
results are for TOPAZ antitagging conditions. 

Fig. 23 The triple differential di-jet cross-section at LEP1 and LEP2, as predicted in LO using 
the WHIT1 parametrization. The steps at pseudorapidities 771 = 772 — 1-5 (for \/s = 90 
GeV) or 2.2 (for yfs = 180 GeV) occur since the ALEPH antitagging condition has 
been used, which limits the scattering angle of the outgoing only if it carries more 
than 50% of the beam energy. 

Fig. 24 The transverse momentum distribution of jets produced in 77 events at two stages of 
the planned JLC collider, as predicted in LO using the WHIT1 parametrization. The 
dotted curves show the total result in the absence of beamstrahlung, while the dashed 
curves show different contributions as indicated, and the solid curves their sum, for 
realistic beamstrahlung spectra. The comparison between the solid and dotted curves 
show that the beamstrahlung spectrum is expected to be significantly harder at \fs = 1 
TeV (b) than at 0.5 TeV (a). 

Fig. 25 The triple differential di-jet cross-section at a 77 collider, based on a 500 GeV e + e~ 
collider with unpolarized beams. Due to the hardness of the photon spectrum, there 
are substantial differences between predictions using WHIT1 and LAC1 even at the 
high value of pt chosen here. 

Fig. 26 Predictions for the total cross-section for cc pair production in 77 collisions at e + e~ 
colliders. The direct contribution (dotted) includes the parametrized form (p8|) of NLO 
corrections, whereas the resolved photon contribution (dashed and solid curves) has 
been computed in LO only. Two curves of a given pattern show the uncertainty of the 
theoretical prediction, as described in the text. 

Fig. 27 The signal for a Standard Model Higgs boson of mass m# = 120 GeV decaying into a 
bb pair, as seen at a 77 collider based on a 350 GeV e + e~ collider, as well as various 
backgrounds. 90% electron beam polarization has been assumed in order to suppress 
the direct background. Both cc and bb pair production processes and c and b excitation 
processes are included in the background calculation. The width of the Higgs signal, 
and of the Z peak, is determined by the resolution, not the intrinsic widths of these 
particles. From Baillargeon et al. |[201|| . 



Fig. 28 The total cross-section for the production of J / ip mesons in 77 collisions, as predicted 
by the colour singlet model. The dotted curve shows a LO prediction for scale /i 2 = 
m j/V> w hile for the dashed and solid curves NLO corrections to the 1-res jg fusion 
contribution have been included using the parametrization of eqs. (|29|) , (|30|) , but the 
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NLO direct and 7g scattering contributions have been omitted. Two curves of a given 
pattern show the theoretical uncertainty of our estimate, as described in the text. 

29 Leading order predictions for the production of a hard, central photon in no-tag two- 
photon collisions (no antitag condition has been imposed). Fragmentation contribu- 
tions are not included. The long and short dashed curves show 2-res and 1-res contri- 
butions as predicted from the WHIT1 parametrization; they add up to the solid curve. 
For all other parametrizations, only the total sum has been shown. 
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